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ABSTRACT 


An interactive computerized experimental procedure is 


presented for obtaining the elastic-plastic fracture 
toughness (JIce) for a material. The process employs a bend 
type, Single specimen approach by using the measured 


unloading crack opening displacement (COD) compliance to 
obtain the crack length. Data are taken, manipulated and 
plotted during the test to provide the operator with the 
necessary information to control the testing machine. The 
program automatically plots the J versus crack extension 
curve and determines J[Ic¢ upon completion of the test. A 
second program isa presented to assist the operator in 
fatigue precracking the specimen. A third program is used 
to recall, list and plot the data that were stored on 


diskette during the test. 
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I. INTRODUCTION 


A. GENERAL 

The toughness of a material is a measure of the energy 
per unit volume required to cause fracture. It can be 
quantified by simply taking the area under the true stress 
strain curve from a aimple tension test (Ref. 11]. The 
practical usefulness of this measurement is limited however, 
because it provides no information on the material’s ability 
to resist crack propagation. Since all real atructures 
contain notches or cracks of some sort, by design or 
otherwise, which act asa asatress concentrators, another 
measurement of toughness was developed called fracture 
toughness. 

The fracture toughness (K.-) of a material is 4a saeasure 
of the resistance of the material to fracture in the 
presence of a defect. K ia defined aa the plane strain 
fracture toughness and is used to quantify a aaterial’s 
resiatance to crack propagation (Ref. 2]. The American 
Society for Testing and Materials (ASTM) developed 


Standard E399-83 (Ref. 3] which establishes precise 


details for accepted testing procedures. The standard 
also provides rigid criteria for the specimen size and 
geometry to ensure plane atrain conditions exist at the 


tip of the crack. 





The J-integral method, developed by Rice in 1968, 
provides a saeana of including plastic deformation at the 
crack tip (Ref. 4]. J ia a measure of the intensity of the 
entire elastic-plastic satresa-atrain field that surrounds 
the crack tip. JIe ia defined as the value of J at the 
onset of cracking. J-integral teating has gained acceptance 
as a valuable means of evaluating a metallic material’s 
resistance to crack initiation as deaonstrated by 
eee ahlietwent eof ASTM Standard £E813-81 (Ref. Si]. The 
neasured value of JIce has practical applications in 
measuring the effects of smetallurgical variables, heat 
treatments, and weldments and is being used by researchers 


today to make such quantitative distinctions (Ref. 5-7]. 


B. JIe VERSUS KIe TESTING 
The procedures eatablished by ASTM Standard E399-83 for 

obtaining Kie impose ainimum specimen size requirements 
that cause difficulty when teating materials that possess 
high toughness. ‘The specimen thickness must be sufficient 
to cause plane strain conditions to exist (Ref. 3]. 
Hertzberg (Ref. 2] presents the empirical equation 
developed by Brown and Srawley (Ref. 8]: 

t and a> 2.5 * (KIe¢/Ya) 
where 

t = sapecinen thickness 


@ = crack length 
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KIe = plane atrain fracture toughness 
Ya = yield atrength 
A AY wider analyaia for J-integral teat specimens reveala 
that: 
bor B > 25 # (JIe/Fa) 
where 
b = initial uncracked ligament 
B = specinen thickness 
JIe = anticipated value if JIc 
recheck when test is corupleted 


Fa = flow atreas ((Ys+Uta)/2) 


Hertzberg takea this reault one step further by showing that 
Jte equala (KIe)“=/E (where E equals the material’s elastic 
modulus) for Jie. By dividing the two equations one sees 
that: 
CKIe specimen size) >» 20 * (Jf[e specimen size) 
thus, because the J-integral allows plastic deformation, one 
of ita most aignificant advantages ia that a satisfactory 
epecimen can be obtained at one-twentieth (1/20) of the 
thickness of a Kj{ce specimen. This not only reana aaaller 
specimens ray be used but, also, smaller teating machinea 
are required and therefore coats are reduced. (Ref. 2] 
Although JIe and Kfe provide a different reassure of 


toughness, JI-« can be used to eatinmate KIc. ASTM Standard 
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E813-81 (Ref. SJ provides limited guidelines on the 


relationship between the two. 


C. PURPOSE AND SCOPE OF RESEARCH 

The purpose of this thesis is to develop the procedures 
required to obtain J[Ic from J-integral teat resulta using 
equipment currently available within the Materials Science 
group of the Mechanical Engineering Department at the Naval 


Postgraduate School (NPS). A list of the major equipment 


follows: 
a. Series 810 MTS rachine 
b. HP-9826 coauputer with internal disk drive 
and 320K random access menory 
c. HP-3497A data acquisition aystem 
d. HP-3437A voltmeter 
e. HP-2671G thermal printer 
£. HP-722SB graphics plotter 
Figure 1 shows the system layout. - 


The test method developed ia to comply with ASTM 
Standard E813-81 (Ref. 3), be easy for others to use, and 
require a minimum of operator interaction. it is 
specifically designed for use with three point bend 
specimens using a single vice rnultiple specinen approach. 
The method developed and presented here covers the entire 
teat, including specimen preparation, precracking, testing, 


and validation of results. 


le 





II. NATURE OF THE PROBLEM 


A. J-INTEGRAL 
The J-integral ia a mathematical expression used _ to 
characterize the near tip atress-atrain environment of 
cracked elastic-plastic bodies (Ref. 9]. It is a path 
independent line integral that surrounds the crack tip for 
non-linear plaatic metallic materials (Ref.5]. The equation 
delveloped by Rice is (Ref. 4]: 
J oh dy - T°(€au/dx)> da) 
where 
W = loading work per unit volume, or for elastic 
bodies, atrain energy 
f = the path of the integral, which encloses the 
crack tip 
da = increment of the contour path 
T = outward traction vector on ds 
u = diaplacement vector at ds 
X,yY,zZ = rectangular coordinates 
T:(@u/axdda = the rate of work input frorma the satressa 
field into the area enclosed by[ . 
Figure 2 providea a pictorial representation of this 
equation. 
The mathematical expression developed by Rice (Ref. 4] 


has been amanipulated and reduced to aimpler forma for 
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specific sapecimen geometries. Rice, Paris, and Merkel [Ref. 
10] reduced the J-integral expression for a three point bend 


specimen to: 


where _ 

J = resistance to crack initiation 

A = area under load versus load line 

displacement curve 

b = uncracked ligament 

B = specinen thickness 
This thesis ia not intended to redevelop the J-integral but 
rather apply it, therefore only the result is presented 
here. 

More recent work by Ernst, Paris, and Landes (Ref. 11] 
corrects the calculation of J due to crack growth by means 
of the following equation: 

Jcie1) = CJyrCn/B) 4 RAy,ie4) * C1-CY/D) # Cageyz-ayz)) 
where 

J(i+1) = current value of J 

J; = last J coaputed 

n = 1 for three point bend type specimens 

Aj, ie*l = area enclosed by the actual test 
record and lines of constant 
displacement dj and dj+ji 


Y = 1 for three point bend type specimens 
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b = current uncracked ligament 

W = specimen width 

aj = last crack length 

@i+1 = current crack length 
Correcting for crack growth during the teat is necessary 
when using a single specimen approach as used by this 
method. 

The blunting line defined aa: 


J =2*s* Fa * deltaia 


where 
Fa = flow atress 
delta. a = change in crack length 
approximates apparent crack advance due to crack tip 


blunting in the absence of slow-asatable crack tearing. It is 
plotted on the J versus crack extension plot. A linear 
regression fit through the applicable data points ia then 
plotted. The intersection of the two lines is by definition 


JIe ‘(see Figure 3).(Ref. 5] 


B. PREVIOUS WORK 

The early experimental methods used to obtain Jf{ec 
required the preparation of several apecinmens. Each 
specimen received a different load. The specimens were then 
heat tinted (to identify crack length) and broken apart. 
Each apeciaen provided one data point on the J veraus crack 


extension curve, commonly referred to as the J-R curve [Ref. 
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9 and i1lj. Checks were then performed on the validity of 
the data. A minimum of four valid data points were required 
to obtain a linear least squares fit and hence J[{c (Ref. 
Sprorand 12). 

Researchers over the last five years have developed a 
much improved test method whereby a single specimen can be 
used. Today’s computerized data acquisition equipment and 
data manipulation are ideally suited for the task. By 
relating the compliance Chow much a specimen bends under a 
given load) to the crack isngth, all the data required to 
obtain the J-R curve can be collected by using a asingle 
specimen [Ref.11-16]. Clarke {fRef.16] states that prior to 
the use of computers the time required for data reduction to 
obtain the J-R curve was four to five hours. Using 
computerized, single specimen, unloading compliance methods 
it is now possible to to obtain the J-R curve in a matter of 
twenty to thirty minutes after the test starts. The obvious 
savings in specimen preparation and data reduction have made 
the J-integral a valuable means of analysis where no 


practical meana previously existed. 


C. TEST CRITERIA 

The test method deveioped in this thesis ia specifically 
designed for use on athree point bend specimen. fe 
utilizes the single specimen approach by relating the crack 


opening displacement (COD) compliance to the crack length. 
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The following equations are currently used by DTNSRDC. For 
bend specimens with a span to width (S/W) ratio of 4.0, Dr. 
J. Joyce (Ref. 171) of the United States Naval Academy 
determined: 

a/W = 0.998265 - 3.81662 * U - 1.80596 * U2 + 32.31038 * US 


- 44.15665 « U4 - 52.67876 *# US 





where 
@ 2 crack length 
W = sapecimnen width 
af 
U = 
+> Pe Bao Wr sonen yes 
where 
Be = effective thickness = Bmax - (Bmax - Bain)2/Brax 


E = material’s elastic modulus 

S = span between supports 

W~ = specimen width 

compl 2? compliance = i/slope of the load versus COD line 
for the specific unloading. 

The reason for using the single specimen approach has 
already been presented. The reasons a bend type  specianaen 
and COD type unloading compliance are used are: 

a. Three point bend specimen georetry is better suited 
than coapact tension specimens for testing weldments. 
b. Three point bend specimen geometry is also easier to 


machine than compact-tension specimens. 
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COD voltage readings are less susceptible to noise 
than linear variable differential transformer (LVDT) 
voltages used to measure load line displacement. 

The COD compliance equation affords 0.1% fitting 
accuracy and allowa comparison of results with 


DTNSRDC test data (Ref. 17]. 
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III. TEST METHOD 


A. OVERVIEW 

A J-integral teat can be divided into three distinct 
ateps, 1) pre-test, 2) test, and 3) post-test. Each of 
these are described in greater detail in this chapter and 
presented as a satep by step procedure in Appendix A. 

After the asapecirmen is machined it must be  precracked 
using cyclic fatigue loading to minimize residual satresses 
and create a sharp crack tip (Ref. 5]. A computer program 
was developed to assist the operator in accomplishing this. 
The program computes the loading parameters for precracking 
in accordance with ASTM Standards E399-83 and ASTM £E813-81 
CRef. 3 and 5]. The operator uses this information to run 
the MTS rnachine. 

During the J-integral teat a second corputer program 
collectja and records the load, COD, and load point 
displacement (deflection). Theae data are used by the 
program to determine J and the amount of crack extension for 
each operator controlled unloading cycle. An unload cycle 
consiatsa of decreasing the load by 10x of the expected 
maximum load and reloading to the original load point. The 
prograrm provides the operator with a real-time graphical 
representation of the any two parameters (load versus COD or 


deflection) during the test. It will also plot the J-R 
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curve and compute JIc upon completion of taking all the test 
data. Data are stored on diskette to serve as a record and 
allow future recall listing or plotting. 

The post-teat process includes heat tinting the specimen 
to aid in visual verification of crack extension upon 
cormpletion of the test as required by ASTM Standard £813-81 


(Ref. SJ. 


B. SPECIMEN PREPARATION 

Figure 4 ahowsa apecifically the dimensions used in 
preparing eight specimena fror two different one inch thick 
walded plates of HY-80. One plate was preheated while the 
other was not. It was necessary to reduce the width to 
approximately 0.95 inches to remove the distortion Cbow) 
created by the weld. Four specimens were prepared from 
each plate, the machined notch was placed in the weld netal 
for two and in the base retal for the other two. The clip 
gage ‘(used to measure COD) ia attached to the top of the 
g@pecirmen by means of two cut razor blades. The bladea are 
held down by small straps which allow adjustment of the 
initial gap. The initial gap is set to 0.475 inches, the 
calibrated zero édteing for the clip gage, by means of a 
machined gage tool that rests in the asapecimen’s machined 
notch. This method ia also used by DTNSRDC. 

Two additional holes are drilled and tapped for 


retaining a samall bracket on the face of the specimen. The 
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bracket is used to provide a sensing point for a linear 
variable differential transformer (LVDT). The LVDT signal 
is used to calculate the area under the load versus 
diaplacement curve. At present the external LVDT is unable 
to control the MTS machine so the current method uses the 
MTS’ internal LVDT to measure deflection via actuator 
movement. 

Clarke (Ref. 16] states that side grooves improve the 
accuracy of calculating the crack length when using 
unloading compliance procedures. The side grooves reduce 
plastic deformation on the surface and aid in propagating a 
straighter crack front. DTNSRDC has’ found 20% side 
grooving produces the best results, hence the same was used 
for this work. 

After the specimen is prepared a test record sheet (see 
Appendix B) is completed. Notch and side groove depths can 
be measured using a dial indicator with a needle tip. The 
record provides the operator with ready access to the 
information necessary to run the precracking and test 


programs. 


C. PRECRACKING 

Before the J-integral test can be performed the 
specimen must be precracked to establish a sharp crack tip. 
ASTM Standards E813-81 and E399-83 establish parameters that 


must be met during precracking to prevent over-stressing the 
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apecimen or propagating the crack through the entire width. 
Precracking is accomplished with the aid of a prograa called 
“PRECRACK’ which will compute loading requirements based 
upon these paramenters. 

The program, written in HP-Basic, affords flexibilty by 
allowing the use of different bend specimen sizes’ and 
materials. It ia easy to use due to its interactive design 
and easy to modify should auch a need ariae. A liating of 
“PRECRACK”™ ia provided in Appendix D. 

The operator is requeated by the program to enter 
apecific apecimen data. The program will compute the initial 


maximum and ainimum loada and change in the aaaociated 


atressa intenaity factor (delta K). If this delta K exceeds 
40 kasisinl/2 then new rmaximum and sninimum loads are 
cornputed. The operator is responsible for ensuring the MTS 


machine is operating within the prescribed limits. The 
operator can obtain a hard copy of the entered and computed 
data for record purposes (see Appendix G). 

The MTS machine producea cyclic loading by means of an 
inverted (coaupreasive load) haver sine wave. 

The operator rust periodically interrupt the cyclic 
loading and measure the current crack length by amanually 
loading and unloading the apecimen (determine the specinmnen’s 
coapliance). The program coraputes the new crack length and 
new loading limita. During the last 50x of precracking: 


Pmax < .4PL 
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where 
Prmaax = maximum load 
Pr = (€(4/3) * (B © b2 * Fa / (2 * Ha))]) 
where 
B = net thickness (Bmin if specimen has side grooves) 
b = current uncracked ligament (W-a) 
Fa = flow stress {(CYa+Utsa)/2) 
Ha = half span 
and Delta_K/E < 0.005 manl/2 (0.001 inl/2) whichever, 
produces the least load (U[Ref. Sy) The prograr will 
automatically begin using these parameters when the an,ost 
recently calculated crack length exceeds this 30x of the 
required extension to reach a/W = 0.6. This value is used 
because it ia the optimum crack length with which to begin 


the J-integral teat (Ref. 5). 


DBD. J-INTEGRAL TEST 

The J-integral test is performed by slowly applying a 
load to the speciren and periodically unloading te by 10x of 
the expected maximum load. The MTS machine is set up in 
atroke control with a very srall inverted rarp function 
resulting in a loading rate of 0.005 inch/minute. The 
operator aust again interrupt the MTS program and manually 
unload the specimen. Each unloading cycle Cunload = and 
reload) provides the neceasary information to compute the 


current crack length. The crack length is then used _ to 





determine the value of J up to that point. The computer 
program “J_INTGL*" provides the operator with pertinent 
information during the test. Once the test is completed, 
the operator can obtain the J-R curve and value of JIce (see 
Figure 4 and Appendix H). The program will use only data 
points between the offsets, (blunting line + O.15mm) to 
(blunting line + 1.5mm) as required ASTM Standard £813-81 
when coaputing the linear regression line. A minimum of 
four data points must be within the offsets and one nust lie 
near (within 1/3 the overall range of delta_a) the blunting 
line. Specific acceptability criteria for data is provided 
in ASTM Standard E£813-81 (Ref. 5]. 

The progran “J_INTGL” is the foundation of the test 
method, the single specimen approach would not be possible 
without it. The program collects, manipulates, diaplays and 
records all pertinent data during the test. A flow chart of 
“J_INTGL” is provided in appendix C. A listing of “J_INTGL”™ 
ia provided in Appendix E. 

Deflection measurements are presently taken using the 
LVDT installed in the MTS _ smachine. A correction is 
therefore required to account for the system’s compliance 
while loading a specimen. The correction was obtained by 
loading the system with a two inch thick piece of HY130 and 
recording the stroke reading for various’ loads. This 


process was checked three timea with consistent (within 5x) 
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results. A plot of the data and the 4th order ploynomial 


curve fit is included as Figure 5. 


E. POST TEST VERIFICATION 

The results rust be verified upon completion of the 
test. ASTM Standard E813-81 presents the following criteria 
that muat be met before the test can be considered valid: 

1. Total crack extension (visually measured) must be 
within 15x of that recorded during the test. 

2. by and Bmin > 25 #* JIe / Fs (Call terms are as 
defined earlier except bj = post-fatigue uncracked 
ligament). 

3. The slope of the linear regression fit to the data 


must be less than flow stress (dJ/da < Fs). 
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IV. RESULTS 


A. TEST RESULTS 

Specimen HY80-5B, prepared from the base metal of the 
preheated welded plate, was precracked,and tested resulting 
in a Jtc value of 142 kJ/m@. Figure 8 is a tracing of a 
photograph of the final Pear eareraurtace after heat tinting. 
All post test checks were completed and yielded satisfactory 
results. The JIc value is in good agreement with previous 
tests of HY80 plate as conducted by DTNSRDC resulting in an 
average JIc of 144 kJ/m2 (Ref. 6]. 

Plots of J-R curve, load versus COD, load versus 
deflection and the are included as Figures 5, 6 and 7, 
respectively. The load versus deflection curve was obtained 
using the program "J_DATA”™, thus proving that the data is 
stored on disk, and is able to be recalled. 

Data was taken at a nearly constant rate of 2.5 data 
points per second (150 per minute) while the program was 
executing. 

Test reults were verified following the guidelines 
provided by ASTM Standard E813-81. The fracture surface was 
enlarged by taking a photograph, a tracing of which is 
included as Figure 8. The observed crack length was 14.4% 
longer than the last recorded crack length. This is within 


the acceptable limit (recall ASTM requires < 15%) but 
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because the clip gage reached its limit detector and caused 
the test MTS to stop between unloadings, it isa possible some 
crack extension occurred after the last unloading cycle. 
This would help to explain the difference. 

Additional information on this specific test is provided 
by being used in the example calculations shown in Appendix 
A. The exemplary computer output for each of the progrargs 


is also from specimen HY80O-SB. (See Appendixes G, H and I). 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The method developed and presented by this work produces 
accurate and reliable results that meet the specifications 
established by ASTM Standard £813-81. Precracking specimen 
HY80-SB took approximately 2 hours and the J-test required 
approximately 2.5 hours including several pauses to arganually 
verify several calculations. The time required for the 
processes will most likely decrease as operator experience 
and confidence increases. 

All three computer programs have run successfully and 


perform as intended. 


B. RECOMMENDATIONS 

As with all first versions of computer programs as 
extensive as “J_INTGL", changes and improvements are to be 
expected. The following programming changes ray be of value 
to future operators: 

1. Display pertinent information (i.e. current crack 
length in both English and metric units (currently 
only metric is provided). 

2. Provide a means of storing the data arrays during the 


test without being required to rerun the entire 





program (re-entering specimen data, conversion data, 
etc.) 

3. Provide a means of eliminating a specific J-R data 
point from being used in obtaining the linear 
regression fit. This may be programmed to be 
accomplished automatically by applying Chauvenet’s 
Criterion to the data. 

Improvements in the current means of conducting the test 
can also be made, such as making the MTS control both 
loading and unloading using the same rate, requiring the 
operator to do no more than press a button or throw a 
switch. Lastly, pursue the use of the external LVDT, this 
will eliminate the need for correcting the deflection 


reading and speed up the rate of data acquisition. 
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Figure 1. System Layout 
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J-Integral 





Figure 2. Pictorial Representation of J-Integral 
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Figure 4. Specimen Drawing 
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Figure 8. Tracing of HY80-5B Fracture Surface 
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APPENDIX A 


TESTING METHODS 


A. INTRODUCTION 
The steps required to perform a J-integral test using a 
single asapecimen approach can be divided into three smajor 
categories, pre-test, test, and post test. Each of these 
categories may be further divided as shown below. 
1. Pre-test 
a. Machine specimen 
b. Precrack sapecimen by fatigue loading 
2. Test 
a. Run J-integral test 
3. Post-test 
a. Heat tint specimen and break apart 
b. Measure crack extension and verify results 
Fach process is described in greater detail in thia 
appendix. A working knowledge of the equipment is 
presupposed and awareness of all safety precautions is 


required before operating the NTS. 


B. PRE-TEST 
a. Machine Specimen 
ASTM Standard E813-81 provides very specific guidance on 


specimen dimensions. All dimensions are related to the 





width «W). Fatigue precracking ia enhanced if the radius at 
the tip of the notch is kept to approximately 0.003 inches. 
This waa accomplished by machining with a apecial cutter 
which was used to make only the final passea in the notch. 
Figure 2 can be consulted for the actual dimensiona used to 


machine the specimens used in this work. 


b. Precracking 

The following step by atep guide is provided to aasist 
operators. 
1. Before turning on the MTS machine ensure the proper 
range card is insatalled in the strain signal conditioner for 
the clip gage to be used. 
2. Turn on power to all equipment with the computer being 
switched on laat. 
3. Load and run "“PRECRACK® program on the HP-9826 coaputer. 


4. Select range scales on MTS; 


example: 10 v 4000 lbs 20% of full load Crange #3) 


10 v 0.02 in 10% of full clip gage (range #4) 


10 v = 2.5 in 50% of full stroke (range #2) 
S. Preasa kO-key on computer and enter specimen data as 
requested. 
6. Preasa k2-key and enter conversion data asa requeated. 
7. Preas kS-key to obtain a table of COD readings for 
various loada and crack lengths. 


8. Set READOUT knoba on the MTS Controller to be: 
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top = load 


middle = satrain 


bottom = stroke 
Set all to (+) polarity. See Figure Al. This ensures 
channels of the MTS machine are properly aligned with the 
data acquisition systen. 
9. Set MTS limit detector switches based on the information 
just generated by the progran. Limit switch settings are 
determined by the simple formula: 
Limit switch = L #* 10 / FSR 
where 
L = desired limit in English units 
FSR = ful) scale range of current operating 
being used in English units 
example “PRECRACK” determined: 
Pmax = 1668 1bs set load limit to 1600 lbs 
for specimen HY80O-S5B 
load limits: 
upper limit --- small positive load for example 40 lbs 
(should never be positive) 
+40 lba = 10 v / 4000 lbs = +0.1 v 
therefore set upper limit knob to +0.1 v 
lower limit --- 1600 lbs maximuna load 
-1600 lbs * 10 vy / 4000 lbs = -4.0 v 


therefore set lower limit knob to -4.0 v 
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strain limits: 
upper limit --- 0.01 in beyond maximum expected during 
precracking 
#+0.01 in * 10 v / 0.02 in = +5.0 v 
lower limit -~- some small value for example -0.1 v 
Cshould never be negative) 
stroke limits: 
upper limit --- 1.5 in to allow clearance for installing 
specimen and yet stop the MTS if this 
limit is reached 
*1.5 in * 10 v / 2.5 in = *6.0 v 
lower limit --- 0.25 in of compressive actuator raovenent 
is more than sufficient for precracking 
“0.25 in #* 10 v / 2.5 in = -1.0 v 
10. Adjust the upper roller supports as required to obtain a4 
span-to-width ratio (S/W) of 4.0. 
11. Install specimen (with the clip gage already installed) 
using alignment collars to properly center. 
12. Zero the load and strain voltages. 
13. Reset and clear all limit switches (all green lights on 
controller panel should be showing). 
14. Using load control energize hydraulics and apply a low 
€<100 1bs)d compressive load to the specimen by SLOWLY 


turning the SETPOINT knob counter-clockwise. 
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15. Slowly decrase the compressive load on the specimen to 

between 5 and 10 1ba by turning the SETPOINT knob clockwise 

and then zero the stroke voltage reading. 

YOU ARE NOW READY TO OBTAIN AN UNLOADING COD COMPLIANCE 

MEASUREMENT OF THE CRACK LENGTH 

16. Press k4-key on the computer and proceed with taking 

data. 

17. Slowly, and as smoothly as possible increase the load to 

approximately SOO Ilbs and return to approximately 20 lbs 

ending the loading cycle with a slight (10 lbs) increase to 

trip the computer out of the data-taking loop. 

The computer will diaplay all pertinent crack length 
information and ask if you want another measurement. You 
may repeat this process as often as you like, just be aure 
the maximun load is lower than that required to extend the 
crack. 

18. Make necessary adjustaents to the MTS such that the 
strain limit switch will trip the MTS when the crack extendas 
to SOx of the deaired extension. This isa determined by 


interpolating results from the table obtained earlier. (See 


Step 7). 
NOTE: As the crack propagates the acceptable maximum 
load decreases. For the specimen, HY80-SB, used to test 


this program a Prax = 1000 lbs was suitable to initiate the 
crack and propagate it to SOx of the desired extension 


(a = 0.506 in). 
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Example: 


ier 


50% crack extension # crack length (a) of 0.506 in 


from table: 


COD = 1.7 volts at 1000 lba and a = 0.506 in 
Therefore set UPPER STRAIN LIMIT SWITCH to +1.7 v. 
YOU ARE NOW READY TO BEGIN CYCLIC LOADING 


Setup MTS: 


DISPLAY: 


a. Setup to read maxinum and rmaininum load in 1l1ba. 


FUNCTION GENERATOR: 


a. CONTROL MODE: Renotea 
b. RATE 1: 10 Ha 

c. RATE 2: N/A 

qd. HOLD BUTTON: In 

e. BREAKPOINT: N/A 

£. HAVER SINE: In 

ge. INVERT: In 


CONTROLLER : 


a. SPAN 1: 90.0 
b. SPAN 2: 0.0 


COUNTER: 


20. 


2i. 


the 


22. 


a. COUNTER INPUT: Oasacillator 

b. COUNTER MULTIPLIER: X100 

c. Press RESET 

Slowly increase the SETPOINT to Pain. 

Preas the RUN button on the CONTROLLER and then release 
HOLD button of the FUNCTION GENERATOR. 


Slowly increase the span of the cyclic loading using 


SPAN 1 knob to the desired Pmax. 
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The MTS will continue cyclic loading until it is 
manually atopped or a limit detector ia tripped. Most 
likely the COD limit detector will trip at which time the 
machine will will atop and unload the specimen. 

23. A new crack length measurement can now be obtained 
ugaing the same manual process as described in Step 16 and 
Step i7. 

24. Readjust limit awitches and resume fatigue loading as 
done in Steps 18 through 22. 

25. Continue this process until the crack reaches’ the 


desired length. 


C. J-INTEGRAL TEST 

The following ateps are required to run a J-integral 
test. They are very similar to those used in precracking. 
The major change is that a ramp function under  asatroke 
control is now used to control the MTS machine. 
1. Before turning on the MTS machine ensure the proper 
range card is inatalled in the atrain signal conditioner for 
the clip gage to be used. 
2. Turn on power to all equipment with the computer being 
switched on last. 
3. Load and run “J_INTGL” prograrm on the HP-9826 computer. 
4. Select range scales on MTS; 


exarnple: 10 v = 4000 lbs 20% of full load (range #3) 
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10 v = 0.02 in 10% of full clipgage (range #4) 

10 v = 0.5 in 10% of full atroke (range #2) 
3S. Preas kO-key on computer and enter specimen data aa 
requested. 
6. Preas ki-key and enter conversion data as requested. 
7. Preaa k3-key and follow instructions for creating a data 
file on an initialized data diskette. See HP-Mannual on how 
to initialize a diakette if necessary. 
8. Preas k3-key to setup the plotter. 
9. Press kS-key and obtain the maximum load you may apply 
while measuring a0. 

YOU ARE NOW READY TO SETUP THE MTS MACHINE 


10. Set READOUT knobs on the MTS Controller to be: 


top = load 
middle = strain 
bottom = stroke 


Set all to (+) polarity. See Figure Al. 
11. Set MTS limit detectors switches. 
example: 
Pmax = 895 lbs for HY80-5B (see Appendix H) 
load limita: 
upper limit --- small positive load (40 lbs) 
(should never be positive) 
+40 lba # 10 v / 4000 lba = +0.1 v 
lower limit --- 895 lbs conpreasion maxiaum 


-895 lba * 10 v / 4000 lbs = -2.125 v 
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strain limita: 


upper limit --- want to be able to read full scale 
on range #2 (0.1 ind 


*O.lin = 10 v / 0.1 in = +10 v 
lower limit --- a very small value set to -0.1 v 
€ should never be negative) 
stroke limits: 
upper limit --- full scale (*10.0 v) on this naost 


sensitive scale 


lower limit --- a maximum compression of 0.25 in is 
sufficient for the test and will atop 
the actuator before it can cause any 
serious damage to the specimen holder 
-0.25 in * 10 v / 0.5 in » -5.0 v 

12. Adjust the upper roller supports as required to obtain a 

span-to-width ratio of 4.0. 

13. Install the specimen (with the clip gage already 

installed) using the alignment collara to properly center. 

14. Zero the load and strain voltages. 

15. Reset and clear all limit awitched ‘all green lights on 

controller panel should be showing). 

16. Using load control energize hydraulics and apply a low 

€<100 lbs) compressive load to the specimen by SLOWLY 

turning the SETPOINT knob counter-clockwise. 

17. Decrease the load to between 5 and 10 1ba and zero the 


atroke voltage. 
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18. Adjust the MTS to the following settings: 
DISPLAY: 
a. Setup to read load in lbs and atroke in inches. 
FUNCTION GENERATOR > 
a. CONTROL MODE: Remote 
be. RATE 1: 6000 (sec) 
as determined by: desire 0.005 in/min loading rate 
therefore; 
0.5 in (full scale) / 0.005 in/mim * 60 sec/min = 6000 sec 
c. RATE 2: N/A 
d. HOLD BUTTON: In 
e. BREAKPOINT: N/A 
£. RAMP BUTTON: In 
ge. INVERT: IN 
CONTROLLER: 


a. SPAN i: 10.0 (full scale) 
b. SPAN 2: 0.0 


COUNTER: N/A 
YOU ARE NOW READY TO MEASURE THE INITIAL CRACK LENGTH 

19. The computer will begin taking data after the ‘CONTINUE’ 
key is preased. The MTS will begin loading the specimen 
after the MTS’ RUN switch is pressed and the HOLD button isa 
released. 

20. When the desired maximum load (approximately SOO lbs) is 
reached; depress the HOLD button, awitch to load control, 
and manually decrease the load (to approximately 20 lbs) by 
SLOWLY turning the SETPOINT knob clockwise. As with 
precracking you must end the unloading with a small increase 


in load (10 lbs from where you stop) to signal the computer 


you have corapleted the load/unload cycle. The computer will 
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plot the data as you go, ensure you remain on the linear 
portion of the loading curve. 

21. The coaputar will diasaplay all the pertinent crack length 
information. You may repeat this operation by switching 
back to atroke control and repeating Steps 19 and 20 as 
often as you like as long as you remain on the linear 
portion of the loading curve. 

22. After the computer has stored the desired initial crack 
length (aO) measurement on data diskette, you are ready to 
begin propagating the crack. 

23. Enaure you are in atroke control with near zero (<10 
lbs) load and the RUN light ia on with the HOLD button 
depressed. 

24. Press the k7-key on the computer follow instructions and 
proceed to collecting data (data acquisition system will be 
continuosly aswitching channels and the plotter will be 
plotting), now release the HOLD button on the MTS. 

25. When you reach a point where you wish to unload the 
specimen; press the HOLD button, switch to load control, and 
Slowly decrease the load. The computer will beep when you 
have unloaded the specimen by 10x of the expected maximum 
load. Switch back to stroke control and release the HOLD 
button to reload the specimen. 

26. When the load returns to the point from which the 


unloading began the coauputer will display and printout all 
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the information it haa determined with regarda to the 
current unloading cycle. 

27. Thia process of switching control modes (Steps 25 and 
Step 26) are repeated for every unloading cycle until the 
you decide to atop by pressing the k9-key on the computer or 
one of the limit detectors shuta off the MTS machine. 
Remember, you must have atleast four data points between the 
offsets one of which liea near the blunting line on the J-R 
curve. 

28. The k9-key interrupts the conputer program and tella it 
to atop taking data and write the data it has collected on 
the data diskette. 

29. The computer will then display a table of the collected 
J versua delta_a and ask if you deaire a plot. Simply follow 
the directions from the coaputer and you will obtain the J-R 


curve and the value of JIc. 


D. POST-TEST 
The poat-teat includes verifying the final crack length, 
and validating the other parameters as discussed in ASTM 
Standard E£813-81 auch as: 
by and Bain > 25 * J[Ic / Fa 

where 

by = uncracked ligament after precracking 

Bain = minimum apecimen thickness 


Fa = flow stress 
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and that the slope of the linear regression fit to the data 
is less than the flow stress (dJ/da < Fa). This data is 
obtained directly from the "“J_INTGL” program output. 

Heat tinting the specimen, 10 minutes at 300C for most 
steels, aids in the visual identification of the amount of 
crack extension. After the specimen is removed from the 
oven and air cooled, place it in liquid nitrogen to assist 
in obtaining a brittle fracture. After breaking the 


apecimen apart, place the pieces in alcohol while returning 


to roorn terperature. This will prevent water vapor from 
condensing and freezing on the surface and mininize 
oxidation. 


A photographic enlargement of the surface can now be 
obtained and the crack length can easily be verified. The 
final crack extension must be within 15% of that recorded 


during the test. 
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APPENDIX B 


TEST RECORD SHEET 


THREE POINT BEND SPECIMEN 


TEST REPORT 


SPECIMEN 1.5. MATERTAL ORTENTATION 


Date Operator Test Teap 
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End _ test 
<9} bi 
y__ 4 ae 
orrec 
E [ese [ORE 
Setup_plotter 


NO jig 


External plotter 


Main_prog 
Line # 9380 












Line # 9400 


Line # 9480 





Prt_results 


Lod2 .98*Lod9? ; 















Line #10430 
(calc area) 


Line #10440 
(reassign 
variables) 












No crack ext 
I1=0 
% 


a od2 Lod9-LodS* 
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APPENDIX E 


PRECRACK LISTING 


PRECRACKING PROGRAM FOR THREE SOINT 
BEND SPECIMEN. THIS PROGRAM INCORPORATES 
THE USE OF THE FOLLOWING EQUIPMENT: 
tee =9826 COMPUTER 
2. HP-3497A DATA ACQUISITION SYSTEM 
3. HP-2671G THERMAL PRINTER 
4. MIS MACHINE 


THE PROGRAM USES BOTH ENGLISH AND METRIC 
UNITS. ALL METRIC VARIABLES ARE DENOTED 
WITH AN m IN THEIR NAME ie Hm=WIDTH IN 
METRIC UNITS (mm IN THIS CASED. 


’ 
’ 

' 

' 

' 

' 

' 

' 

’ 

’ 

’ 

' 

' 

! THIS PROGRAM WAS WRITTEN BY: 
! LCDR W. K. TRITCHLER, USN 
VERSION 1.0 DATE 11/33 

' 

' 

C 

' 

? 

’ 

' 

9 

: 


CM /J_intgl/ FI 


Dim Lodm(2000) ,Deflm¢ 2000) .Codm( 2000) these arrays will be dimensioned 
required in the New_a routine. 

Lodm(I) = LOAD IN NEWTONS 

Deflm(I) =~ DEFLECTION IN mn 

Codm¢(I) = CRACK GPENING DISPLACEMENT IN mm 


oIM Coeff (6) .G¢(6) .AND¢ 100) 
DUMP DEVICE [S 706 


Dvm=709 

CLEAR Dum 

PRINT CHRS(¢(12) ' CLEARS CRT 

PRINT "°° 

PRINT “ PRECRACKING PROGRAM” 

PRINT “* ° VERSION 1.0" 

PRINT “°° 

AY “ts WRITTEN BY LCDR W.K. TRITCHLER™ 
N 9@ 6@ 

PRINT PLEASE SELECT A k-KEY”™ 

' BDO DBDDD SH DH OS OP BOD @SOSSRHP DROP DGSTDEVDPIOOCBHPLeE SBO 

Begin: !?! 


ON KEY 0 LABEL “SPEC_DATA” GOSUB Spec_data 
ON KEY 2 LABEL “CONV_DATA” GOSUB Conv_data 
ON KEY 4 LABEL ~NEW _A” GOSUB son a 

ay KEY S LABEL “a_TABLE” GOSUB A_table 


’ DD DPRBHBDPDPHRPPBDSHMBEPSOMBSPH VPSwBDSHoPSOOOVPR ROHR DHSS OLD 
’ 

‘ WAITS FOR KEY INPUT 
Gato Soin 


2D D-DD WD 4D 1D 4D 1-48 4D 1D 1D 1D 10 10 1D ED-10 2D 0D ED 1D 1D 20 20 1D 20 1D 1D 1D 26 1D 1D 1 1D 1D 1D 1 2D 


: 


Spec_data: ! INPUTS REQUIRED eo DATA 


PRINT CHRS(12) ! CLEARS 

PRINT “INPUT SPECIMEN IDENTIFICATION (ID):" 
INPUT DS 

PRINT “* 
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630 PRINT "J-INTEGRAL TEST for °;D0S 
510 PRINT “’ 
a BRENT "INPUT SPECIMEN DATA USING ENGLISH UNITS (LBS, PNGHES sete.) 
630 PRINT ‘*"* 
of 0 PRINT “WHAT IS THE SPECIMEN’S WIDTH. NOTCH LENGTH, Bmax, Bmin (THICKNESS) 
65) INPUT W,A.Bmax.3min 
660 PRINT "°°" 
670 PRINT “WHAT IS THE SPECIMEN’S HALF SPAN AND POISSON’S RATIO?” 
680 INPUT Hs,Nu 
630 PRINT °* 
700 PRINT “WHAT IS THE SPECIMEN’S ELASTIC MODULUS. YIELD STRENGTH AND ULT 
IMATE STRENGTH?” 
7140 INPUT E.Ys.Uts 
720 PRINT CHRS(12) ' CLEARS CRIT 
730 Fs=(Y¥stUts)/2 
740 PRINT “USING ENGLISH UNITS Cin. lbs.psi)” 
750 PRINT “THE FOLLOWING HAS BEEN ENTERED for: “;05 
760 PRINT “WIDTH =";W;” NOTCH LENGTH =":A 
770 PRINT “Bmax *"3Bmax;" Bmin =";3Bmin 
780 PRINT "HALF SPAN =":Hs;" PONISSON’S RATIO ="*;Nu 
730 PRINT “YIELD STRENGTH ="3Y¥s:° ULTIMATE TS =";Uts 
800 PRINT “ELASTIC MODULUS =";E;" FLOW STRESS ="sFs 
810 PRINT “°° 
82 PRINT “THIS PROGRAM ALLOWS NO EASY RECOVERY, IF ANY DATA IS INCORRECT YCU 
MUST RE-ENTER "; 
830 PRINT “EVERYTHING. ee 
840 PRINT ° 
850 PRINT “IS ALL DATA CORRECT?” 
860 INPUT HS 
870 IF HS="NO" OR HS="N* THEN 550 ' REPEATS Spec_data 
880 PRINT °* 
830 PRINT “DO YOU WANT A HARD COPY?" 
’ 900 INPUT HS 
ao of HS="N" QR HS="N" THEN 1060 
aap ence. IS 706 * TOGGLES PRINTER TO THERMAL PRINTER 
94 PRINT 
350 PRINT “USING ENGLISH UNITS ¢in.Ilbs.psi) THE FOLLORING INFORMATION” 
360 PRINT “HAS BEEN ENTERED FOR: °*;DS 
970 PRINT 
980 PRINT "WIDTH ea oe NOTCH LENGTH =";A 
990 PRINT “Bmax =° *sBmaxs" Bmin ="s3Bmin 
10060 PRINT “HALF SPAN ="*sHs:"* POISSON’S RATIO ="3Nua 
1010 PRINT “YIELD STRENGTH ="°¥s;"* ULTIMATE TS ="3Uts 
1020 PRINT “ELASTIC MODULUS =";E;" FLOW STRESS =":Fs 
1030 PRINT ** 
on PRINTER Is 1 ' TOGGLES PRINTER BACK TO CRT 
ee cia CHRS( 12) ' CLEARS CRT 
1080 ¢' CONVERSION TO METRIC UNITS 
1090 Conlen=25.4 ' jin = 25.40mm 
1100 Km=W#Conlen 
1110 Am=A#Conlen 
1120 Bmaxm=Bmax*Conlen 
1130 eae men conten 
1140 Hsm=Hs*Conle 
1150 Conpress=b. 39475293 ' tpsi = 6.89475293kPa (HP PROGRAMING BOOK) 
1160 Em=E*Conpress 


Sl 





1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 


Ysm*Y¥s*Conpress 
Utsm=Uts*Conpress 


D2=1 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 


IF HS="N*° OR HS="NO" THEN 2450 


PRINT 
PRINT 


Fsm=Fs*Conpress 

9 

PRINT “CONVERTING THIS DATA TO METRIC UNITS (mm.N,KPad for. 0S*7 YIEL 
PRINT “HIDTH =#":Wm;”° NOTCH LENGTH 7'°:Am 

PRINT “Bmax =":Bmaxm:” Bain =" ;Bminm 

PRINT “HALF SPAN ='°:Hsm;” POISSON’S RATIO =";Nu 

PRINT “YIELD STRENGTH =";Ysm:" UTILMATE TS =";Utsm 

art "ELASTIC MODULUS =";Em;" FLOW STRESS =";Fsm 

N tees 

PRINT “OO YOU KISH TO HAVE A PRINT QUT OF THIS DATA?" 

INPUT GS 

IF GS="N" CR G3S="NO" THEN 1410 

2 

PRINTER [S 706 TOGGLES TO THERMAL PRINTER 

PRINT “CONVERTING THIS DATA TO METRIC UNITS (mm.N,KPa) FOR: “:0S 
PRINT “WIDTH =":Wm;" NOTCH LENGTH =";Am 

PRINT “Bmax =";Bmaxm:” Bmin =";Bminm 

PRINT "YIELD STRENGTH *";:Ysms" ULTIMATE TS =";:Utsm 

BRENT "ELASTIC MODULUS =";Em;" FLOW STRESS ©";Fsm 

N e009 

PRINTER IS 1 ' TOGGLES PRINTER TO CRT 

PRINT CHR3(12) TGLECREERT 

' 

Flag3=0 

Al=A 

GOSUB Params 

! 

RETURN 

’ 

; 9D 2D 2D 0 0 DO - 2D 1 2D 4D 20 1 2 2D 2D 20 2 > DD 4 4D 4 DO OS 0 2 WD 0 2 De Oe we 

! 
New_a: * THIS ROUTINE COMPUTES NEW CRACK LENGTH AND NEW PARAMETERS 
! :¢. 

PeiPacPeaxaa(inaxmetwiem) 2/Bmaxm  ' B EFFECTIVE FROM ASTM E£813-81 


' INITIALIZES COUNTER 
CHRS (12) ' CLEARS CRT 
“THIS ROUTINE WILL OBTAIN @URRENT CRACK LENGTH” 
* AND NEW PARAMETERS FOR PRECRACKING.” 
“BEFORE YOU CAN RUN THIS ROUTINE YOU MUST:” 
1. ENTER SPECIMEN DATA” 
“ 2. ENTER CONVERSION DATA” 
3. SETUP MTS FOR MANUAL LOADING OF THE”; 


cw SPECIMEN.” 
" ARE YOU READY?” 
HS 


! RETURNS TO k-KEY SELECTION 
CHRS( 12) CLEARS CRT 
"YOU MAY LOAD AND UNLOAD SPECIMEN AS OFTEN AS YQU LIKE -- BUT YOU 


UST REMAIN "; 


1690 
1700 
1710 
1720 
1730 
1740 


PRINT 


“BELOW CURRENT PMAX OF’: Pmax; “Lbs” 


DIM Lod(2000) .Def1¢2000) ,Cod¢2000) 
DIM Lodm(2000) ,DefIlm¢(2000) ,Codm(2000) 


PRINT 


"WHEN YOU ARE READY TO BEGIN PRESS THE ‘ENTER’ OR ‘CONTINUE’ KEY.” 


INPUT HS 


PRINT CHRS(12) 


PSCEEARS seRT 
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1750 
1760 
1770 
1780 
1730 
1800 
1810 
1820 
1830 
1840 
1850 
1360 
13790 


Kai ' USED AS COUNTER 
[=| ' USED AS INDEX 
Filagl=0 
Flag2=0 
WHILE Flagt=0 ' ESTABLISHES ONE LOAD-UNLOAD DATA TAKING LOOP 
GOSUB Data_acq 
END WHILE 
Flag2=0 ' RESETS FLAG2 
WHILE Flaq2=0 
GOSUB Data_aca 
END WHILE 
' THIS SECTION CAN BE USED TO VIEW THE COLLECTED DATA 3Y REMOVING THE REMA 


RK CHARACTER ¢!). 


1880 
1890 
1900 
1910 
f= 20 
1930 
1940 
1950 
1960 
1970 
1980 
fo30 
2000 
2010 
2020 
2030 
2040 
29050 
2060 
2070 
2080 
2030 
2100 
2110 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 
2200 
2210 
2220 
2200 
2240 


2280 
2290 
2300 
2310 
2320 


' FOR [=1 TO K 
’ PRINT USING 1960;I,Lodm¢I) ,Defim¢(I) ,Codm(T) 
' IMAGE X.20,4X,60.0,4X%.30.7D.7X.5D./7D 
f NEXT I 
' WAIT S 
Ki=1 
K3=K 
GOSUB Slope SLOPE = {/COD COMPLIANCE 
GOSUB Crack_length ! USES COD COMPLIANCE 
p 
A0(D2)=CrackIm/25.4 ’ CONVERTS TO INCHES 
AlL=A0«(D2>) 
? 
' FIND 95% CONFIDENCE INTERVAL FOR CRACK LENGTH 
Stiffm=StiffimeTQa4 
GOSUB Crack_length 
Amax=Crack1lm/25.4 * CONVERTS TQ INCHES 
Stiffm=Stif fm-2*198U4 
GOSUB Crack_length 
Amin=Crack im/25.4 ' CONVERTS TG INCHES 
PRINT CHRS(12) * CLEARS CRT A. 
PRINT “THE FOLLOWING INFORMATION HAS BEEN OBTAINED:” 
PRINT ‘°° 
PRINT “CRACK LENGTH =":A0¢(D2):"in FROM MEASUREMENT #°:D2 
Nt "WITH 95% CONFIDENCE INTERVAL FROM’;Amin;” TO’; Amax 
PRINT “USING COD COMPLIANCE OF:”";Compl:"“in/Lb =""-Complm; “mm/N” 
aN "WITH A CORRELATION OF :*;Correl;” FROM’;:K;°DATA POINTS” 
IN oe ee 
PRINT “DQ YOU WANT A HARD COPY?” 
INPUT GS 
IF GS="N" OR GS="NO" THEN 2340 
PRINTER IS 706 * TOGGLES TO THERMAL PRINTER 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


PRINTER i> 1 


“THE FOLLOWING INFORMATION HAS BEEN OBTAINED: 

"CRACK LENGTH =";A0<(D2);:"in FROM MEASUREMENT #°°;D2 

“WITH 9SZCONFIDENCE INTERVAL FROM’ ;:Amin:” TO" :sAmax: in” 
“USING COD COMPLIANCE OF :“":Compl;:"in/Lb =";:Comolm: ’mm/N” 
"WITH A CORRELATION OF:":Correl;”"°FROM’:K;°* DATA POINTS.” 


! TOGGLES PRINTER TO CRT 


a8, 





2330 
2340 
Za50 
2360 
2370 
2380 
2390 
2400 
2410 


BRUNT "* 

PRINT “* DO YOU WANT ANOTHER MEASUREMENT?” 
INPUT FS 

D2=D2+! ' INCREMENTS COUNTER 


IF FS="Y" OR FS="“YES" THEN 1670 

! 

IF AQ(D2-1)>Qpta THEN GOSUB End_precrack ! PRECRACKING IS COMPLETE 

’ 

IF AQCD2-1)>.4*DelatA THEN F lag3=! ' SETS FLAG INDICATING IN LAST 60% O 


F PRECRACKING 


2420 
2430 
2440 
2450 
2460 
2470 
2430 
2490 
2500 
2510 
2520 
2530 
2540 
2550 
2560 
2a 0 


’ 
GOSUB Params 
, 


RETURN 


’ 32 2 4D 10 WD WD DD WD Oe 1 OD DDD DOOD DHE DHD DBT HH HD HD 
? 


Data_aca: ! DATA ACQUISITION ROUTINE 


OUTPUT Dym:"“VRS AFI AL3" 
' 


QUTPUT Dvm;"AI! VTI" ' READS LOAD 
ENTER Dvm;Lod¢I) 

Lodm(I)=Lod(I)#Slm¢t) 

? 


QUTPUT Dvm:"AI2 YTI" * READS COD 
Biter Dum:Godt 1? 

Codm(I)=Cod¢ I) #Slm(2) 

? 


OUTPUT Dvm:"AI3 VTI" t READS DEFLECTION 

BHIPER Dum: Defl¢1) 

Def lm¢I) =Def l¢ TI) #Slm¢3) 

? 

' TEST DATA TO BREAK OUT OF DATA ACG LOOP 

; 

IF I=! THEN 2720 - !¢ SKIPS IF FIRST DATA POINT 

Lodif fm=Lodm(I)-Lodm¢(I-1) 
ialadirim<-25.0 THEN Flagt=t t MUST BE UNLOADING 
IF Lodiffm>25.0 THEN Flag2=1 ' MUST BE RELOADING. 

Lodiffm=0 

i) 

[=I+t t INCREMENTS INDEX 

K=K+1 * INCREMENT COUNTER 

ij 

RETURN 

‘ 


’ 90 WD TE WD WD 10 TE TE OD HE WD ED ED WD 2D 1 WD 1 DD OD OD EE ED DD WD 2 wD OD 1 1 


Crack_length: ' CALCULATES CRACK LENGTH USING HUDAK ‘DITNSRDC PROG) 


' NOTE: THIS ROUTINE USES COD COMPLIANCE METHOD VICE DEFLECTION COMPLIA 


FOR I=! TQ 6 
READ Coef¢I) 
NEXT I 


Comp lm=1/Stiffm 
C=2Wm=Bef fmeEm=Complm/ (Hsm#1000) 
C=1/¢SQR(C) +1) 
G(t)=C 
FOR [=2 70 5 
GCI) =GCI-1)#GC¢1) 
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2910 NEXT [ 

22920 F=Coefit) 

2930 ##$>FFOR I! TQ 5 

2940 FeF +Coef(I+1)#GCI) 

2950 NEXT I 

2960 Crack Im=F «Wm 

2370 RESTORE 

2980 DATA J.998265,-3.81662,-1.80596,32.31038,-44.15665.-52.67376 
' 


3000 RETURN 
t 


3010 ' 

3020 I Re PHHHSHHHHH HHH HH HHH KSHHHORHO SOOO 

3030 ! 

3040 Conv_data: ' INPUTS CONVERSION DATA FOR MEASURING DEVICES 
3050 PRINT CHRS(12) SEREARS CRI 


3060 PRINT “INPUT: LOAD RANGE, STRAIN RANGE, STROKE RANGE” 
3070 PRINT “" 

3080 PRINT * (BY PERCENT ie. 10.20.50 OR 100)” 

30390 INPUT Lr,Cr,Sr 

3100 PRINT “°° 

3110 PRINT “LOAD RANGE = “:Lr 

3120 PRINT “STRAIN (COD) RANGE = “:Cr 

3130 PRINT “STROKE (DEFLECTION) RANGE = “:Sr 


3140 Lmax=20000 ' MAX LOAD RANGE IN Clbs). 

3150 Cmax=.200 ! MAX COD RANGE IN Cin). 

3160 Smax25.0 ' MAX STROKE RANGE IN (in). 

ae V=-19 ' MAX VOLTAGe IN (volts). 
180 ' 


3:90 ' MAXIMUM OPERATING POINTS ARE: 
3200 Olmax=Lmax*Lr/100 
3210 Ocmax=Cmax«Cr/100"#(-1.0) ' -1 CHANGES SIGN 
3220 Osmax=Smax*Sr/100 
9 


3240 ' CONVERSION CONSTANTS ARE: CENGLISH UNITS) 
3250 $1¢1)=Olmax/V . 
3260 $1¢2)=Ocmax/V 
3270 $1¢3)=UOsmax/V 
3280 ' CONVERTING TO METRIC: 
3290 Conlen=25.4 ' 25.4mm = Jin 
3300 Conib=4.448221615 ' 4,44822161SN = 1I1b (HP PROGRAMING BOOK) 
3310 SIm€1)#S1¢1)*Vonilb 
3320 Slm(2)=S1¢2)#Conlen 
3330 Slm¢(3)2S1(3)#Conlen 
' 


3350 PRINT °” 
3360 PRINT °™ 3 
3370 PRINT “YOUR ENTRIES HAVE PRODUCED THE FOLLOWING CONVERSION DATA: 


See0 PRINT “* 

3390 PRINT “FOR LOAD: “;Slm¢1);"NEWTON/VOLT C"';S1C1) 3s "Lbs/V)” 
3400 PRINT “* : 

3410 PRINT “FOR COD: “°;Slm(2);"mm/V Ch sS1CZ) 3" in/V)" 

3420 PRINT “°° 

3430 PRINT “FOR DEFLECTION: “;Silm(3)3"%mm/V (°3S163) 3" in/V)” 
3440 PRINT °* 

34350 PRINT “WOULD YOU LIKE A HARD COPY?” 

3460 INPUT HS 

34790 [IF HS="N° QR HS="NO" THEN 3570 

3480 ! 

3490 PRINTER IS 706 ! TOGGLES PRINTER 
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3500 


3990 
4000 


4040 
4050 
4060 
4070 
4080 
4090 


Tay SN) 

PRINT “ CONVERSION DATA” 

PRINT “FOR LOAD: "*;Slm¢1):"N/V Coss Gis bs/V)" 

PRINT “FOR COD:":Sim¢(2) :"mmsV (eral C2) 3 in/V?” 

PRINT “FOR DEFLECTION: ";Sim(3):"mm/V ees (3) 5 7 insV) 
PRINTER [S 1 ! TOGGLES PRINTER 

' 


PRINT CHRS(12) 
 } 

RETURN 

? 


s 
’ DRM ARAAHHHASHBHHHDHSHSAHHMPSAHSHHKHAHVRHHHHHRHHBHAHHH Hwee 
’ 


Slope: ! LINEAR LEAST SGUARES FIT 


PRINT CHRS(12) ' CLEAR CRT 

! INITIALIZE VARIABLES 

Sumx=0 ! SUM OF X VALUES 

Sumx2=0 ’ SUM OF X SGUARED VALUES 

Sumy =9 ' SUM OF Y VALUES 

Suny2=0 ! SUM OF Y SQUARED VALUES 

Sumxy #0) ! SUM OF THE PRODUCT 

FOR [=K1 TO K3 ' KI=INDEX OF FIRST, K3=INDEX OF LAST DATA POINT USED. 
' 
X=Codm(T) 
Y=Lodm( 1) 
' 


Sumx =Sumx +X 

Sumy =Sumy+Y 

Sumx2=Sumx2+X#X 

Sumy2=SumyZ2+Y¥*¥ 

Sumxy =Sumxy+X*/Y 
NEXT I 
K3=K3-P2 f USED TG CALCULATE NUMBER OF DATA POINTS WITHIN SOUNDS. 
K4=K3-K1+1 ' K4=TOTAL « DATA POINTS USED IN CALC 
Qat =Sumx2-Sumx*Sumx/K4 
Qa2=Sumxy-Sumx*Sumy/K4 
ee ey gos Sumy 7S .2> 


M=Qa2/Ga! ! SLOPE 
} 


IF M<O THEN 4190 ' IF SLOPE < 0, BAD DATA 
' FROM DITNSRDC PROGRAM 
Qt=(K4"Sumx2-Sumx~2)/(K4#(K4-1)) 
Q2=(K4*"Sumy2-Sumy ~2)/(K4"(K4-1)) 

Q3=(K 4-1) #(Q2-M#N#Q1)/(K4-2) 
Q4=SQR(Q3/(Q1"#(K4-1))) 


t 
B=(Sumy—M*Sumx) /K4 ! INTERCEPT 
Correl=M*SQR(Q1/Q2) 

] 


Qa4=SGRCABS( (Sumy2-B*#Sumy-M#Sumxy) /(K4-2)))? 
Sigqmab=Qa4/SQR(ABS (Gal >) 

Sigmam=S igmab#SQR(Sumx2/K4) 

Stiffm=" ‘ STIFF = SLOPE 

Comp lm=!/M ' COD COMPLIANCE = 1/SLOPE 

Comp | =Comp lm#(4.448221615)/25.4 !'! CONVERTS TQ in/LB 
9 


{ STUDENT t TEST 
U4=04 
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4')0 T9=!.96 

4itv IF K4>120 THEN 4179 ' RETURNS 

4120 To=2.] 

4130 IF K4<20 THEN 4170 * RETURNS 

4140 N9=K4-2 

4150 T92N9/(0-.6171094+N9# 510239) 

4160 ' 

4170 RETURN 

4180 ! 

41930 PRINT “POOR DATA-- SLOPE < oO” 

4200 mms GG933~.3 

4210 Correi=0 

4220 ! 

4230 RETURN 

4240 t 

4250 P20 0 ao A ab 10 Sb AD SP 20 90 Sh AD AD AO 1 HD AE 4D AD 20 DA AD AO aD AD AD A OO DD 10 aD OO OO 4D 0 OD 

4260 ' 

4270 Param ® THIS ROUTINE WILL DETERMINE THE CURRENT MAX AND MIN LOAD FAR 
a PRECRACKING 

4280 

4290 ’ NOTE THE USE OF ENGLISH UNITS DUE MTS READOUT BEING IN POUNDS. 
4300 ' INITIALIZE VARIABLES 

4310 Faw=0 

4320 Pmax=0 

4330 Pmin=0 

4340 Delp=0 

4350 Kmaxp=0 

4360 Kminp=f 

4370 Delk=0 

4380 De lke=0 

4390 ' 

4400 Opta=.6«W ' OPTIMUM CRACK LENGTH 

4410 Dela=Opta-Al ¢ EXTENSION REQUIRED 

rie Dela50-=Dela#(.5)+A ' CRACK LENGTH AT SOZ OF EXTENSION 

4430 t 

4440 Pmax=(4/3)#(Bmin#(W-Al)*2*Fs/(2*Hs) ) ' FROM ASTM E813 

4450 IF Flag3=1 THEN Pmax=(.4)#Pmax ’ MAX LOAD DURING LAST 50% OF PRECR 
ACKING (ASTM E813) 

4460 Pmin=Pmax-.95*Pmax ' FROM ASTM E813 

4470 Delp=Pmax-Pmin 

4480 ’ 

Hep0 DETERMINE Delk BASED ON P USING ASTM E399 

ae Fauna 39(AL/W)™ .S#¢1.S9-CAL/W)@(T-AL/H) #(2. 15-3. 93eAL/W +2. 7*(AL/W) 2)) 
4530 Fawd=2"¢(1+2"AL/W)#(1-AL/W) 71.5 

4540 t 

4550 Faw=F awn/F awd 

4560 ' : 

4570 t 

4580 Kmaxp=(Pmax#2*Hs/(BmineW*1.5)) Faw ' FROM ASTM E399 

4530 Kminp=(Pmin#2*#Hs/(Bminew 1.5))*F au ' FROM ASTM E399 

ea Deik=Kmaxp-Kminp 

4620 ' DETERMINE Delke BASED ON MODULUS FROM ASTM E813 

4630 a Flag3=0 THEN Delke=40000 § DURING EARLY EXTENSION 
4640 F Flag3=! THEN Delke=¢(.001)#E ' DURING LAST S0% OF PRECRACKING [CRE 


F. ASTM "£399 & £813] 


4650 


PRINT ‘°° 
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4660 
4670 
4680 
4630 
4700 
4710 
4720 
4730 
4740 
4730 
4760 
4770 
4780 
4790 
4800 
4810 
4820 
4830 
4840 
4850 
4860 
4870 
4880 
4890 
4900 
43910 
4920 
4930 
4940 
4950 
4960 
4970 
4980 
4990 
5000 
3010 
3020 
3030 
5040 
5050 
3060 
3070 
3080 
3090 
5100 
5110 
3120 
Sio0 
3140 
3150 
3160 
5170 
5180 


3190 A 


aR "SascD ON CURRENT INFORMATION: 

INT 9@e¢ 

PRINT “CURRENT CRACK LENGTH =";sAl:"“in 

PRINT “OPTIMUM a0=";Geca:"in' 

PRINT “EXTENSION REQUIRED =";:De!a 

PRINT “CRACK LENGTH AT SOZ% OF EXTENSION a" -DelaS0;° 
PRINT “Pmax ai? :Pmax : “Lbs Pain =' :Pmin: “lbs” 
PRINT “Kmaxp ="3Kmaxp:” Kminp =";Kminp 

PRINT ‘DELTA_K =" Delk:"in’.5” 

PRINT “DELTA_KE =";Delke;:"in*~.5” 


a “FAW =">Fau f SERVES AS A CHECK 

P IN eeee 

PRINT “IF Delk > Delke. Pmax MUST BE DECREASED” 
IF Delk<Delke THEN 4900 ' SKIP IF QKAY 


* COMPUTE NEW Pmax 
Delpi=(Delke/Faw/(2"*Hs) #*BmineW7 1.5) 

Pmax!=Delp!/.95 

Pmint=Pmax1-.95#Pmax1 
Kmax!=¢(Pmax!*2*Hs/(BmineW~1.5))*F aw 

PRINT “THE NEW Pmax BASED ON Delke IS:°;Pmax1;"“Lbs” 
PRINT “NEW Pmin =";Pmin!; “Lbs” 

PRINT “NEW DELTA_P =";Delpi;"Lbs” 

PRINT “NEW Kmax =";Kmax! 

’ 


PRINT * DO YOU WANT A HARD COPY?’ 
INPUT HS 
IF HS="N" OR HS="NO” THEN 5080 
PRINTER IS 706 ' TOGGLES PRINTER 
BRINE “BASED ON CURRENT INFORMATION FOR °;DS;°°;" 
N eeee 
PRINT “CURRENT CRACK LENGTH #":ALl;“in” 
PRINT “OPTIMUN CRACK LENGTH =";Opta:"“in 
PRINT “EXTENSION REQUIRED =";sDela:in” 
PRINT “CRACK LENGTH AT 5S0Z OF EXTENSION =" DelaS0;“in” 
PRINT “Pmax =";Pmax:"Lbs Pmin =" hie PS 
PRINT “Delp =";Delp;"Lbs Kmax 2*":Kmax 
PRINT “DELTA_K =";Delk:"in*.5S wae “KE ="-Delke;s"in~™.5” 
IF Delk<Delke THEN 5090 
PRINT °°” 
PRINT “NEW LOADS BASED ON DELTA_KE ARE:” 
PRINT “Pmax 2";Pmaxi;:"Lbs Pmin =";Pmint;“Lbs” 
PRINT “DELTA_P =";Delp!;"Lbs Kmax =" ;Kmax] 
PRINT °° 
Nee 1 ' TOGGLES PRINTER TQ CRT 
PRINT “MAKE THE NECESSARY ADJUSTMENTS BEFORE CONTINUING.” 
PRINT “WHEN COMPLETED PRESS “CONTINUE’.” 
INPUT cs 


RETURN 

’ 

s 

if DE DD AD ED BD ED ED DD D2 A DD DD 2 2D DD OD 2 DD OD OD 0 0 2 2 2 


! 
table: ! PRODUCES TABULATED VALUES FOR LOAD vs COD VALUES AT VARIQU 


S CRACK LENGTHS 


3200 
3210 
9220 
3230 


! 
PRINT “BEFORE YOU CAN RUN THIS SUBROUTINE YOU MUST ENTER:” 


PRINT °° 1. SPECIMEN DATA” 
PRINT * 2. CONVERSION DATA” 
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e230 
5300 
a2 10 
5320 


PRINT °° 

PRINT “ARE YOU READY?” 

INPUT HS 

IF HS="N" OR HS="NO"° THEN S550 ' RETURNS TO k-KEY SELCETION 
’ 

Bef f =Bmax-(Bmax-Bmin) *2/Bmax ' IAW ASTM E813-81 


PRINT -BEFES =; Bert 
’ 
! USING CRACK LENGTH EQUATION FROM DTNSRDC PROGRAM CALC CRACK LENGTH FOR 


VARIOUS VALUES OF U 


A00=.998265 
Al=-3.81662 
A2=-1.80596 
A3=32.31038 
A4=-44.15665 
AS=-52.67876 
FOR U=.07 TO .15S STEP .005 
Ar =A00 tA1 *U+A2*U"2+A3 "US tA4GHU* 4+A52U"5 
Ck l=Ar«W 
PRINTER IS 706 ‘ TOGGLES TO THERMAL PRINTER 
PRINT “U =": 
PRINT “a/W =":Ars"in CRACK LENGTH =":Ck1l:"in” 
PRINT “LOAD (lbs) COD (in) CGUDeCVaolts?: 
FOR Lod!=200 10 1800 STEP 100 
C=BeffeWet*2/Hs 
Cod! =Lodt#(€(1-U)/U)7~2/C 
Codv=Cod1/S1¢(2) 
PRINT Lod! ,Cod!,Codv 
NEXT Lod! 
NEXT U 
PRINTER rst ' RETURNS TG USING CRT AS PRINTER 


RETURN 


' WD WD Wh Wh WD ID OE I DE ED D-DD I 10 I 8 8 DD 


ij 


End_precract.: f END OF PRECRACKING 
? 


PRINT PRECRACKING IS COMPLETED!” 
PRINT “FINAL CRACK LENGTH =:A0(D2-1) 
PRINT "CORRELATION OF DATA =";Correl 
PRINT * REMOVE SPECIMEN FROM MTS.” 
q 


END 
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APPENDIX E 


J_INTGL LISTING 


J-INTEGRAL FRACTURE TOUGHNESS PROGRAM 
SOR ninco “Ot bEND SPECIMEN. [His 
PROGRAM INCORPORATES THE USE OF THE 
FOLLOWING EQUIPTMENT: 

ne 3925 CONRUTER 

wine soda Weta acOULS LI LON SYSTEM 
wae soso mn SYSTEM VOLIMETER 

2 HP>26716 eeRMAtPR INTER 
Meiiga 22 oo er eu heis 


THIS PROGRAM WAS WRITTEN BY: 
BGDReneet. IRITCHLER. USN 
VERSION |.0 Dele 11/33 


UI BOIN = 


OM /.J_intgl/ Ft 


DIM Lod‘8000) .Def1(8000) .Cod(8000) 

NOTE: Lod, Defl AND Cod ARE ACTUALLY DIMENSIONED IN [nitiai_a ROUTINE 
Lod¢I) = LOAD 

Defl¢I) = DEFLECTION 

Coa(I) = CRACK OPENING DISPLACEMENT (CdD) 


DIM AQ(SO0) .Amin0(50) .Amax0(50) .Como10¢(S0) ,Correl0(50) 
DIM Area9(100).AL¢Ci00) .Coef (6) .G¢5) 

DIM J9¢100) .Delta_a(100) ,Ndefl( 100) 

DIM DataS(1:4)(10) 

DUMP DEVICe IS 706 ' THERMAL PRINTER ADDRESS 


er ee ee 2 a Cle A ee i ee ae ee a | 


Dvm=709 ! DATA ACQUISITION ADDRESS 

CLEAR Dvm 

PRINT CHRS¢(12) ’ CLEARS CRT 

PRINT ‘°° 

PRINT * . J- INTEGRAL TEST PROGRAM” 

PRINT * VERSION 1.0" 

PRINT °° WRITTEN BY LCDR W.K. TRITCHLER” 

PRINT ‘* 7 

PRINT - PLEASE SELECT A k-KEY’ 

‘ ioe ee fe ef ee eee ee ee ee se eee eS ee Ye ee eee eon ole ol ae aes 
Begin: ° 


ON KEY 0 LABEL “SPEC_DATA’ GOSUB Spec_data 
ON KEY 1 LABEL “CONV_DATA’ GOSUB Conv_data 
ON KEY 2 LABEL “DAT_FILE® GOSUB Data_file 

ON KEY 3 LABEL “PLOTTER’ GOSUB Setup_plotter 
ON KEY 4 LABEL “GCLEAR® GOSUB Gelr 

ON KEY S LABEL “a0” GOSUB Initial_a 

ON KEY 7 LABEL “START GOSUB Main_prog 

'(#) ON KEY 8 LABEL “R-CURVE” GOSUB R_curve 
{ (SEE NOTES ON ¢(#) IN R_curve ROUTINE 

on KEY 9,4 CALL Stop_test 


’ SD 2b 2D WD SD 2b 2D OD 2b 2D 20 2D 2D 2D 2 2D Uh OD 2D OD 2D 2 2D WD 2 OD 2S 2 OD OD 2D 2D 2D ED D-DD 2D 2D 2D 2D OD OD 


! 


Spin: ! WAITS FOR KEY INPUT 


COTO Spin 


t bh Ub th BY 2D Ub 2D ID WD 2D 2 2D 2 OD 2D 2D OD OD UD 2D WD 4D OD 1D OD HD OD 1D 2D << 2 OD 2D OD 2 OD OD OO OD OD 
t 
sd 
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600 Gelr: 


' GRAPHICS SCREEN CLEAR 


610 * THIS ROUTINE ALLOWS SINGLE KEY STROKE TO REMOVE A GRAPHICS DISPLAY "ROM 
Cr! 
629 GCLEAR 
Bro RETURN ¢! RETURNS TO «-KEY SELECTION 
40 ' 
650 YEAHS REREEEDE HEE EON ODN ESON DEDEDE DEED EOS 
650 ! 
679 Spec_data: ' INPUTS REGUIRED SPECIMEN DATA 
6&0 PRINT CHRS(12) ' CLEARS CRT 
630 PRINT “INPUT SPECIMEN IDENTIFICATION ¢ID):” 
700 INPUT DS 
7°04 PRINT °° 
720 PRINT “J-INTESRAL TEST for “;DS 
730 PRINT “* 
cas ae * INPUT SPECIMEN DATA USING ENGLISH UNITS (LBS, INCHES, etc.)” 
f Ni #066 
Bee PRINT “WHAT IS THE SPECIMEN’S WIDTH. NOTCH LENGTH, Bmax, Bmin (THICKNESS) 
770 INPUT W,A,Bmax.Bmin 
780 PRINT ‘* 
720 PRINT “WHAT IS THE SPECIMNEN’S HALF SPAN AND POISSON’S RATIO?" 
860 INPUT Hs.Nu 
810 PRINT ‘** 
320 PRINT “WHAT IS THE SPECIMNEN’S ELASTIC MODULUS. YIELD STRENGTH AND ULT 
IMATE STRENGTH?" 
Bac iipst €. %si,Uts 
840 PRINT CHRS(12) ' CLEARS CRI 
850 Fs=(YstUts)/2 
860 PRINT “USING ENGLISH UNITS Cin. lbs.ps1)" 
870 PRINT “THE FOLLOWING HAS BEEN ENTERED for: “;DS 
88) PRINT “WIDTH ="*sW;" NOTCH LENGTH =";A 
890 PRINT “Bmax 2";Bmax:” Bmin =":Bmin 
300 PRINT “HALF SPAN =":Hs;” POISSON’S RATIO =";Nu 
910 PRINT “YIELD STRENGTH #="3Y¥s;:" ULTIMATE TS =";Uts 
na ial “ELASTIC MODULUS =";E;" FLOW STRESS 2";Fs 
940 PRINT "THIS PROGRAM ALLOWS NO EASY RECOVERY, IF ANY DATA IS INCORRECT YOU 
Muot RE-ENTER “: 
950 PRINT “EVERYTHING.” 
960 PRINT °°° 
970 PRINT “IS ALL DATA CORRECT?" 
3980 INPUT HS 
990 IF HS="NO" OR HS="N”" THEN 670 ' REPEATS Spec_data 
1000 PRINT “* 
1010 PRINT “DO YOU WANT A HARD COPY?” 
1020 INPUT HS 
hate IF HS="N° OR HS=#"N" THEN 1180 
ae eee Is 706 ' TOGGLES PRINTER TO THERMAL PRINTER 
N eeeg 
1070 PRINT “USING ENGLISH UNITS (in,.lbs.psi) THE FOLLOWING INFORMATION” 
1080 PRINT “HAS BEEN ENTERED FOR: °:DS 
1090 PRINT ‘*” 
1700 PRINT “WIDTH =":W:" NOTCH LENGTH =";A 
1110 PRINT “Bmax =";:Bmax;" Bmin 2"3Bmin 
1120 PRINT “HALF SPAN =":Hs;"* POISSON’S RATIO ="3Nu 
1130 PRINT “YIELD STRENGTH #":Y¥s:° ULTIMATE TS #=":3Uts 
eG a “ELASTIC MODULUS =";E:" FLOW STRESS #";Fs 
N eeee 


67 





1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1530 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 


PRINTER IS 1 t TOGGLES PRINTER BACK TO CRT 
’ 

PRINT CHR3¢12) ' CLEARS CRT 

9 


* CONVERSION TO METRIC UNITS 

Conlen225.4 ’ Jin = 25.40mm 

WeWeConlen 

A=A#Conlen 

Bmax *Bmax*Conien 

Bmin=Bmain#Conlen 

Hs7Hs*Conlen 

Conpress76.39475293 ' Ypsi = 6.89475293kPa (HP PROGRAMING BOOK) 
E=E#Conpress 

Ys7Ys*Conpress 

Uts2Uts*Conpress 

Fs=Fs*Conpress 

, 

PRINT “CONVERTING THIS DATA TO METRIC UNITS (mm.N.,KPas fot) 0S. 9) (ee 


PRINT “WIDTH =";h:" NOTCH LENGTH =";A 

PRINT “Bmax =":Bmax:" Bmin =";Bmin 

PRINT "HALF SPAN =";Hs;” POISSON’S RATIO =";Nu 
ie) «66 Y CED STRENGTH =~"; Ys: UTILMATE TS =";:Uts 
PRINT “ELASTIC MODULUS =";5€;"° FLOK STRESS =":Fs 


PRINT 

PRINT "DO YOU WISH TO HAVE A PRINT QUT OF THIS DATA?” 

INPUT GS 

IF GS="N"° OR G3S="NO" THEN 1540 

a 

PRINTER IS 706 * TOGGLES TO THERMAL PRINTER 

PRINT “CONVERTING THIS DATA TO METRIC UNITS <(mm.N.KPa) FOR: ";DS 
PRINT “WIDTH =":sW;" NOTCH LENGTH =";A 

PRINT “Bmax 2”;Bmax:” Bain =" :Bmin 


PRINT “HALF SPAN =";Hs:" POISSON'S RATIO =";Nu 
PRINT “YIELD STRENGTH =";Y¥s;" ULTIMATE TS =";Uts 
Pen “ELASTIC MODULUS =";E;" FLOW STRESS =";Fs 


PRINT '** 

PRINTER IS 1 ' TOGGLES PRINTER TO CRT 
9 

PRINT CHRS(12) ' CLEAR CRT 

i) 

RETURN ' RETURNS TO &-KEY SELECTION 

' 


’ WD D> Tb WD Tb 2D TOD TD OD 4 4 4D 4D 4 4D 10 1 4 OOO Oe ee 
’ 


Setup plotter: ! SET UP INTERACTIVE PLOTTER 


' YOU CAN CHOOSE BETWEEN THREE DIFFERENT PLOTS 

' (P vs COD. P vs DEFLECTION. OR J vs DELTA_A)D 

PRINT CHRS(12) ' CLEARS CRT 

DIM TS(801,XS(80),YS(80] 

GCLEAR 

GINIT 

GRAPHICS ON 

PLOTTER IS 3,° INTERNAL” ' PLOTS ON CRT FOR VIEWING BEFORE SENDING TO 


Mee? 225B PLOTTER 


1690 
1700 
1710 
1720 
1730 


mov.  0YPE OF PLOT:” 

PRINT * 1. LOAD vs COD” 

ea ae 2. LOAD vs DEFLECTION” 

PRINT " 3. J vs CRACK EXTENSION” 

pon ENTER THE 2 9Ure tHe PLOT YOU DESIRE” 
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1740 [INPUT Ce 

175G PRINT "” 

1760 PRINT “* 

1770 IF C2=! THEN 1810 

1780 IF C2=2 THEN 1860 

1790 IF C2=3 THEN [910 

1800 PRINT “” 

1810 TS="PLOT of LOAD vs COD for "&DS 

1820 PRINT TS 

1830 xXS="COD (mm)” 

1840 YS="LGAD (N)” 

1380 GOTG 19350 

1860 TS=""PLOT of LOAD vs DEFLECTION for “&DS 
1870 PRINT TS 

1880 XS="DEFLECTION (mm)" 

1890 YS$="LOAD (N)” 

13900 GOTO 1350 

1910 TS="PLOT of J vs CRACK EXTENSION for "4053 
1920 PRINT TS 

1930 XS="DELTA_a (mm)" 

1940 YS="J  (kJU/m~2)” 

1950 PRINT “" 

1960 PRINT "USING METRIC UNITS (mm,N.xJ/m'2,etc):” 
1970 PRINT “ENTER HERE THE Xmin,Xmax.Ymin. Ymax"’ 
1980 BNeUT  X1¢CC2),X2¢0C2),71 (C2) , Y2¢C2) 

1990) X3=X2(C2)-X1(C2) 

2000 Y3=Y2(C2)-Y1¢(C2) 

2010 VIEWPORT 13,.100*RATIO.10,90 

2020 WINDOW X1¢(C2),.X2(C2),Y1¢C2) , Y2¢C2) 

2030 PRINT ° 

2040 PRINT “INPUT DELTA_X AND DELTA_Y LABEL INTERVAL” 
2050 PRINT “FOR THE TWO AXES: 

2060 INPUT X4,Y4 

2070 PRINT “” 

2080 PRINT “DO YOU WISH- TO HAVE A GRID OVERLAY?” 
2030 INPUT GS 

2100 IF GS=""NO" OR GS="N” THEN 2140 

2110 PRINT 

2120 PRINT “INPUT GRID LINES PER DELTA LABEL:”™ 
2136 INPUT BO 

2140 PRINT “PRESS ‘ENTER’ TO PLOT” 

2150 INPUT HS 

2160 PRINT CHRS(12) ® CLEARS CRI <------------ 
2170 IF GS="NO" OR GS="N”" THEN 2310 

2180 FOR [=X4/B0+xX1(C2) TO X2¢(C2) STEP X4/B0 
21390 MOVE I,Y¥1(C2) 

2200 DRAW I,Y2(C2) 

2210 NEXT I 

2220 FOR I=Y4/B0+Y1(C2) TO Y2¢C2) STEP Y4/B0 
2230 MOVE X1(C2),I 

2240 DRAW X2(C2) ,I 

2250 NEXT I 

2260 GOTO 2320 

22/70 IF GS="Y" OR GS="YES* THEN 2320 

2280 MOVE X2(C2),Y1(C2) 

2290 DRAW X2(C2).Y2(C2) 

2300 DRAW X1(C2).Y2¢(C2) 

2310 CLIP ON 

2320 AXES X4,Y4,X1(C2),Y1(C2) 

2330 CLIP OFF 
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2340 
2350 
2260 
2370 
2330 
2390 
2400 
2410 
2420 


2850 
2860 
2870 
2880 
2330 
2900 
2910 
2920 
2330 


VIEWPORT 13.100*RATIV.S,.90 
WMONDGH X1(Cl).X2°C2),¥1(Ce), YoCCo) 
LORG 4 
CSEZE3..0 ' LABELS X DIGITS 
FOR X=X1(C2) TO X2(C2) STEP X4 
GeiP OFF 
MOVE X.Y1«(C2) 
Gale gre 
LABEL USING “K";X 
NEXT X 
VIEKPORT 4,.1/00*RATIO, 10,90 
BOENDOW X1¢(C2) ,X20C2) .Y10C2) , Y2(C2) 
LORG 2 ' LABELS Y DIGITS 
Bore reyi¢C2Z) Th ¥26C2) STEP Y4 
Cele OFF 
MOVE Xi(C2),Y 
ier OFF 
LABEL USING "K";Y 
NEXT Y 
VIEXPORT O.RATIO#100,0,100 
aunoun X10C2),X2(C2) , ¥1¢(C2) , Y2¢C2) 
GoIZe 5 
LORG 6 Pero loam TEE 
Mave X2(C2)/2,%2¢C2) 
Babely i> 
DEG 
LDIR 390 
SSIZE 4.7 
mie A1NCE2) , Y20C2)/2 
LABEL YS Pepibtss Y-AXIS LABEL 
LORG 4 
LDIR 0 
MOVE X2(C2)/2,¥1¢C2) 
LABEL XS ! PLOTS X-AXIS LABEL 
f 


PRINT "" 

PRINT "" 

PRINT "* 

PRINT "* 

PRINT 7 IS THE PLOT (AXES AND LABELS) OK?" 

PRINT " NOTE: IF YOU ANSKER YES OR Y™ 

PRINT “ PLOT WILL BE PRODUCED” 

PRINT " ON EXTERNAL PLOTTER.” 

INPUT FS 

IF FS="N" OR FS="NO" THEN GOSUB Setup_plotter 1 REDO PLOTTER SETUP 


GCLEAR * CLEARS PLOT FROM CRT 
GOSUB Ext_plotter |! SENDS PLOT TO PLOTTER 


PENUP 
MOVE 0.0 
? 


RETURN ! RETURNS TO k-KEY SELECTION OR R curve ROUTINE 


9 WD WD WD AD WD AD WD HD 4D WD 2D OD DD DD OD DD 2D DD 2D DD 1D DD OD 1 
’ 


Ext_plotter: ! DURPS PLOTTER@SETUP TO EXTERMAL HP-722SB PLOTTER 
‘ 
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aT4U 
2950 
2360 
2970 
2280 
230 
3000 
3010 
3020 
3030 
3040 
3050 
3060 
3070 


3530 


gem to Ose aP Gl 
VIEWPORY 13.700*°RATIG.10.90 
MDOW X1(C2) .X2¢C2) ,YIECO). YOCCO) 
IF GS="NO0" OR GS="N" THEN 3030 
FOR [=X4/B0+X1(C2) TO X2(C2) STEP X4/B80 
MOVE = 1,11(Ca> 
ORAW I,Y2(C2) 
NEXT I 
FOR [=Y4/BO+Y1(C2) TO Y2«C2) STEP Y4/B0 
MOVE X1(C2),I 
DRAW X2(C2),I 
NEXT [ 
GOTO 3129 
IF GS2="¥" OR GS="YES” THEN 3120 
MOVE X2(C2),Y1¢C2) 
DRAW X2(C2) , Y2¢C2) 
DRAW X1(C2) , Y2¢C2) 
CEEP ON 
AXES X4,Y4,K1(C2) , Y1¢C2) 
Cre OFF 
VIEWPORT 13.100"RATIO.5.90 
WINDOW X1¢(CO) .X2(C2), YI(C2) , Y2¢C2) 
LORG 4 


Goze 3..5 ’ SIZE OF DIGITS 
FOR X=xX1¢C2) TOD X20C2) Steer X4 

CLIP OFF 

MOVE X.Y1(C2) 

CLIP OFF 

LABEL USING “K*;X f LABELS X-AXIS 
NEXT X 


VIEKPORT 4,100*RATIO. 10.90 
WINDOW X1(C2) ,X2(C2) ,Y1(C2) , Y2(C2) 
LORG 2 
FOR Y=Y¥1(C2) TO Y2(C2) STEP Y4 
CLIP OFF 
MOVE X1(C2),Y 
CLIP OFF . 
LABEL USING "“K";Y ! LABELS Y-AXIS 
NEXT Y 
VIEWPORT 0.100*RATIO,O.100 
WINDOW X1(C2) ,X2(C2),Y1(C2) .Y2¢C2) 
CSIZE 5S ' TITLE SIZE 
LORG 6 
MOVE X2(C2)/2,Y2(C2) 
LABEL TS ' PRINTS TITLE 
DEG 
LDIR 90 
CSIZE 4,.7 
MOVE X1(C2) ,Y2(C2)/2 
LABEL YS ! PRINTS Y-AXIS TITLE 
LORG 4 
LDIR 0 
MOVE X2(C2)/2.Y1(C2) 
LABEL XS ! PRINTS X-AXIS TITLE 
PENUP 
HOVE 0,0 


RETURN ' t RETURNS TO Setup plotter 


.. 
’ Sb TD TD TD TD TD TD ID TE TD WS 4 OD TD 4D DD ED DD ED EE 0 DD TD WD 4D AD TD 4D 4 OD 4 ed 2 a 


iA 





sit 
3550 
3560 
3570 
3580 
3590 
3600 
3610 
3620 
3630 
3640 
3650 
3660 
3670 
3680 
3690 


. 


[Initiai_a; 
Bet f =Bmax-(Bmax-Bmin) 


D2=! 

PRINT 
PRINT 
PRINT 
ral 
PRINT 
PRINT 
PRINT 
PRINT 


INPUT H 
PeeHS= "N* OR HS="HO” THEN 5076 


PRINT 
PRINT 


-CRACKING.” 


3700 
3710 


PRINT 
PRINT 


PleswinAtes UNL TAL CRACKS LENGsH 

2/Bmax vor omoterrec IVE) FROM ASTM £8: 5-31 
P INTTISELZES COUNTER 

CHRS(12) URCECANS Cn! 

“BEFORE YOU CAN RUN THIS ROUTINE YOU MUST:” 

a Ven CER OPEC I MEN DATA” 

= 2. ENTER CONVERSION DATA” 

" 3.0.58) UP SPeO ITER] 

Se CREA eset or Lee 


"ARE YOU READY?" 

z 

' RETURNS TO KEY SELECTION SCREEN 
CHRS(12) ' CLEARS CRI 

“THIS ROUTINE WILL OBTAIN THE INITIAL CRACK LENGTH AFTER FATIGUE PRE 


"YOU MAY LOAD AND UNLOAD THE SPECIMEN AS OFTEN AS YOU LIKE -- BUT Y 


OU MUST REMAIN "; 


3729 
3730 
3740 
3750 
3760 
3770 
3780 
3730 
3800 
3810 
3320 
3830 
3840 
. 3850 
3360 
3870 
3880 
3£.30 
2900 


Bold 


3920 


4090 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
’ 


PRINT 


INPUT E 


PORTIGN OF THE CURVE!" 


"ON THE LINEAR 

BiOte: HESATS SHOULD BE IN STROKE” CONTROL 

- AT A VERY SLOW INVERTED RAMP FUNCTION’ 

is RATE OR YOU MAY DO THIS WITH MANUAL LOAD” 
Z CONTROL.” 


“WOULD YOU LIKE AN ESTIMATE OF Pmax (Lbs)?* 
5 


IF ES="N* GR E€S="NO" THEN 4060 ! SKIPS Pmax CALC 
' 


PRINT 
PRINT 
INPUT 


"YOU MUST PROVIDE AN ESTIMATE OF CURRENT CRACK 
“INPUT YOUR ESTIMATE IN INCHES.” 
Aest 


LENGTH” 


t CONVERT ALL DIMENSIONS BACK TO ENGLISH UNITS. 
' a0 IS SUFFIX ON ALL VARIABLES FOR THIS CALC. 
Bef fa0=BefFf/Conlen 

Na0=W/Conlen 


Blig=(Wa0-Aest) 


f ESTIMATED LENGTH OF UNBROKEN LIGAMENT IN INCHES 


Fsau=Fs/Conpress 9 
Hsa=Hs/Conlen 


Piab= 
PRINT 
PRINT 
PRINT 
' 


PRINTER IS 706 


.533"Beffal#Blig*2«Fsa0/(2*Hsad) 


"AN ESTIMATE OF Pmax [S";P1a0 
"Pmax ESTIMATES THE POINT AT WHICH THE CURVE” 
“BECOMES NONLINEAR, IAW ASTM E£813-81" 


' TOGGLES PRINTER 


PRINT ‘*” 
PRINT “GIVEN a=";Aest:"“inches THEN:” 


PRINT 
PRINT 


“ Pmax®";Pla0;"Lbs FROM ASTM E813-81" 


PRINTER IS 1 
! 
Ay Lod(8000) ,Def1¢8000) ,Cod(8000) 


! ASSIGN DEFL CORRECTION COEF USED IN Data_acq AND Tfake_data ROUTINES ; 
! THIS USES A 4th ORDER POLYNOMIAL FIT OF RECORDED DATA ON THE MTS SYSTEM 


S COMPLIANCE 


4100 


' 


ve 





stig 
a2 0 
4130 
4140 
4150 
4160 
4170 
4180 
4190 
4200 
4210 
4220 
4230 
4240 
4250 
4260 
4270 
4280 
42930 
4300 
4310 
4320 
4330 


Corrt!)#5./98E-35 
Corr(2)24,.7419E-6 
Corr(3)2-2.803E-9 
Corr(4)#1.471£-12 
Corr(S)2-2.7741E-16 
 ) 


PRINT "** 
PRINT “WHEN YOU ARE READY TO BEGIN PRESS THE ‘ENTER’ KEY." 
INPUT HS 


PRINT CHRS(12) f CLEARS CRT ; 
K=0 f USED AS A COUNTER 
T= * USED AS INDEX 


F lag! =0 ’ NOTE: FLAG! AND FLAG2 ARE RESET 
Flag2=9 ‘ INe bata aca AOUTINE 


WHILE Flagi=0 ' ESTABLISHES ONE LOAD-UNLOAD DATA TAKING LOOP 
GOSUB Data_acq 

END WHILE 

Flag2=9 ' RESETS FLAG2 


WHILE Flag2=0 
GOSUB Data_aca 
END WHILE 
PENUP 
* THIS SECTION CAN BE USED TO VIEW THE COLLECTED DATA BY REMOVING THE REM 


ARK CHARACTER (!). 


4340 
4350 
4360 
4370 
4380 
4390 
4400 
4410 
4420 


~ 4430 


4440 
4450 
4460 
4470 
4480 
4490 
4500 
4510 
4520 
4530 
4540 
4550 
4560 
4570 
4580 
4590 
4600 
4610 
4620 


mroR. t=! TO K 
' PRINT USING 3780; .[.Lod¢I) .Derl(I) .Cod¢TI) 
tf IMAGE X,3D,4X,6D.0,4X,30.7D,7X,3D.70 
' NEXT I 
' WAIT S 
F=f 
ho aK 
GOSUB Slope ' SLOPE = 1/COD COMPLIANCE 
Correl0(D2)=Correl ' SAVES CURRENT CORRELATION COEFFICIENT IN AN ARRAY 
' 


Comp10(D2)=Comel f SAVES CURRENT COMPLINACE VALUE IN AN ARRAY 

? 

GOSUG Crack_length ' USES COD COMPLIANCE 

e oC, 

AD(D2)=Crek i ' SAVE THE LAST a0 MEASUREMENT OBTAINED IN AN ARRAY 
9 


‘ FIND 95% CCNFIDENCE INTERVAL FOR a0 
Stiff=Stifft+T9=u4 
GOSUB Crack_length 
Amax0(D2) =Crek | 
StifF=Stiff-2"T9#U4 
GOSUB Crack length 
Amin0(D2)=Crck | 
j 
PRINT CHRS(12) ' CLEARS CRT 
ae “THE FOLLOWING INFORMATION HAS BEEN GBTAINED:” 
PRINT “INITIAL CRACK LENGTH ¢a0) =°°;3A0¢D2):"*mm FROM MEASUREME 


NT #°;D2 


4630 
4640 
mm” 

4650 
4660 
4670 


PRINT 
PRINT “WITH 95% CONDIFENCE INTERVAL ON a FROM";AminO(D2);" TO" ;Amax0(D2) 5" 


PRINT. > 


PRINT “USING COD COMPLIANCE OF :";Comel10¢02) :""mm/N” 
PRINT “WITH A CORRELATION OF :";CorrelO(D2);:"FROM":K;"DATA POINTS.” 


us 





4530 
4690 
4700 
4710 
4720 
4730 
4740 
4750 
4760 
4770 
4780 
4790 
4800 
4810 
4820 
4220 
4840 
4850 
4860 


5140 
5150 
5160 
5170 
5180 
5130 
5200 
a2 0 
5220 


0 
5240 


ee livl. . 

PRINT “DQ YOU WANT A HARD COPY?” 
INPUT HS 

IF HS="N" OR HS="NO" THEN 4830 

4 


PRINTER IS 706 * TOGGLES PRINTER 
PRINT ‘* 
PRINT “THE FOLLOWING INFORMATION HAS BEEN OBTAINED: ” 
PRINT “INITIAL CRACK LENGTH =":A0(D2):"mm FROM MEASUREMENT #02 
PRINT “WITH 95% CONFIDENCE INTERVAL FROM: Amin0(D2):" TO':Amax0¢D2) ; "mm" 
PRINT “USING COD COMPLIANCE OF :":Compl0(D2) :"mm/N” 
PRINT “NITH A CORRELATION OF:":CorreiQ(D2):"FROM’:K:"DATA POINTS.” 
N eeo0¢ 
PRINTER IS 1 ' TOGGLES PRINTER 
q 


b2-D2+1 ' INCREMENT COUNTER 
N oe ee 

PRINT “DO YOU WANT ANOTHER a0 MEASUREMENT?” 
INPUT GS 

IF G3b="Y" OR G3="YES" THEN 4060 

 ] 


0 PRINT “INPUT THE NUMBER (#) OF THE a0 MEASUREMENT YOU WISH 70 USE AND S 


Reon, NOTE: ONLY THIS SET OF DATA ARE” 
BrorNg STORED GN DISK.” 

INPUT D2 

’ 


—Al=A0¢D2) t ASSIGNS At TO DESIRED A0«D2) VALUE 


' NOTE: At IS USED IN Main_prog AS INITIAL CRACK LENGTH VALUE 
PRN) °° 

PRINT “YOU HAVE SELECTED MEASUREMENT #° ;D2 

PRINT * ARE YOU SURE?” 

INPUT PS 

NOTE: ONLY THIS aQ MEASUREMENT IS STORED ON DISK 


IF PS="N" OR PS="NO" THEN 4890 - ' ALLOWS CORRECTION 
K2=0 ' K2 IS UNLOADING COUNTER THIS WILL ENSURE CORRECT DAT 


STORED 


 GOSYB Stor_data !  ¢--------- see WRITES ON DISK 
PRINT CHRS(12) 
9 


! 
RETURN { RETURNS TO k-KEY SELECTION 


t 2 028 OH BOD DDR SHRED SOSH BDHDHHDDDD DDD 
4 


Data_acq: ! DATA ACQUISITION ROUTINE 


OUTPUT Dym;"VRS AFIT ALS’ 
' 


OUTPUT Dum:"AIt vTI" ' READS LOAD 
ENTER UGAelod¢ 1) 

Lod(L)=Lod(I)#S1(1) 

9 


OUTPUT Dvm;"AI2 VT1" ! READS COD 

ENTER Dvm;Cod( 1) 

Cod¢(1I)=(Cod¢1I))#SI(2) 

' NOTE: COD IS CALIBRATED ONLY FORM 0.475 TO 0.625 INCHES (MAX COD IS 0.1 


S INCHES) 
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5570 
5580 
9530 


5600 
5610 
5620 
5630 
5640 
5650 
E TO 
5660 
5670 
5680 
56390 
5700 
5710 
5720 
$730 
5740 
S750 
5760 
5770 
5780 
mm) 

5730 
5800 
$310 


Co 


SP 


' 

QUTPUT Dym:"AI3 VTt” Rens Deere TOn 
ENTER Dum: Der lc.) 

Def 1l¢(l)=Def l(I>#Ss1¢3) 

¢ 


GaSUB Correct_derl ' OETERMINE DEFL CORRECTION FACTOR 


Def 1(I)=Defl(L)-Totcor 
’ 


] 
! PLOT THE DATA POINT 
VIEWPORT 13.1090#RATIO.10.90 
Voie «1 (Ce) .X2¢C2) , ¥!¢C2) . Y26C2) 
IF C2=2 THEN 5410 

PLOT Cod¢I),Laod¢I) 
GOTO 5420 

PLOT Def1l¢I),Lod¢I) 


; 

! TAKE DATA FOR ONE COMPLETE LOADING AND UNLUADING 

Peeet) (HEN 5520 SKIES FEST ON THE FIRST DATA POINT 
t 


' TEST DATA 
Lodiff=Lod¢(I)-Lod([-1) 
IF Lodiff<-25.0 THEN Flagt=! *‘ MUST BE UNLOADING 


IF Loaiff>25.0 THEN Flag2=! ! MUST BE LGADING 
Lodiff=0. 

I-E+ ! INCREMENTS INDEX 

K=K+1 t INCREMENT COUNTER 

RETURN ! RETURNS TO Initial_a ROUTINE 

$ 50 20 20 20 40 40 2D 0-00 1D 20 40 40 HD 40 20 20 40 40 4D aD 40 40 2D 40 40 40 Ob 4D SD 40 40 40 40 40 OD ob 0 wD DO wD 


rrect_derl: THIS ROUTINE CALCULATES THE DEFLECTION CORRECTION FACT 


CORRECTION FACTOR [S REQUIRED TO ACCOUNT FOR SYSTEM COMPLIANCE 
A 4TH ORDER POLYNOMIAL FIT WAS PREVIOUSLY DETERMINED. 
4 
NOTE: &D 40 20 1D 40 50 40 50 40 40-40 40 40 1D 40 40 40 39 TD 10 20 
THE S CORRECTION COEFFICIENTS CORR(I) WERE ASSIGNED IN [nitial_a RUUTIN 
ED UP PROG. 


9B 2b Tb 0 20 10 10 ED 10 20 10 20 ODD BD 1D 1D 00 00 20 ED ED 2D 10 10 0 0 


<o om [T) + <0 o<w 0 eww ow 0 


Lodp=Lod¢I)/Conlb ' CONVERTS CURRENT LOAD TO Lbs FOR CALC 
Def 1c(1)=Lodp 
FOR J=2 TO 4 
Def le( J) =Def lco(J-1) Def le (1) 
NEXT J 
Totcor=Corr(1) ' PREVIOUSLY ASSIGNED AS NOTED. 
FOR J=1 TO 4 
Totcor=Totcort+tCorr(J+!) Def lel) 
NEXT J 
f 


Totcor=Toteor=Conlen ' CONVERTS TOTAL CORRECTION FACTOR FROM <in) TO ¢ 


: 
' 


RETURN ' RETURNS 70 Data_acq OR Take data ROUTINE 
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5820 
$830 
5840 


¢ 
’ SSeOee,ea eens asvtdesseweeenseananeanneeennvancdanainaa 


: 


S850 Crack_length: ! CALCULATES CRACK LENGTH USING HUDAK (DINSRDC PROG) (PROF. J 
peouyCce S EQUATION] 


5360 
NCe. 
5370 
5880 
$890 
53300 
53310 
3920 
39330 
3940 
5950 
5960 
3370 
5980 
2290 
6000 
6910 
6020 
6930 
6040 
6050 
6060 
6070 
6080 
6990 
6100 
6110 
6120 
6130 
6140 
6150 
6160 
6170 
6180 
6190 
6200 
6216 
6220 
6230 
6240 
6250 
6260 
6270 
6280 
6290 
6300 
6310 
6320 
6330 
6340 
6350 
6360 
6370 
6380 
6390 
TAGE 


* NOTE: THIS ROUTINE USES COD COMPLIANCE METHOD VICE DEFLECTION COMPLIA 


FOR I=t TO 5 
READ Coet(I) 
NEXT I 
Comol=1/Stif Ff 
C=2"WeRBef ft *Compl/Hs/1t000 
C2=1/(SQR(C)+1) 
Gt) =C 
FOR I=2 10 5 
Gil) =GCI-1)#GC1) 
NEXT I 
F=Coef(t) 
FOR [=1 TO 5 
F=F+Coef (I+!) G(T) 
NEXT I 
Crek l=F ew 
RESTORE 
DATA 0.998265.-3.81662,-1.80596,22.31028,-44.15665,-52.67876 


RETURN f RETURNS 70 Initial_a OR Main_prog 
: 


? DDD DE DHPPPPHPBSHPHPPHPPHPPBSHPPPDHHHDDPHDD Dee DD 
’ 


Conv_data: t INPUTS CONVERSION DATA FOR MEASURING DEVICES 


PRINT CHRS(12) £ CLEARS CRT 

PRINT “INPUT: LOAD RANGE. STRAIN RANGE, STROKE RANGE” 

et ia CENTER IN PERCENT ie. 10,20,50,100)" 
INPUT Lr,Cr,Sr 

PRINT °°" , 

PRINT “LOAD RANGE = “;Lr 

PRINT “STRAIN (COD) RANGE = “:Cr 

PRINT “STROKE (DEFLECTION) RANGE = "3Sr 

Lmax=20000 ! MAX LOAD RANGE IN (lbs). 

Cmax=.200 ! MAX COD RANGE IN (in). 

Smax=5.0 ! MAX STROKE RANGE IN (in). 

eal ! MAX VOLTAGE IN (volts). 

! MAXIMUM OPERATING POINTS ARE: 

Olmax=Lmax*Lr/100 

Gcmax=Cmax*Cr/100 

Osmax=Smax#Sr/100 

t 


f CONVERSION CONSTANTS ARE: (ENGLISH UNITS) 

Sle(1)=Olmax/V 

Sle(2)=Ocmax/V 

Sle(3)=Osmax/V 

‘ CONVERTING TO METRIC: 

Conlen=25.4 { 25.4mm = lin 

Conib=4 448221615 { 4.44822161SN = It1lb (HP PROGRAMING BOOK) 


' Sle(1,2,3) IS USED TO DENOTE CONVERSION OF VOLTAGE TO ENGLISH UNITS. 


$1¢(1)=Sle(1)#Conlb#¢-1.0) ' COMPRESSIVE LOAD IS NEGATIVE VOTLAGE 
S1(2)*Sle(2)*Conien 


$1(3)*Sle(3)*Conlens(-1.0) ! COMPRESSIVE STROKE SUPPLIES A NEGATIVE VOL 
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541) 
S4su 
6329 
6430 


6440 
6450 
6460 
6470 
6480 
6490 
6500 
6510 
6520 
6530 
6540 
6550 
6560 


6860 
6870 
6880 
6890 
63900 
6910 
6920 
6930 
6340 
63950 
6960 
6970 
6980 
6390 


] 
rode ok 
PRINT "* 
PRINT “YOUR ENTRIES HAVE PRODUCED THE FOLLOWING CONVERSION DATA: 
PRINT "" 
PRINT “FOR LOAD: “:SLCI):"NEWTON/VOLT c(*:Sle(t):"Lbe/sy)" 
PRINT "" 
PRINT “FOR COD: “:S1(2)3"mmsY (Sled?) "ans" 
PRINT “** 
a ON “FOR DEFLECTION: “:S1(3):"mm/V (":Sle(3):"in/V)" 
2 oee00 
PRINT “WOULD YOU LIKE A HARD COPY?" 
INPUT HS 
IF HS="N”° OR HS="NO" THEN 6630 
q 


PRINTER Ke WIS ! TOGGLES PRINTER 

NT ee ee 

PRINT * CONVERSION DATA” 

PRINT “FOR LOAD: ";SI(1):"N/V ("sSle(1) :"Lbs/V)" 

PRINT “FOR COD: ";:S1¢(2)3"mm/¥ Ce Slet2Z): in/V) 

PRINT “FOR DEFLECTION: °S102)2"mm/VC" SLet3) 2" in/V)" 


PRINTER: IS ! TOGGLES PRINTER 
' 

PRINT CHRS(12) 

’ 

RETURN ! RETURNS TO k-KEY SELECTION 
9 


o 
’ WW OD OD DOGG DH DDD 
] 


Slope: ! LINEAR LEAST SQUARES FIT 


PRINT CHRS(12) * CLEAR CRT 

t INITIALIZE VARTABLES 

P2=0 ! USED AS A COUNTER 
Sumx=0 * SUM OF X VALUES 

Sumx2=0 ' SUM OF X SQUARED VALUES 
ec * SUM OF Y VALUES 

Sumy 2=0 ' SUM OF Y SQUARED VALUES 
Sumxy #0 * SUM OF THE PRODUCT 


FOR T= K1 TO K3 ' Kft=INDEX OF FIRST, K3=INDEX OF LAST DATA POINT USED. 
' TYPE OF PLOT 
its C2=1 THEN 6840 ‘ LOAD vs COD 
IF C2=2 THEN 6880 { LOAD vs DEFL 
IF C2=3 THEN 6920 ' J vs DELTA_a 


X=Cod(I) 
Y=Lod(I) 
GOTO 6990 
? 


X=Def ICT) 
Y=Lod( I) 
GOTO 6990 


X=Delta_a(I) 
Y=J9(T) 
q 


' CHECK IF DATA POINT IS WITHIN BOUNDS 
IF Delta_a(I)<J9¢I)/(2*Fs)+.15 THEN 7050 ' CHECKS LOWER LIMIT 
IF Delta_a(I)>J9(I)/(2Fs)+1.5 THEN 7050 * CHECKS UPPER LIMIT 


Sumx=Sumx+X 


Afi 





7900 
7010 
7920 
7030 
7040 
7050 
7060 
7070 
7080 
7690 
7100 
7110 
7120 
7130 
7140 
7150 
7160 
7170 
71890 
7130 

200 
7210 
7220 

230 
7240 
720 
7260 
7270 
7280 
7290 
7300 
7310 
7320 
7330 
7340 
7350 
7360 
7370 
7380 
7390 
7400 
7410 
7420 
7430 
7440 
7450 
7460 
7470 


Sumy *Sumy +7 


Sumx2=Sumx2+X*X 

Sumy22Sumy2+Y¥ 

Sumxy =Sumxy +X 

GOTO 7060 * GET NEXT DATA POINT 

a ' INCREMENTS COUNTER EACH TIME A POINT IS ELIMINATED. 
NEX 
K3=K3-P2 ' USED TO CALCULATE NUMBER OF DATA POINTS WITHIN BOUNDS. 
K4=K3-K1 +] ' K4=TOTAL #« DATA POINTS USED IN CALC : 


Qal =Sumx2-Sumx*Sumx/K4 
QaZ=Sumxy-Sumx*Sumy /K4 
pa = suny co Uby =Sumy “A 


M=Qa2/Qal ' SLOPE 
? 


IF M<O THEN 7430 ' IF SLOPE < 0, BAD DATA 
' FROM DITNSRDC PROGRAM 
Q1=(K4#Sumx2-Sumx*2)/(K4e(K4-1)) 
Q2=(K4*#Sumy2-Sumy ~2) /(K4"(K4-1)) 
Q3=«(K4-1)#(Q2-M=M#QI1)/(K 4-2) 
Q4=SQR(Q3/(Q1#(K4-1))) 

' 


B=(Sumy-M«Sumx)/K4 t INTERCEPT 
? 

Correll =M#SQR(Q1I/Q2) 

] 


Qa4=SGRCABS( (Sumy2-B*Sumy-M*Sumxy)/(K4-2))) 
Sigmab=Qa4/SGRC(ABS(Qai > 
Sigmam=S:.gman*#SOR(Sumx2/K4) 

Stiff=m * STIFF = SLOPS 
Comp|1=1/M * COD COMPLIANCE = !/SLOPE 
) 


! STUDENT t TEST 

U4=Q4 

T3=1.96 

“a ae THEN 7410 ! RETURNS 
IF K4<20 THEN 7410 ! RETURNS 
N9=K4-2 

Beet: 617 103EN9~ 510239) 9 


RETURN 

’ 

PRINT “POOR DATA-- SLOPE < 0” 
M=333333.3 

Correl=0 

, 


RETURN ! RETURNS TO Initial_a, Main_prog or R_curve ROUTINE (WHICHE 


VER CALLED IT) 


7480 
7490 
7500 
7310 
7320 
7330 
7540 


? 
» 
? > Wh DD 20 OE 1S 0 OE 1S OS OS DD 4D 4D OE OD DD 4D 4D 4 DD DE 4D SD 1 OS 4 4 
] 


Data_file: t ESTABLISHES DATA FILE 


! THE FILE YOU CREATE CAN BE READ BY THE "“J_DATA” PROGRAM 
! 
! NOTE: ANY CHANGES MADE TO THIS ROUTINE MAY REQUIRE SIMILAR CHANGES TO 


J_DATA” PROGRAM DUE 


7520 


7560 
7570 


t TO DISK FORMATTING CONSULT HP MANUAL FOR ASSISTANCE. 
PRINT CHRS(12) ' CLAERS CRT 
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7580 


PRINT “THIS ROUTINE WILL CREATE A DATA FILE ON OISK FOR THE TEST.” 


AS 


7390 PRIN?; “°° 

7600 PRINT ° "INSERT AN INITIALIZED DATA DISK AND ENTER FILE NAME (SAMPLE I 
D iS SUGGESTED). 

7610 INPUT NameS 

7620 PRINT °** 

7630 PRINT “IS THIS A NEW FILE?” 

7640 INPUT 05 

7650 IF QS="N" OR QS="NO”" THEN 7680 

7660 ’ 

7670 CREATE BDAT NameS.6,40000 

7680 ASSIGN @File TO Names 

7690 : 

7700 PRINT ‘* 

oe ae “FOR RECORD PURPOSES THE FOLLOWING DATA [S ALSO RECORDED :” 
7730 PRINT "ENTER THE DATE (MM-DAY-YEAR)D:” 

774) INPUT DataS¢1) 

7750 PRINT °°" 

7769 PRINT “ENTER THE TIME (HH:MM):” 

7770 INPUT DataS(¢2) 

7780 PRINT ‘** 

77930 PRINT “ENTER THE OPERATOR'S NAME: 

7800 PRINT ° (UP TQ 10 CHARACTERS)” 

7310 INPUT DataS¢3) 

7820 PRINT ° 

7830 PRINT “ENTER SPECIMEN ID:" 

7840 PRON, © (UP TO 10 CHARACTERS)” 

7850 INPUT DataS¢4) 

7860 ' WRITE ON DISK 

7870 =68FOR [=1 TQ 4 

78380 OUTPUT @File;DataS(I)(1;10] 
*- 7890 NEXT I 

7900 ! 

7910 Beer 

7920 PRINT °° 

7930 PRIN]? “DATA FILE HAS BEEN CREATED.” 

7940 t 

7950 ! 

aap e RETURN * RETURNS TO «-KEY SELECTICON 

(380 9 20 20 00 90 20 20 30 30 20 20 20 IO ID OD 10 10 IO OO-TE TO TE Te Hh 2 AD UE UE OD WD 1D 0 4D AD AO OD ab 2O 20D ODD 

7990 ! 

ape Stor_data: ! WRITES DATA ON DISK 

1 ! ' 

8020 ' THE ITEMS STORED BY THIS ROUTINE CAN BE 

8030 ! RECALLED USING THE “J_DATA*’ PROGRAM 

3040 ! CARE MUST BE TAKEN IF ANY CHANGES ARE MADE TO THIS ROUTINE SUCH THAT SI 
MILAR CHAGES ARE 

eh MADE TO THE “J_DATA” PROGRAM. 

3070 IF K2>0 THEN 8160 ' K2 IS UNLOADING COUNTER 
oeeG OUTPUT @File;A0(D2) ,AminO¢(D2) ,Amax0(D2) ,Comp10(D2) ,Correl0(D2) ,K 
8100 PRINT °** 

8110 PRINT **** 

a PRINT “INITIAL CRACK LENGTH ¢a0) INFORMATION FROM MEASUREMENT 2"; 
pe PRINT “IS NOW IN DISK FILE “;NameS 

149 ! 





3150 IF K22)0 SHEN 3220 ' SKIP [F THIS WAS ad MEASUREMENT 


3160 OUTPUT sFile:K2,AlL(K2) Delta_a(K2) ,K.L.Correl.Compl .Amin.Amax .AreaQ(K2 
») JO(K2) 

3170 OUTPUT 3File:Mdef 1(K2) ,Codiff,K0.K! 

8180 ! 


3190 PRINT “™ 
3200 PRINT “UNLOADING INFORMATICN IS NOW IN DISK FILE ":Names 
8210! 


3220 BEEP ; 
Bene ! WAIT 3 ! ALLOKS A PAUSc SO MESSAGE IS PRESENTED ON THE CRT 
<4 ! 
250 ! 
Beoe RETURN ' RETURNS TO Initial_a OR Main_prog 
al : 
8280 Fae 50 50 10 4b 20 40 4b 1b HD 40 40 A 40 10 40 Sb 10 10 10 ab 4 SD 40 0 40 4 SE 40 10 10 10 40 10 40 20 10 10 40 OD 2 40 
3290 ! 
ln Take_data: ! THIS ROUTINE IS USED TQ TAKE DATA DURING THE ACTUAL TEST. 
4 1 ' 
aio 2 [=« ' INDEX 
333 OUTPUT Dym:"AI1 VT1” ! READS LOAD 
3340 ENTER Dvm:Lod(I) 
3350 Lod(I)=Lod(I)#S1(1) 
3360 ! 
8370 OUTPUT Dym:"AI2 VIt" ' READS COD 
8380 ENTER Dvm:Cod(i) 
8330 Cod(I) =Cod(I)#S1(2) 
8400 ! 
8410 OUTPUT Dym:"AI3 YTt” ' READS DEFLECTION 
8420 ENTER Dum: Def 1¢I) 
te Def l<fo=Def ll) *Sit3) 
4 ’ 
ae GOSUB Correct_def 1 ! DETERMINE CORRECTION FACTOR FOR SYSTcM COMPLIA 
; GS 
8460 ! 
3470 Def l¢(I)=Defl¢I)-Totcor 
8430 ’ 


8490 4 PLOT THE DATA POINT 

8500 VIEWPORT 13,100«RATIO,10,90 

3510 ° WINDOW X1(C2) ,.X2°C2),¥1¢C2) , Y2¢0C2) 
3520 IF C2*2 THEN 8550 


6530 PLOT Cod(¢I) ,Lod¢I) 

3540 GOTO 8560 

3550 PLOT Def l¢I) ,Lod(I) 

8560 ! 

8570 ! 

8580 RETURN * RETURNS TG Main_prog 

85390 ! 

8600 eo ee ee ee ee eee ere 
8610 ! 


ooaey * CALCULATES J WITH CRACK GROWTH CORRECTION ERNST MODEL (USED BY D 
8630 ! 

8640 ES=2 

3650 ' CALCULATE J FOR FIRST (Ist) UNLOADING 

8660 IF K2>1 THEN 8690 

8670 J62J7+ES/CW-ALCK2) ) *AreaS(K2)/Bmin 

8680 GOTO 8730 

3690 JS2J7+ES/CW-ALCK2) )*(AreaS (KZ) -Aread(K2-1))/Bmin 

3700 ' IF DELTA_a < 0.254mm DO NOT USE CRACK GROWTH CORRECTION 

3710 IF Delta_a(K2)<.254 THEN 8730 
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8720 
Sit 631 
8730 
8740 
8730 
8760 
8770 
8780 
37390 
8800 


J6*J62(1-1/CW-ALCK2)*CALCK2) -ALCK2-1))) t CORRECTS rOR CRACK GROWTH (A 


3-81) 

J7=J6 

A] 

JS(K2)=J6 t STORES CURRENT J VALUE IN J9 ARRAY. 
t 

RETURN t RETURNS TO Main_prog 

’ 


, SEM He IE Ee I TE TE TE Ee De TE EE EE IE EE EE IE EE 2 
’ 


8810 Main_prog: ! THIS IS THE MAIN PROGRAM FOR’ THE J-INTEGRAL TEST. 
’ 


8820 
8830 
8840 
8850 
EDING:” 
3860 
8870 
8380 
8890 
83900 
8310 
83920 
83930 
8940 


8980 
8990 
9000 
OKE’ CQ 
9010 
9020 
3030 
9040 
TS MACH 
9050 
3060 
3070 
3080 
30390 
9100 
39110 
9120 
93130 
9140 
9150 
9160 
9170 
9180 
UNLOADI 
9190 


3200 

<10 
9220 
9230 
9240 
9250 


PRINT CHRS(12) * CLEARS CRT 
R N eee 
PRINT “ENSURE YOU HAVE COMPLETED THE FOLLOWING ITEMS BEFORE PROCE 


PRINT “ 1. ENTERED SPECIMEN DATA.” 
PRINT " 2. ENTERED CONVERSION DATA.” 


PRINT " Cer oe Ure Invencl IVE PLOTTER.” 

Nt : 4, OBTAINED INITIAL CRACK LENGTH.” 
N tees 

PRINT * ARE YOU READY?” 

INPUT HS 


Sa OR HS="N" THEN 10530 ! RETURNS TQ WAIT FOR k-KEY SELECTION 


PRINT °°” 

PRINT “WHAT IS YOUR EXPECTED MAXIMUM LOAD IN NEWTONS?” 

INPUT Pmax : 

LodS=.1*Pmax * USED TQ SIGNAL TURN ARGUND FROM UNLOADING 
7 


PRINT “THE TEST SHOULD START WITH NEAR ZERO (0) LOAD WITH THE MTS IN “STR 


NTROL.” 

PRINT °°" 

PRINT "ENSURE THIS IS DONE.” 
PRINT eeoe 


PRINT “WHEN YOU ARE READY PRESS ‘ENTER’ OR ‘CONTINUE’ KEY AND START THE M 
INPUT HS 
PRINT CHRS(12) t CLEARS CRT 


¢ 


' INITIALIZE VARIABLES 
Area=0 t AREA UNDER LOAD vs DEFLECTION CURVE  DTNSRDC (18) 
Area9<1)=0 ' AREA BENEATH SPECIFIC LOD vs DEFL (DTNSRDC I9(K2)) 


K=K ’ CURRENT INDEX (AFTER OBTAINING INITIAL CRACK LENGTH) 
T=K t INDEX OF FIRST TEST DATA POINT 

KO=K ! WHERE TEST DATA STARTS 

K1=0 { LOWER INDEX FOR SLOPE CALC 

K2=0 t A COUNTER - (K2) NUMBER OF UNLOADINGS 

K3=0 ! UPPER INDEX FOR SLOPE CALC 

Ne ' A COUNTER ODTNSRDC (L) <r" DATA POINTS TAKEN DURING 
Lod!=0 ! TEMP LOAD . (P1) USED IN AREA CALC 

Lod2=0 ! TEMP LOAD ‘si (P2) = : aes 
LodS=.1*Pmax ! : (PS) USED FOR TURN AROUND SIGNAL 
Lod3=0 ! MAX LOAD . (P93) 

Cod1= lent cop (DID 

Cod2=0 ' TEMP COD i (D2) 

Cod4=0 ! TEMP COD . (C4) 
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3260 


9830 


Cod9=9 ! MAX COD 7 (D9) 

Codi 20)? CODRDIFFERENCE ™ CUS) 

De+11=0 Sehr CRF es (D4) USED IN AREA CALC 
Def 12=0 J eal tam (DS) a i 
J6=0 ft JI FROM CURRENT UNLOAD 

J7=0 ! JI FROM PREVIOUS UNLOAD DTNSDRC ¢J7) 

I1=0 ! FLAG (Q=LOADING, 1t=UNLOADING) 

a Y FLAG (Fr te9S9 ENDS TEST) 

! 

PRIN 

? 


IF F1=999 THEN GOSUB End_test 
? 


GOSUB Take data 
Def 16=Def1l¢I)D 
Lod/=Lod¢I) 
Lod2=Lod( I) 

Def 12=Def lI) 
Cod2=Cod (1) 


TEST CONDITIONS OF DATA POINT JUST TAKEN 
IF (K-KO)<20 THEN 10420 f CALC AREA & TAKE MORE DATA 
Codif f=Cod2-Cod! 


IF Codiff<-3.0E-4 THEN 9540 f IF CODIFF IS ¢-) CHECK IF LOADING OR UNLO 
IF I1=0 THEN 10430 f IF LOADING CALC AREA & TAKE MORE DATA 

IF Cod2>Cod3 THEN 3608 ' IF TRUE CHECK IF GOOD UNLOADING 

GOTO 10340 DISREGARD UPPER 2% OF UNLOADING CURVE 

IF 11<>0 THEN 10340 IF UNLOADING CHECK IF ON UPPER 2% OF CURVE 
Lod3=Lod2 f IF LOADING MAKE CURRENT LOAD = MAXIMUM LOAD 
Cod9=Cod2 f MAKE CURRENT COD = MAX COD 

Def 19=Defl2 ' MAKE CURRENT DEFL = MAX DEFL 

[1=1 f MUST BE UNLOADING 

GOTO 10340 t TEST IF UPPER 2% OF UNLOADING 

IF L<S THEN 10290 ! BAD UNLOADING, SET I1t=1 

Kt=K-L ! SETS LOWER INDEX FOR SLOPE CALCULATION. 
K3=K1+L ! SETS UPPER INDEX FOR SLOPE CALCULATION. 

PENUP 

! 

——- ! INCREMENT UNLOADING COUNTER 

GOSUB Slope f IF GOGD UNLOADING 


Stiff2=Stif f 

GOSUB Crack_length ! CALC CRACK LENGTH 

Al (K2)=Crcek l ! CURRENT CRACK LENGTH 
Delta_a(K2)=Al(K2)-Al * DETERMINE CRACK GROWTH 


‘ NOTE: Al = INITIAL CRACK LENGTH; ASSIGNED IN Initial_a ROUTINE. 
' 


Crek 14=Crckl 

PRINT " ---- UNLOADING 7":K2:" ---” 
TPR ENT.***’ 

PRINT "CRACK LENGTH =" FALCK2) 5mm" 

PRINT "“DELTA_a =";Delta_a(K2) ;"’mm' 

PRINT “TOTAL ¢# DATA POINTS = "3K 

PRINT “# DATA POINTS THIS ilaapane ao 
PRINT "CORRELATION =";Correl 

PRINT “COD COMPLIANCE =";Compl1:"mm/N" 


82 





10020 
10030 
10040 
10050 
10060 
10070 
10080 
10090 
10100 
10110 
10120 


 ] 
AreaS(K2)=Area 

? 
* CALCULATE 95% CONFIDENCE INTERVAL 

Stiff=Stiff+T9*U4 

GOSUB Crack_length 

Amax=Crcek | 

StiffF=Stif F-2*T9*U4 

GOSUB Crack _length 

Stiff=Stiff2 

Amin=Crcek l 

PRINT "™* 

PRINT "95% CONFIDENCE INTERVAL ON CRACK LENGTH FROM “;Amin; "TO" ;Amax: 


GOSUB J.cale ! CALC J (USES J-ERNST MODEL) 
Ndef 1 (K2)=Def 1 (Kk) ' MAX DEFLECTION 


PRINT '" 

PRINT "AREA ="sAreaS(K2) ;"N*mm” 

PRINT "J ="3J9¢€K2) 3:"KJ/m*2" 

PRINT “LOAD LINE DEFLECTION ="; sMdef 1(K2) "mm" 
PRINT "CHANGE IN COD ="; sCodiff;" mm'' 

PRINT “KO#"sKO;:" Kt1=": Kt 

’ 

i 
GOSUB * PRINTS RESULTS ON PRINTER 


! DATA POINTS: CAN BE PRESENTED FOR REVIEW BY SIMPLY REMOVING THE REMARK C 


Prt_results 


HARACTER (!). 


10130 
10140 
10150 
10160 
10170 
10180 
10190 
10200 
10210 
10220 
10230 
10240 
10250 
10260 
10270 
10280 
102930 
10300 
10310 


K 
10320 
10330 
10340 


RE DATA 


10350 
10360 
10370 
10380 
16390 
10400 
10410 
10420 


PRINT °°" 


! PRINT “ PT # LOAD(N)D DEFLECTION(mm) 
' FOR I=(Kt-1) TO K 


£ PRINT USING 9170:I,Lod¢I) ,Def1l¢I) ,Cod¢I) 
! IMAGE 4D,4X,5D.D,6X,DD.6D.7X,DD.6D 


PENEAT 1 
f PAUSE TO ALLOW VIEWING 


COD<mm)" 


WAIT 5 
GOSUB Stor_data 


t 


IF Crcek14-A1>0 THEN 10290 
Pine 

PRINT " NGO CRACK EXTENSION THIS POINT.” 
PRINT "* 

I1=0 

ar 

! IF t=? THEN Kt#K 


GOTO 10360 
IF Cod2<Cod1 THEN 10440 f GET MORE DATA 
IF Lod2>.98*Lod93 THEN 10440 ! DISREGARD UPPER 2% OF UNLOADING & GET MO 


=e f INCREMENT DATA POINT COUNTER 
IF Lod2>Lod9-LodS THEN 10440 ! IF TRUE GET MORE DATA 
ae ! SENDS GUT BEEP 


1 AS TURN AROUND 
BEEP 
GOTO 10440 


' STORES INFORMATION ON DISK 


{ SETS LOWER INDEX FOR SLOPE CALCULATION AT CURRENT 


f SIGNAL TO BEGIN RELOADING 
f RESUME TAKING AND TEST DATA 


CALCULATE AREA UNDER LOAD vs DEFLECTION CURVE 
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10430 
10440 
10450 
10460 
10470 
10480 
10490 
10500 
10510 
10520 
10530 
10540 
10550 
10560 
10570 
10580 
10590 
10600 
10610 
10620 
10630 
10640 
10650 
10660 
10670 
10680 
10690 
10700 
10710 
10720 
10730 
10740 
10750 
10760 
10770 
10780 
10790 
10800 
10810 
10820 
10830 
10840 
10850 
10860 
10870 
10880 
10890 
10900 
10910 
10920 
10930 
103940 
10950 
10960 
103970 
10980 
es l0 
11000 
T1010 
11020 


Area=Areart (Def l2-Defll)>*(Lod2+Lod1)/2 


Lodt=Lod2 

K=K+1 ! INCREMENT TOTAL #* DATA POINTS COUNTER 

Cod1=Cod2 

Derive te 

GoTo 9380 ’ LOOP BACK CHECK IF TEST IS 

! OVER, IF NOT TAKE MORE DATA. 

! 

! 

RETURN ’ NEVER REACH THIS BUT ROUTINE WOULD NOT BE VALID WITHOUT IT 


’ FOE EE BG EE IE TS BE IO BE BE 06 TE 1G GE EE 8 EE IE EO 26 IE 26 3 EF 2 OE 


Prt_results: 
PRINTER IS 706— 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
' 


PRIN 
PRINTE 
° 
RETURN 


! PRINTS TEST INFORMATION WHILE PROGRAM IS RUNNING 
R ! TOGGLES PRINTER . 


s ---- UNLOADING #";K2;" ----" 


"CRACK LENGTH =":AL(K2) ;"mm DELTA _a =""3Delta_a(K2) ;“*mm" 
"95% CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM"’;Amin;” TO";Amax 


"TOTAL * DATA POINTS =";kK;" # DATA POINTS THIS UNLOADING ="3L 
“CORRELATION =";Correl 

"COD COMPLIANCE =";Comp! ;"mm/N” 

"AREA =":AreaQ9(K2) 3°°"N*mm” 

eS) et UOCK29 KS /m° 2" 

"MAX LOAD LINE DEFLECTION ="3;Mdef1l(K2) 3mm” 


"CHANGE IN COD =""sCodiff;"mm” 
RIS 1 ! TOGGLES PRINTER TO CRT 
! RETURNS TO Main_prog 


i Ce ee ee el 


' 


Stor_arrays: ' STORES ARRAYS (Lod,Defl,Cod) ON DISK 


! THE ARRRAYS CAN BE RECALLED LATER USING 


Ee [int 
! 


"J_DATA” PROGRAM. 


OUTPUT @File;Lod(*) ,END 
OUTPUT @File;Defl(*=),END 
OUTPUT @File;Cod(*) ,END 


' 
BEEP 
! 


? 


RETURN 


! RETURNS TO End_test 


’ Ce ee ee ee ee ee 


End_test: 
' 


PRINT 


’ 


GOSUB 


1 ENDS EST 
CHRS¢ 12) 


Stor_arrays 
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COWOOnNAMSAONM—-Ow 
COOCOOCOCOOoCoOoO 


PRINT “DATA HAS BEEN STORED ON DISK IN FILE “;NameS 
! 

GOSUB R_curve 

' 


PRINT "" 
PRINT ° 
PRINT "* 

PRINT “" 

PRINT """ 

PRINT “* 

PRINT "" 

PRINT “ TEST IS OVER!" 
PRINT "" 

PRINT "" 

PRINT """ 

PRINT “" 

PRINT "" 

PRINT " TEST IS OVER!” 


LES?) loedVent! © 


! 
GOTO 12590 ' GO TO END STATEMENT 
' 


RETURN ! NEVER REACH THIS POINT BUT ROUTINE WOULD NOT BE VALID WITHOUT 


! 

t Pe TE TE TE EE HE Ee IE Ee EE I EE EE IE I EE EE I EE I EE AE 0 EE I EE 8 

! : 

R_curve: ' J vs DELTA_a DATA IS PRESENTED FOR PLOTTING 


PRINT “" 
if 


f ALL (*) STATEMENTS WILL ALLOW DIRECT INPUT OF Delta_a(I) vs J9¢I) DATA 
{ YOU MUST REMOVE [(f¢(=)] FROM BEGINNING OF EACH LINE TO ACTIVATE 
'(*) PRINT “INPUT THE # OF UNLOADINGS” 
'(*) INPUT K2 
'(*+) PRINT "ENTER DELTA_ACI),J9¢ID" 
PRINT “UNLOADING # DELTA_a(mm) JS(k J/m*2)" 
FOR IT=1 TQ K2 
{¢*) INPUT Delta_a(I),J9¢(T) 
PRINT USING 11420;:I,Delta_a(TI) ,J9¢I) 
IMAGE 6X,3D,4X,DD.6D,5X,6D.D 


NEXT I 

BRUNI "7 

PRINT “DO YOU WANT A HARD COPY?” 

INPUT HS 

IF HS="N" OR HS="NO" THEN 11680 

PRINTER IS 706 ' TOGGLES TO THERMAL PRINTER 


PRINT “UNLOADING « DELTA_a(mm) JSCk J/m>2)"" 
FOR I=1 TO K2 
PRINT USING 11420;1,Delta_a¢I) ,J9(T) 
NEXT I 
PRINTER iS-t ! RETURNS TO CRT AS PRINTER 


'¢*) PRINT “ARE ALL OKAY?" 

'¢*) INPUT GS 

'c=) IF GS="Y" THEN 11582 

'(*) PRINT “INPUT I CINDEX) OF BAD #” 

'C#>) INPUT Ll 

'¢#) PRINT “INPUT VALUES OF DELTA_ACI) AND JSC)" 


85 





730 


12190 


'¢#) INPUT Delta_atLl),J9<LID 

'€*) PRINT "ARE YOU OK" 

'(*) INPUT EeeS 

'(») IF EeeS="N" THEN GOTO 11574 

' END OF (*) CHANGES REQUIRED FOR DIRECT INPUT OF DELTA_ACI>) vs J9¢I) DAT 


' 
PRINT *"" 
PRINT “DO YOU WANT A PLOT OF J vs DELTA_a?” 
INPUT RS 
vee OR RS="NO" THEN 12550 
N oeee 
PRINT “WHEN THE PLOTTER IS READY (NEW GRAPH PAPER aN OREEED BRESS ENTER 


OR “CONTINUE’ KEY" 


INPUT CS 
! 


GOSUB Setup_plotter 
) 


' PLOT THE DATA THE DATA POINTS 
VIEWPORT 13,100*RATIO,10,90 
WINDOW X1(C2) .X2(C2) , Y1 (C2) , Y2¢C2]) 
FOR I=1 TO K2 

X=Delta_a(I) 


YeJ9CT) 
MOVE X-X3*.005,Y+Y3*.005 ! DRAWS A 
DRAW X+X3*.005,Y-Y3*.005 ! SMALL X 
MOVE X+X3*.005,Y+Y3*.005 * FOR DATA 
DRAW X-X3*.005,Y-Y3*.005 ! POINT 
NEXT I 
’ 
! PLOT THE.»BLUNTING LINE 
VIEWPORT 13,100*RATIO,10,90 
WINDOW X1(C2) ,X2¢(C2),Y1(C2). Y2¢C2) 
Bly=Y2(C2)-Y3/8 ’ TAKES LINE TO 7/8 OF Ymax 
Blx=Bly/(2*Fs)*1000 ! DELTA_a = J/2*Fs 
MOVE 0,0 


Blminx=Blx+.15 
MOVE .15,0 

LINE TYPE 4 
DRAW Blminx,Bly 
’ 


Blmaxx=Blx+!.5 


' DRAWS LOWER BOUNDS 


! DRAWS UPPER BOUNDS 


MOVE 1.5,0 

LINE TYPE 4 

DRAW Blmaxx,Bly 

PENUP 
K1=1 f LOWER INDEX FOR LEAST SQUARES FIT 
K3=K2 f UPPER INDEX 


GOSUB Slope 
* PLOT THE LINE OBTAINED 
VIEWPORT 13,100*RATIO,10,90 
WINDOW X1¢(C2) ,X2¢C2) ,Y1¢C2) , Y2¢C2) 
Xp lt=X2(C2)-X3/8 ' TAKES LINE TO 7/8 OF Xmax 
Yplt=Xplt*M+B 
MOVE 0,B 
LINE TYPE 6 
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12200 
12210 
12220 
12230 
12240 
12250 
12260 
12270 
12280 
12290 
12300 
12310 
12320 
12330 
12340 
12350 
12360 
12370 
12380 
12390 
12400 
12410 
12420 
12430 
12440 
12450 
12460 
12470 
12480 
12490 
12500 
12510 
12520 
12520 
12540 
12550 
12560 
12570 
12580 
12590 
12600 
12610 
12620 
12630 
12640 
12650 


ROGRAM 


12660 


DRAW Xplt,Yplt 


! k9-KEY IS NOT LABELED, BUT IS ASSIGNED A HIGHER PRIORITY TI INTERRUPT P 


! 


PENuP 
PRINT °"* 
PRINT "“GNLY DATA WITHIN BOUNDS IS USED" 
PRINT "FOR LEAST SQUARES FIT.” 
PRINT "" 
PRINT “SLOPE=";M;" INTERCEPT=";B 
PRINT “CORREL=":Correl 
PRINT "NUMBER OF DATA POINTS USED:";K4 
PRINTER US 706 ' TOGGLES PRINTER 
PRINT 
PRINT ‘'*" 
PRINT “DATA FROM J vs DELTA_a PLOT” 

“PRINT "" 

PRINT "SLOPE =": ;" Y-INTERCEPT =";B 
PRINT “CORRELATION =";Correl;” 

PRINTER IS 1 ! TOGGLES PRINTER 
’ 

' FIND THE INTERSECTION. HENCE JCI1c). 

' SOLVE BOTH EQUATIONS SIMULTANEOUSLY 

! J = M*(DELTA _a)+B = (DELTA _a)*(2"Fs) 
J1c=(B/M)/¢€1/M-(1000/(2*Fs))) 

PRINT “" 

PRINT " iC w"e SUC 

PRINTER IS 706 ! TOGGLES PRINTER 
PRINT ** : 

PRINT * Jie ="3Jic 

PRINT "" 

PRINT °"* 

PRINT “ TEST IS OVER!” 

PRINTER fst ' TOGGLES PRINTER 
GOTO 12100 ‘ GO TO END STATEMENT 
pean f RETURNS TO End_test ROUTINE 
$ 3 36 38 3G 3G 6 6 36 6 6 36 2G 8 6 EF 6 EF 28 0 6 36 6 6 6 6 6 38 6 8 8 6 8 36 6 8 8 3G 6 8 8 8 
! 

END ' END OF TEST AND PROGRAM 
t DE HE IE GE HG GE BG OE GE ES EG OS GE EG 3G OG GG GG 0 0 0 St 0G 8 Et 
' 

SUB Stop_test feonORS TEST 


k9-KEY ACTS AS PROGRAM INTERRUPT 


12670 COM /J_intgl/ FI 


j2os0 Ft=999 


12690 


! 


a SUBEND 


$ 9G OE HG 3G HG HG HG EG 3G HG HG HG 36 ab HG HG 1G 2b 2 0G EE 0G 1G OE 1G Ab ED 1 GG 1 36 8 1G 36 1G 


2710 


Poe oo hUP SELAG 
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GJ) 
Ws) 
io) 

1 


APPENDIX F 


J DATA LISTING 


J_DATA PROGRAM 
THIS PROGRAM WILL ALLOW YOU TM READ A DATA 
FILE CREATED BY “J_INTGL’. YOU MUST KNOW 
HOW MANY UNLOADINGS WERE RECORDED DURING 
THE TEST TO USE THIS PROGRAM. 


PROGRAM WAS WRITTEN BY: 
LCDR W.K. TRITCHLER 
VERSION 1.0 DATE 11/83 


SD em SD ee OD om OR omm ' RB ow tw 


DIM DataS(1:4)(10] 

DIM Lod(8000) .Def 1¢8000) ,Cod(3000) 
DIM ALC 100) .Area9( 100) ,Delta_a¢100) 
DIM Mdef1¢(100),49¢(100) 

, 


PRINT CHRS(12) * CLEARS CRI 
9 


DUMP DEVICE IS 706 
PRINT “WHAT DATA FILE DO YOU WISH TO READ?” 
INPUT Names 
ASSIGN sFile TO NameS 
FOR I=1 T0 4 
ENTER SFile;:DataS(I)(1:10] 

NEXT [ 
PRINT "**° 
ae " J- INTEGRAL TEST OF “:DataS(4);"BY: *;Datas(3>d 
PRINT “DATE: “;DataS¢!);” TIME: “;DataS¢2) 

ENTER @File;Crcek!,Amin.Amax,.Comol,Correl ,K 
PRINT °°" 
PRINT °°” 
PRINT ** 
aT “a0="5Crcekl, “WITH 95% CONDIFENCE INTERVAL FROM’;Amin;” TO” ;Amax 


oat "FROM COMPLIANCE = “;Comp! ve 

RINT ee ee 

eT “CORRELATION=";Correl; "FROM" 3K; “DATA POINTS” 

PRINT "DO YOU KANT A HARD COPY OF THIS INFORMATION?” 

INPUT GS 

IF GS="N" OR GS="NO" THEN 450 

DUMP ALPHA 

PRINT "°° 

PRINT “INPUT THE NUMBER OF UNLOADINGS THAT WERE RECORDED DURING THE TEST.” 
INPUT C1 ‘ WILL BE USED AS COUNTER 


FOR I=! TO Ct 
ENTER @File;K2,A1(K2) ,Delta_a(K2) ,K,L,Correl,Compl ,Amin,Amax,Aread(K2), 


ENTER ®File:Mdef1l(K2) ,Codiff,KO.K! 

GOSUB Prt_results ’ DISPLAYS INFO 

PRINT ‘°° 

PRINT “DO YOU WANT A HARD COPY?” 

INPUT QS 

IF QS="N" OR QS="NO” THEN S80 

PRINTER IS 706 
nee Prt_results t WILL NOW PRINT INFOQ 
PRINT "DO YOU WANT A LISTING OF THE DATA POINTS?" 
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ES 
Tie GR E5=—N0" THEN 320 


oe ies LOAD (N) DEFLECTION (mm) COD (mm>" 
aay soe error Coats) 
| TO K 


NT USING 670:I.Lod(I) .Defi¢I),Cod(I) 


IMAGE X,40.3X,5D.D,8X,D.6D.3X,D.6D 


"DO YOU WANT A HARD COPY?” 
HS 
"NN" OR HS="NO" THEN 830 


PRINTER IS 706 


“DATA FROM a0 FOR: ";DataS(4);"BY: “;Datas(3) 


“TESTED ON: “sDataS¢1);"TIME: “;:DataS(2) 
“PT iz LOAD (ND) DEFLECTION (mm) COD (mm)" 
1 TO K 


NT USING 670;:I.Lod¢I) Def i¢I) .Cod¢ I) 


PRENTER IS | 


CHRS( 12) Choos CRI 

"DO YOU WISH TO PLOT ANY DATA?” 

ES 

"¥" OR ES="YES” THEN 880 ! OTHERWISE MUST NOW READ DATA POINTS OFF 


ow ile:Lod(#) Def l(*) ,Cod(*) 

"YY" OR PS="YES” THEN GOSUB Setup_plotter 
"D0 YOU WANT ANOTHER PLOT?’ 

AS 

"YY" OR AS="YES” THEN GOSUB Setup_plotter 
ate YOU WANT ANOTHER PLOT?” 

s e 

"¥" OR FS="YES" THEN GOSUB Setup_plotter 
PROGRAM IS OVER!!!" 


ae PRGGRAM IS OVER!” 
a PROGRAM IS OVER!" 
290 


, 3 2D 2D 2D 2b 2 2b 2b 2D TD 4b TD OS OD OD 2D OD OD OO OP OO HE OD OD 2S OP OD 2D DD OO 2D OD 2D OD Ob 2 2 OP OO 20 


Prt_results: ! DISPLAYS OR PRINTS UNLOADING INFORMATION 


600 INPUT 

6:0 ile ae 
620 PRINT 

650 PRINT 

640 ENTER 

fou FOR f= 
650 PRI 
670 

680 NEXT I 
690 PRINT 

700 =PRINT 

710 INPUT 

720 IF HS 
730 

740 = PRINT 

750 = =PRINT 

760 = PRINT 

770~—=—s PRINT 

7380 3 =6PRINT 

ZOU) FOR I= 
800 inset 
Bi0 NEXT I 
829 

830 PRINT 

849 PRINT 

850 INPUT 

860 IF ES= 
BIPSKET TE 

870 ENTER 

880 ie a= 
830 PRINT 

900 INPUT 

310 IF AS= 
220 PRINT 

920 INPUT 

940 IF FS= 
So) PRINT 

960 PRINT 

970 =6©PRINT 

980 PRINT 

Sou PRINT 

1000 PRINT 

1010 PRINT 

1020 PRINT 

1030 PRINT 

1040 PRINT 

1050 PRINT 

1060 PRINT 

1070 PRINT 

1030 =! 

1090 GOTO } 
1100 ! 

1110 

1i20 ! 

1130 

1149 

HisD0 PRINT 

1160 PRINT 

70) 60P RENT 

1180 PRINT 


“FROM UNLOADING #":K2;"RECORDED";L: "DATA POINTS” 

“COD COMPLIANCE =";Comel;:"mm/N —§ CORRELATION =";Correl 
“CRACK LENGTH =";3AL(K2);"mm DELTA_A =";Delta_a(K2) ; "mm" 
“9524 CONFIDENCE INTERVAL FROM’: Amin; "mm TO"; Amax: "mm" 
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1190 
1200 
t210 
1220 


1240 
1250 
1260 
1270 
1280 
12990 


1310 
1320 


PRINT “AREA =";Areagd(K2):"'Nemm Joe ISK) s KP aem 
PRINT “DATA POINT INDEXES FROM:K1; “TOKE WERE USED." 
PRIME =” 

PRONTER IS | 

° 

RETURN 

’ 

| ennnenaneennnnnnnnenneennenenssensnenannens 

’ 
Setup_plotter: Poe UR elle 

! YOU CAN CHOOSE BETWEEN THREE DIFFERENT PLOTS 

PRINT CHRS(12) PGREARS Chl 

DIM Ts(80)],XS(80].YS(80] 

GCLEAR 

GINIT 


GRAPHICS ON 
PLOTTER [S 3,° INTERNAL” PETS ONVeRT 
eee “TYPE OF PEGT:” 


PRINT “ 1. LOAD vs COD” 
RROUN @. LOAD vs DEFLECTION” 
PRINT “ 3. J vs CRACK EXTENSION” 


PRINT “ENTER THE ¢ OF THE PLOT YOQU DESIRE” 
INPUT C2 

PRINT “* 

PRINT “" 

IF C2=1 THEN 1470 

IF C222 THEN 1520 

IF C2=*3 THEN 1570 

TS="PLOT of LOAD vs COD for “&DataS(4) 

PRINT TS 

XS="COD (mm)” 

YS="LQAD «N)” 

GOTO 1610 

TS="PLOT of LOAD vs DEFLECTION for “&DataS(4) 
PRINT TS 

XS="DEFLECTION (mm) 

YS="LQAD «N)” 
GOTO 1610 
TS="PLOT OF J vs CRACK EXTENSION for “&DataS(4) 
PRINT TS 

X$S="DELTA_a (mm)"™ 

Ys="J3 00 (ke J/m*2)" 

PRINT “* 

PRINT “USING METRIC UNITS (mm,N.kJ/m°2.etc.):” 
PRINT “ENTER THE Xmin.Xmax,Ymin, Ymax” 

INPUT X1¢C2) ,X2¢(C2) ,Y1¢C2) , Y2¢C2) 
X3=X2(C2)-X1¢C2) 

Y3=Y2(C2)-Y1¢C2) 

VIEWPORT 13,100#RATIO,10,90 

eee ce arene) 

PRINT "INPUT DELTA_X AND DELTA_Y LABEL INTERVAL” 
PRINT “FOR THE TWO AXES:” 

INPUT X4,Y4 

PRINT "™" 

PRINT "DO YOU KISH TO HAVE A GRID OVERLAY?” 
INPUT GS 

IF GS="N" OR GS="NO"° THEN 1800 

PRINT ‘**° 

PRINT “INPUT GRID LINES PER DELTA LABEL:” 


°: 


90 





1790 
1800 
1810 
1820 
1830 
1340 
1850 
1860 
1870 
1880 
1890 
1900 
1910 
1320 
1930 


2260 
2270 
2280 
2290 
2300 
2310 
2320 
2330 
2340 
2350 
2360 
2370 
2380 


INPUT Bd 
PRINT “PRESS ‘ENTER’ OR “CONTINUE” KEY TQ PLOT™ 
INPUT HS ' DUMMY VARIABLE 


PRINT CHRS(12) ' CLEARS CRT 
IF GS="N° BOR GS="NO" THEN 1970 
FOR [=X4/B0+X1(C2) TQ X2¢«C2) STEP X4/B0 
Mev ly tices) 
DRAW I, ¥2¢(C2) 
NEXT I 
FOR [=¥4/B0+Y1I«C2) TO Y2¢C2) STEP Y4/B0 
MOVE X1¢C2).1 
DRAW X2(C2),1 
NEXT [I 
GOTO 1370 
[eGSs="y" OR GS="YES” THEN [970 
MOVE X2¢(C2).Y1¢C2) 
DRAW X2¢(C2) , Y2¢(C2) 
DRAW X1(C2),¥2¢(C2) 
CLIP ON 
AXES X4,Y4,X1(C2) , Y1(C2) 
Cie OEE 
VIEWPORT 13,.100#RATIO.S.90 
WINDOW X1(C2) ,X2(C2) ,¥1¢(C2) , Y2¢C2) 
LORG 4 
GoiZze 2, .5 
FOR X=Xt«(C2) TO X2¢C2) STEP X4 
CEP OFF 
MOVE X,Y1I¢C2) 
LABEL USING “K";X f PLOTS X-DIGITS 
NEXT X 
VIEWPORT 4,100*RATIO.10,90 
WINDOW X1(C2) ,X2¢C2) , Y1¢(C2) , Y2¢C2) 
LORG 2 
FOR Y=Y1«(C2) TO Y2¢C2) STEP Y4 
Cele OFF : 
MOVE X1¢(C2),Y 
PRBee USING “Kr 3 Y 
NEXT Y 
VIEWPORT 0,100"RATIO.O. 100 
WINDOW X16°2) ,X2¢C2) ,YI(C2) , Y2¢C2) 
Goze 3 
LORG 6 Seuss PETLE 
MOVE X2(C2)/2,Y2¢C2) 
LABEL TS 
DEG 
LDIR 90 
CS T7ee4.. 7 
MOVE X1(C2),¥2°C2)/2 
LABEL YS ' PLOTS Y-AXIS LABEL 
LORG 4 
LDIR 0 
MOVE X2(C2)/2,Y1¢C2) 
LABEL XS ' PLOTS X-AXIS LABEL 


PRINT “” 
PRINT "™ 
PRINT """ 
PRINT "* 
PRINT IS THE PLOT (AXES AND LABELS) UK?" 
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PRINT °° NOTE: IF YOU ANSKER YES OR Yr” 

PRINT °° PLOT WILL BE SENT TO” 

PRINT °* EXTERNAL PLOTTER” 

INPUT FS 

IF FS="N" OR FS="NO" THEN GOSUB Setup_plotter ! REDO PLOTTER SETUP 
i] 

GCLEAR * CLEARS PLOT FROM CRI 

i] 

GOSUB Ext_plotter 

t 

PENUP 

MOVE 0,0 

’ 

GOSUB Plot_data 

RETURN 

? 

) BSDPDPPPSDPHPDPHHHPPSHADAHDHHAHHPHPHeHPPAPsHPHBsAHRsHB Hse sw 

’ 
Ext_plotter: t DUMPS PLOTTER SETUP ON EXTERNAL HP-7225B PLOTTER 


Bide teR 91S) 7/05, rok” 
VIERKPORT 13.100*RATIO.10.930 
WINDOW X1¢C2) ,X2¢(C2).Y1(C2) , Y2¢C2) 
IF GS="N" OR GS="NO" THEN 2739 
FOR [=X4/B0+X1<(C2) TO X2¢C2) STEP X4/B0 
MOVE 1, ¥1¢Cz) 
_ DRAW I,Y2¢(C2) 
NEXT [ 
FOR I[=¥4/80+Y¥1(C2) TO Y2¢C2) STEP Y4/B0 
HOVE X1(C2),I 
DRAW X2¢C2),I 
NEXT I 
GOTO 2760 
IF GS="Y" OR GS*"YES* THEN 2760 
MOVE X2(C2),¥1(C2) 
DRAW X2(C2) , Y2¢C2) 
DRAW X1(C2), ¥2¢(C2) 
CLIP ON 
AXES X4.Y¥4,X1(C2) , ¥1(C2) 
Cite OFF 
VIEWPORT 13,100#RATIO,S,90 
WINDOW X1¢(C2) ,X2(C2),Y1(C2) , Y2¢C2) 
LORG 4 
CSeze 3..6 
FOR X=X1(C2) TO X2(C2) STEP X4 
Sie OFF 
MOVE X,Y1¢(C2) 
LABEL USING “K"";X 
NEXT X 
VIEWPORT 4,100#RATIO, 10,90 
WINDOW X1(C2) ,X2(C2) ,Y1¢€C2) , Y2(C2) 
LORG 2 
FOR Y=¥1(C2) TO Y2¢«C2) STEP Y4 
CietP OFF 
MOVE X1(C2),Y 
LABEL USING “K";Y 
NEXT Y 
VIEWPORT 0.100#RATIO.O.100 
WINDOW X1¢C2) ,X2¢C2) , Y1(C2) , Y2¢C2) 
LORG 2 


og 





Zeke 


3040 
3050 
3060 
3070 
3080 
3090 
3100 
3110 
3120 
3130 
3140 
3150 
3160 
3170 
3180 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
3260 
3270 
3280 
290 
3300 
3310 
3320 
3330 
3340 
3350 
3360 
3370 
3380 
3390 
3400 
3410 
3420 
3430 
3440 
3450 
3460 
3470 
3480 
3490 


SEs 
LORG 6 
Maeve A2ZCE2)72,¥2¢C2) 
LABEL IS 
DEG 
LDIR 90 
CSIZE 4. =/ 
Mave Xi (Ge) .'2(C2)72 
LABEL YS 
LORG 4 
LDIR 0 
MOVE X2(C2)/2,Y1(C2) 
, LABEL XS 


PENUP 
MOVE 0,0 


RETURN 
? 


? BDU DDDBSBDMBBUMHSBSHSARSHSHSHSHDRBPBHHKHSHBHSHSHHSHS De 


i] 


Plot_data: t PLOTS TO DATA SELECTED 
? 


YVIEHPORT 13.100*RATIO, 10.90 
WINDOW X1(C2) ,X2(C2) ,.Y1(C2) , Y2(C2) 
i} 


IF C2=! THEN 3290 
IF C2=2 THEN 3340 
IF C2=3 THEN 3390 


FOR I=t TO K 

PLOT Cod(I) ,Lod(TI) 
NEXT I 
GOTO 3430 


FOR I=! TO K 
PLOT Def1l¢I) ,Lod(I) 
NEXT I 
GOTO 3430 
t 


FOR I=1 To C1 


PLOT Delta_a(K2) ,J(K2) 
NEXT [I 
! 


RETURN ’ RETURNS TQ SETUP_PLOTTER 


’ 3D WD WD WD WD WD TD WD WD WD 2D UD 4D UE 2D HE OO 10 WD WD TD WD OD 1 2D WD 40 0 I OD 4 2D 2D 08 4 4 2 0 0 


’ 


END 
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APPENDIX G 


PRECRACK OUTPUT 


USING ENGLISH UNITS (in.lbs,psi) TRE FOLLOWING INFCRMAT ION 
HAS BEEN ENTERED FOR: HY80-5B 


WIDTH = .9469 NOTCH LENGTH = .444 

Bmax = .496 Bmin = .394 

HALF SPAN = 1.894 PQISSON’S RATIO = .3 - 
YIELD STRENGTH = 86600 ULTIMATE TS = 100000 
ELASTIC MODULUS = 2.9E&+7 FLOW STRESS = 93300 


CONVERTING THIS DATA TO METRIC UNITS (mm.N.KPa) FOR: HY80-5B 
WIDTH = 24.05126 NOTCH LENGTH = 11.2776 

Bmax = 12.5984 Bmin = 10.0076 

YIELD STRENGTH = 597085.603738 ULTIMATE TS = 689475.293 
ELASTIC MODULUS = 1.9994783497E+8 FLOW STRESS = 643280.448369 


BASED ON CURRENT INFORMATION FOR HY80-SB: 


CURRENT CRACK LENGTH = .444 in 

GPTIMUN CRACK LENGTH = .563814 in 

EXTENSION REQUIRED = .12414 in 

CRACK LENGTH AT 50% OF EXTENSION = .50607 in 

Pmax = 3272.42651647 Lbs Pmin = 163.621325823 Lbs 
Delo = 3108.80519064 Lbs Kmax = 32563.7466931 
DELTA_K = 78435.5593585 in”.S DELTA_KE = 40000 in”.5 


NEW LOADS BASED ON DELTA_KE ARE: 
Pmax = 1668.84843723 Lbs Pmin = 83.4424218615 Lbs 
DELTA_P = 1585.40601537 Lbs Kmax = 42105.2631579 


CONVERSION DATA 
FOR LOAD:-1779.288646 N/V (-400 Lbs/V) 


' FOR COD: .0508 mm/V ( .002 in/V) 


FOR DEFLECTION:-12.7 nm/V (-.5 in/V) 


U = .07 
a/W = ./732186117982 in CRACK LENGTH = .693307035117 in 
LOAD (lbs) * COD (in) COD «Volts? 
200 .00256290169302 1.28145084651 
300 .0038 4435253953 1.92217626976 
400 .00512580338604 2.562901693%2 
500 .00640 725423255 3.20362711627 
600 .00768870507906 3.84435253353 
700 .00897015592557 4.48507796278 
800 .0102516067721 5.12580338604 
900 -0115330576186 S.76652880929 
1900 .0128145084651 6 .40725423255 
1100 .0140959593116 7.0479796558 
1200 .0153774101581 7.68870507906 
1300 .0166588610046 8 .32943050231 
1400 .0179403118511 8.97015592557 
1500 .0192217626976 9.610881348382 
1600 .0205032135441 10.2516067721 
1700 .0217846643907 10.8923321953 
1 ae .0230661152372 11.5330576186 
a/W = .713968763392 in CRACK LENGTH = .676057022056 in 
LOAD (lbs) COD ‘ind COD (Volts) 
200 .00220863052208 1.10431526104 
300 .00331294578312 1.65647289156 
4N0 .00441726104416 2.20863052208 
500 .00552157630521 2.7607881526 
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Sim ENGLISH UNITS ¢in.tbs.esi) THE FOLLOWING INFORMATION 
HAS BEEN ENTERED FOR: HY80-5SB 


WIDTH = .9469 NOTCH LENGTH = .444 

Bmax = .496 Bmin = .394 

HALF SPAN = 1.394 PQISSON’S RATIO = .3 

YIELD STRENGTH = 86600 ULTIMATE TS = 100000 
ELASTIC MODULUS = 2.93E+7 FLOW STRESS = 93300 


CONVERTING THIS DATA TO METRIC UNITS (mm.N.KPad FOR: HY80-5B 
WIDTH = 24.05126 NOTCH LENGTH = 11.2776 

Bmax * 12.5984 Bmin = 19.0075 

YIELD STRENGTH = S597085.603738 ULTIMATE TS = 689475.293 
ELASTIC MODULUS = 1.9994783497E+8 FLOK STRESS = 643280.448263 


BAScD ON CURRENT INFORMATION FOR HYS80-SB: 


CURRENT CRACK LENGTH = .444 in 

GPTIMUN CRACK LENGTH = .56814 in 

EXTENSION REQUIRED = 712414 in 

CRACK LENGTH AT 50% OF EXTENSION = .50607 in 

Pmax = 3272.42651647 Lbs Pmin = 163.621325823 Lbs 
Delp = 3108.80519064 Lbs Kmax = 82563.7466931 
DELTA_K = 78435.5593585 in’.S DELTA_KE = 40000 in*.5 


NEW LOADS BASED ON DELTA_KE ARE: 
Pmax = 1668.84843723 Lbs Pmin = 83.4424218615 Lbs 
DELTA_P = 1585.40601537 Lbs Kmax = 42105.2631S79 


CONVERSION DATA 
FOR LOAD:-1779.288646 N/V (-400 Lbs/V) 
FOR COD: .0508 mm/V ( .002 in/V¥) 
FOR DEFLECTION: -12.7 mm/V (-.5 in/V) 


TRE FOLLOWING INFORMATION HAS BEEN OBTAINE?: 

CRACK LENGTH = .477215195165 in FROM MEASUREMENT = 1 

WITH SSZCONFIDENCE INTERVAL FROM .47792545433 TO .476507410981 in 
USING COD COMPLIANCE OF: 2.841943110S54E-6 in/Lb = 1.62905173472E-S mm/N 
WITH A CORRELATION OF: .999771386472 FROM 106 DATA POINTS. 


THe FOLLOWING INFORMATION HAS BEEN OBTAINED: 

CRACK LENGTH = .477198279382 in FROM MEASUREMENT # 2 

WITH SSZCONFIDENCE INTERVAL FROM .478064054454 70 .476336178981 in 
USING COD COMPLIANCE OF: 2.84076126094E-6 in/Lb = 1.63037750213E-S mm/N 
WITH A CORRELATION OF: .99959257929 FROM 124 DATA POINTS. 


BASED ON CURRENT INFORMATION FOR HY80-SB: 


CURRENT CRACK LENGTH = .477198279382 in 

UPTIMUN CRACK LENGTH = .56314 in 

EXTENSION REQUIRED = .0909417206182 in 

CRACK LENGTH AT SO0Z OF EXTENSION = .489470860309 in 
Pmax = 2854.63728597 Lbs Pmin = 142.731864298 Lbs 
Deip = 2711.90542167 Lbs Kmax = 80310.2753096 
DELTALK = 76294.7615441 in°.S DELTA_KE = 40000 in’.S 


NEW LOADS BASED ON DELTA_KE ARE: 


Pmax = 1496.65606161 Lbs Pmin = 74.8319030804 Lbs 
DELTA P = 1421 .30425353 Lbs Kmax = 42195.2631579 
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Teo TLL SRG INFORMA LON cAS SEEN JBTATNES: 

CRACK LENGTH = ,499315101298 in FROM MEASUREMENT = 1 

WiTH SSZCONFLDENCE INTERVAL FROM .500434044632 TO .498202259284 in 
USING GOD COMPLIANCE OF: 3.23884557474E-6 in/lLb = 1.386189031428E-S mm/N 
WITH A CORRELATION OF: .399477933275 FROM 24 OATA POINTS. ; 


THE FOLLOWING INFORMATION HAS BEEN OBTATHNED: 

CRACK LENGTH = .499773214641 in FROM MEASUREMENT # 2 

WITH SSZCONFIDENCE INTERVAL FROM .50177796077 TO .497787999797 in : 
USING COD COMPLIANCE OF: 3.2477593972936E-6 in/Lb = 1.8769892409E-S mm/N 
WITH A CORRELATION OF: .999919979378 FROM S57 ‘DATA POINTS. 


BASED ON CURRENT INFORMATION FOR HY80-SB: 


CURRENT CRACK LENGTH = .499773214641 in 

OPTIMUN CRACK LENGTH = .56314 in 

EXTENSION REQUIRED = .0633657853593 in 

CRACK LENGTH AT 50% OF EXTENSION = .47818339268 in 
Pmax = 1034.73227996 Lbs Pmin = S1.7236512998 Lbs 
Deip = 982.995665963 Lbs Kmax * 31462.7523142 
DELTA_LK = 29889.6145985 in°.5 DELTA_LKE = 29000 in”™.S 


NEW LOADS BASED ON DELTA_KE ARE: 
Pmax = 1003.9351936 Lbs Pmin = 50.1967596/798 Lbs 
DELTA_P = 952.738433917 Lbs Kmax * 90526.3157895 


TRE FOLLOWING INFORMATION HAS BEEN OBTAINED: 

CRACK LENGTH = .52323832S5069 im FROM MEASUREMENT #& 1 

WITH SSZCONFIDENCE INTERVAL FROM .524952826657 TO .521538109085 in 
USING COD COMPLIANCE OF: 3.74604099481E-6 in/Lb = 2.16181567648E-S mm/N 
WITH A CORRELATION OF: .9991834S56903 FROM 6t DATA POINTS. 


THE FOLLOWING INFORMATION HAS BEEN OBTAINED: 
CRACK LENGTH = .519260644793 in FROM MEASUREMENT # 2 
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WITH Q5%CONFIDENCE INTERVAL FROM .521257634776 TO .51723301U25 in 
USING COD COMPLIANCE OF: 3.65545467572E-5 in/Lb = 2.11309694446—E-5 mm/N 
WITH A CORRELATION OF: .998905179767 FROM 61 DATA POINTS. 


THE FOLLOWING INFORMATION HAS BEEN OBTAINED: 

CRACK LENGTH = .520164834896 in FROM MEASUREMENT * 3 

WITH 3SACONFIDENCE INTERVAL FROM .S52173216005f TO .51860945694 in 
USING COD COMPLIANCE OF: 3.67581279234E-6 in/Lb = 2.11927866293E-5 mm/N 
WITH A CORRELATION OF: .999107091656 FROM 73 DATA POINTS. 


BASED ON CURRENT INFORMATION FOR HY80-58: 


CURRENT CRACK LENGTH = .520164834896 in 

OPTIMUN CRACK LENGTH = .56814 in 

EXTENSION REQUIRED = .0479751651041 in 

CRACK LENGTH AT 502% OF EXTENSION = .467987582552 in 
Pmax = 942.504620395 Lbs Pmin = 47.1252310197 Lbs 
Delp = 895.379389375 Lbs Kmax = 30831.6649108 
DELTA_K = 29290.0816653 in®™.5 DELTA_KE = 29000 in”.5 


NEW LOADS BASED ON DELTA_KE ARE: 
Pmax = 933.170289649 Lbs Pmin = 46.6585144825 Lbs 
DELTA_P = 886.51/1775167 Lbs Kmax = 30526.3157895 


THE FOLLOWING INFORMATION HAS BEEN OBTAINED: 

CRACK LENGTH = .557542180253 in FROM MEASUREMENT - 1 

WITH SSACONFIDENCE INTERVAL FROM .559342565208 TO .555757596637 in 
USING COD COMPLIANCE OF: 4.6522934609E-6 in/Lb = 2.6877228564E-5 mm/N 
WITH A CORRELATION OF: .999026349545 FROM 60 DATA POINTS. 


THE FOLLOWING INFORMATION HAS BEEN OBTAINED: 
CRACK LENGTH = .556985314799 in FROM MEASUREMENT # 2 
WITH SS5ZCONFIDENCE INTERVAL FROM .559603140017 TO .554400760935 in 
«USING COD COMPLIANCE OF: 4.63555555373E-6 in/Lb = 2.69227537518E-5 mm/N 
WITH A CORRELATION OF: .998286279816 FROM 51 DATA POINTS. 


THE FOLLOWING INFORMATION HAS BEEN OBTAINED: 

CRACK LENGTH = .55746232304 in FROM MEASUREMENT # 3 

WITH SSZCONFIDENCE INTERVAL FROM .558944733365 TO .555990641885 in 
USING COD COMPLIANCE OF: 4.64988856512E-6 in/Lb = 2.68079074153E-5 mm/N 
WITH A CORRELATION OF: .999219985095 FROM 70 DATA POINTS. 


BASED ON CURRENT INFORMATION FOR HY80-5B: 


CURRENT CRACK LENGTH = .55746232304 in 

OPTIMUN CRACK LENGTH = .56814 in 

EXTENSION REQUIRED = .0106776769602 in 

CRACK LENGTH AT 50% OF EXTENSION = .449338833848 in 
Pmax = 784.95101444 Lbs Pmin = 39.247550722 Lbs 

Delp = 745.703463718 Lbs Kmax = 29594.1425953 
DELTA_K = 28114.4354655 in*.5 DELTA_LKE = 29000 in*.5 


THE FOLLOWING INFORMATION HAS BEEN OBTAINED: 

CRACK LENGTH = .556873246958 in FROM MEASUREMENT + 1 

WITH 354CONFIDENCE INTERVAL FROM .558630143876 TO .555131398279 in 
USING COD COMPLIANCE OF: 4.63219612139E-6 in/Lb = 2.67532481877E-5 mm/N 


STi 





W(TH A CORRELATION JF: .999174171419 FROM 54 DATA POINTS. 
BASED SN CURRENT [INFORMATION FOR HYSO-SB: 


CURRENT CRACK LENGTH = .$56873246958 in 

GPTIMUN CRACK LENGTH = .563!4 in 

EXTENSION REGUIRED = 0112667530417 Ln 

CRACK LCENGTH AT 50% OF EXTENSIGN » .449623376521 in 
Pmax = 787.327495262 Lbs Pmin = 39.365374763i Lbs 
Delp = 747.961320499 Lbs Kmax 2 29614.5202835 
DELTA_K = 28133.39037694 in°.S DELTA_KE = 29000 1in’.5 


TrE FOLLOWING INFORMATICN HAS BEEN OBTAINED: 

CRACK LENGTH = .560000734025 itn FROM MEASUREMENT = 1 

WiTH 95 .CONF iDENCE INTERVAL FROM .56106384425 [0 .55893721925 in 

USING COD COMPLIANCE OF: 4.72709402318E-5 in/Lb = 2.71897983847E-S mm/i 
WITH A CORRELATION OF: .999396466938 SROM 101 DATA POINTS. 


TRE FOLLOWING INFOQPMATION HAS BEEN OBTAINED: 

CRACK LENGTH = .562208107534 in FROM MEASUREMENT # 2 

WITH SSZCONFIDENCE INTERVAL FROM .56543540602 10 .5609881270658 in 
USING COD COMPLIANCE OF: 4.795529965S5S1E-5 in/Lb = 2.76035025025E-5 mm/h 
WITH 9 CORRELATION OF: .99948124671 FROM 57 DATA POINTS. 


SASED ON CURRENT INFORMATION FUR HY3O-5e: 


CURRENT CRACK LENGTH = .562208107534 in 

GP TIMUN CRACK LENGTH = .56814 in 

EXTENSION REQUIRED = .00593199246594 in 

CRACK LENGTH AT 50% OF EXTENSION = .446965946233 in 
Pmax = 765.936366326 Lbs Pmin + 38.2968183162 Lbs 
Delp = 727.63954801 Lbs Kmax = 29428.399829 

DELTA_K = 27957.4548376 in’.S DELTA_KE = 29000 in°.5 


TRE FOLLOWING INFORMATION HAS BEEN UBTAINED: 

CRACK LENGTH = .5641513471472 in FROM MEASUREMENT « 1 

WITH SSZCONFIDENCE INTERVAL FROM .565674835605 TQ .562639453926 in 
USING COD COMPLIANCE OF: 4.85680462792E-6 in/Lb = 2.801124373Z9E-5 mm/N 
WITH A CORRELATION OF: .999133273453 FROM 72 DATA POINTS. 


Ca] 
rd 


THE FOLLGWING INFORMATION HAS BEEN OBTAINED: 

CRACK LENGTH = .568361917225 in FROM MEASUREMENT ¢ 2 

WITH SSZCONFIDENCE INTERVAL FROM .569305393318 TO .567422811142 in 
USING COD COMPLIANCE OF: 4.99294827909E-6 in/Lb = 2.86883916896E-5 mm/N 
WITH A CORRELATION GF: .99960869725S4 FROM 84 DATA POINTS. 


PAECReCMINGe bs COMPLETED! 


FINAL CRACK LENGTH = .568361917225 
CURRELATION OF DATA = .999601697254 


REMOVE SPECIMEN FROM MIS. 
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APPENDIX H 


J_INTGL OUTPUT 


THREE POINT BEND SPECIMEN 


TEST REPORT 


SPECINEN 1.D. HY8O- 5&8 _ matertar_HYY 9O ORIENTATION 7-S 


Date_/2% 113 8.3 Operator WE Mas CE Teat Teap_2/ C¢ 


= = 


(Load pt.) 





w= 0.9469 nx 0.444 Fe >, ee ee” 


\ 


Ha * 1.g94* Baax ® 0-496G0- Bain # 6.394" 





Ys = 8b66CO pe Uta = /00 COO pie Fa = 937300 fle 


Poiason’a ratio # C. 3 Elaatic Modulua = 2-9 K/0" ae 
Conamentsa: 
wm... a Pre 
Rae “ : J- Tesr . / } > e142 KY,,% 








Se 





WOING ENGLISH UNITS Cin. lbos.asi) TSE FULLUWING INFORMATION 
HAS BEEN ENTERED FOR: HY50-5B 


WIDTH = .9463 NOTCH CENGTH = .444 

Bmax = .496 Bmin = .394 

HALF SPAN = 1.894 POISSON'S RATIO = .3 

YIELD STRENGTH = 86606 ULTIMATE TS = 100000 
ELASTIC MODULUS = 2.39E+7 FLOW STRESS = 93300 


CONVERTING THIS DATA TO METRIC UNITS (mm.N,KPa) FOR: HY80-5B 


WIDTH = 24.05126 MOTCH LENGTH = 11. 2776 
Bmax = 12.5984 Bmin = 10.0076 
HALF SPAN = 48.1076 POISSON’S RATIGO = .3 


YIELD STRENGTH = 597085.603738 ULTIMATE TS * 689475.293 
ELASTIC MODULUS = 1.9994783497E+8 FLOW STRESS = 643280.448369 


. CONVERSION DATA 
FOR LOAD:-1779.288646 N/V ( 400 Lbs/V) 
FOR COD: .254 mm/V ( .Of in/Y) 
FOR DEFLECTION:-1.27 mm/V ( .05 in/V) 


GIVEN a= .568 inches THEN: - 
Pmax= 895.290182804 Lbs FROM ASTM £813-81 


THE FOLLOWING INFORMATION HAS BEEN OBTAINED: 

INITIAL CRACK LENGTH = 14.3875382557 mm FROM MEASUREMENT = 1 

WITH 95% CONFIDENCE INTERVAL FROM 14.4079726404 10 14.3671845862 mm 
USING COD COMPLIANCe OF: 2.81520064535E-5 mm/N 

WITH A CORRELATION OF: .999323843963 FROM 189 DATA POINTS. 


THE FOLLOWING INFORMATION HAS BEEN OBTAINED: 
INITIAL CRACK LENGTH = 14.4326755548 mm FROM MEASUREMENT # 2 


WITH 35% CONFIDENCE INTERVAL FROM 14.4442813846 TO 14.42109582 mm 


USING COD COMPLIANCE OF: 2.84829954475E-5 mm 
WITH A CORRELATION OF: .99984474937 FROM 134 DATA POINTS. 
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INITIAL CRACK LENGTH = 14.4455712237 mm FROM MEASUREMENT - 3 

WITH 95% CONFIDENCE INTERVAL FROM 14.456550679!1 TO 14.434615129 mm 
USING COD COMPLIANCE OF: 2.85784514/764E-5 mm/N 

WITH A CORRELATION OF: .99981652255/7 FROM 176 DATA POINTS. 


---- UNLOADING #1 ---- 


CRACK LENGTH = 14.4487507815 mm DELTA _a = .016075226715 mm 
95% CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 14.4693514291 TO 14.4282322251 


TOTAL # DATA POINTS = 632 # DATA POINTS THIS UNLOADING = 128 
CORRELATION = .999529434569 
COD COMPLIANCE = 2.87555293817E-5 mm/N 


AREA = 711.365042717 N*mm 

J = 14.8049806501 KJ/m"2 

MAX LOAD LINE DEFLECTION = .24034834255 mm 
CHANGE IN COD = .00109982 mm 


---- UNLOADING # 2 ---- 


CRACK LENGTH = 14.4711513219 mm DELTA _a = .0384757670308 mm 
95% CONFIDENCE INTERV@L ON CRACK LENGTH IS FROM 14.4927078135 TO 14.4496847244 


TOTAL # DATA POINTS = 947 # DATA POINTS THIS UNLUADING = 132 
CORRELATION = .999467111958 
COD COMPLIANCE = 2.89307913363E-5 mm/N 


AREA = 1693.69098454 N*mm 

dee 35.2970225265 KJ/m 2 

MAX LOAD LINE DEFLECTION = .3823475560072 mm 
CHANGE IN COD = .0008382 mm 


---- UNLCADING #3. ---- 


CRACK LENGTH = 14.4875164202 mm DELTA _a = .0548408654168 mm 
95% CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 14.5143346851 TO 14.460837176 


TOTAL # DATA POINTS = 1300 * DATA POINTS THIS UNLOADING = 157 
CORRELATION = .999017363069 
COD COMPLIANCE = 2.90942709989E-5 mm/N 


AREA = 2033.0003923 N#*mm 

J = 63.28381022 KJ/m°2 

MAX LOAD LINE DEFLECTION = .543433314464 mm 
CHANGE IN COD = .0009398 mm 


---- UNLOADING # 4 ---- 


CRACK LENGTH = 14.5414951598 mm DELTA _a = .108819604992 mm 
95% CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 14.5601316376 TO 14.5229259715 


TOTAL # DATA POINTS = 1670 * DATA POINTS THIS UNLOADING = 109 
CORRELATION = .999675611283 
COD COMPLIANCE = 2.94442988118E-5 mm/N 


AREA = 5208.32927275 Nemm 


wee 108.998439815 KJ/m' 2 
Poel no CINE VPSEECTION = .7/S231829029= mm 


MO 


CHANGE IN COD = .0007874 mm 


---- UNLOADING #5 ---- 


CRACK LENGTH = 14.7407494376 mm DELTA _a = .308073862773 mm 
95% CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 14.7587196422 [0 14.7228419911 


TOTAL * DATA POINTS = 2060 * DATA POINTS THIS UNLOADING = 95 
CORRELATION = .999731105207 
COD COMPLIANCE = 3.10248861726E-5 mm/N 


AREA = 8216.7400787 N#mm 

J = 169.858689852 KJ/m~2 

MAX LOAD LINE DEFLECTION = 1.10455056629 mm 
CHANGE IN COD = .0008636 mm : 


---- UNLOADING + 6 ---- 


CRACK LENGTH = 14.8152586049 mm DELTA _a = .382583050059 mm 
95% CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 14.8481840789 10 ,14.782543466 


TOTAL # DATA POINTS = 2280 # DATA POINTS THIS UNLOADING = 138 
CORRELATION = .998642018089 
COD COMPLIANCE = 3.1771921946i1E-S mm/N 


AREA = 9001.97373669 N#mm 

J = 185.342172083 KJ/m°2 

MAX LOAD LINE DEFLECTION = 1.18989063161 mm 
CHANGE IN COD = .0009906 mm 


soa  CNLORDTNG 4 7 ~----- 


CRACK LENGTH = 14.8695206824 mm DELTA _a = .436845127588 mm 
95% CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 14.8849566722 TO 14.8541311238 


TOTAL # DATA POINTS = 2531 # DATA POINTS THIS UNLOADING = 147 
CORRELATION = .999677856393 
COD COMPLIANCE = 3.20855798284E-S mm/N 


AREA = 9678.11102926 Nemm 

J = 198.876552682 KJ/m"2 

MAX LOAD LINE DEFLECTION = 1.26297067275 mm 
CHANGE IN COD = .0004572 mm 


ooo | UMLBADING * 3 ---- 


CRACK LENGTH = 15.0148445285 mm DELTA _a = .582168973634 mm 
95% CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 15.0416778385 TO 14.9881516089 


TOTAL #« DATA POINTS = 2761 * DATA POINTS THIS UNLOADING = 97 
CORRELATION = .999364204197 
COD CUMPLIANCE = 3.34673392487E-5 mm/N 


AREA = {0586.6765231 N#mm 


J = 215.448776655 KJ/m° 2 
MAX LCAD LINE DEFLECTION = 1.36396114491 mn 
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CHANGE IN COD = .0006604 mm 


---- UNLOADING #9 ---- 


CRACK LENGTH = 15.1216765395 mm DELTA _a = .689000984704 mm 
95% CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 15.1358606171 TO 15.1075318754 


TOTAL # DATA POINTS = 3042 * DATA POINTS THIS UNLOADING = 119 
CORRELATION = .99977636544 
COD COMPLIANCE = 3.43343078736E-S5 mm/N 


AREA = 11709.4655627 Nemm 

J = 237.699077286 KJ/m°2 

MAX LOAD LINE DEFLECTION = 1.49020183752 mm 
CHANGE IN COD = .0005842 mm — 


---- UNLOADING # 10 ---- 


CRACK LENGTH = 15.2624877891 mm DELTA _a = .829812234314 mm 
95% CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 15.2760292626 10 15.2489823604 


TOTAL # DATA POINTS = 3394 # DATA POINTS THIS UNLOADING = 167 
CORRELATION = .999700915055 
COD COMPLIANCE = 3.56785340653E-5 mm/N 


AREA = 12983.4083914 N#mm 

J = 262.394826944 KJ/m"2 

MAX LOAD LINE DEFLECTION = 1.63760345064 mm 
CHANGE IN COD = .0007366 mm . 


eas UNEGRUUNG ¥ 11. ---- 


CRACK LENGTH = 15.3992548961 mm DELTA _a = .966679341294 mm 
95% CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 15.4107739959 TO 15.3879615278 


TOTAL # DATA POINTS = 3867 * DATA POINTS THIS UNLOADING = 204 
CORRELATION = .999734089719 
COD COMPLIANCE = 3.703512357S5S9E-5 mm/N 


AREA = 14557.9277896 N=mm 

J = 294.038000741 KJ/m*2 

MAX LOAD LINE DEFLECTION = 1.82262379631 mm 
CHANGE IN COD = .0004318 mm 


---- UNLSADING # 12 ---- 


CRACK LENGTH = 15.669999544 mm DELTA _a = 1.2373239892 mm 
35% CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 15.6847004127 TQ 15.6553416287 


TOTAL # DATA POINTS = 4267 # DATA POINTS THIS UNLOADING = 154 
CORRELATION = .999652892868 
COD COMPLIANCE = 4.00050024785E-5 mm/é 


AREA = 16330.4160222 N=mm 


see 32D. 842585525 Ki/m 2 
MAX LOAD LINE DEFLECTION = 2.03825676401 mm 
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CHANGE IN COD = .0013462 mm 


---- UNLOADING # 13 ---- 


CRACK LENGTH = 15.8781418346 mm DELTA _a = 1.44546627973 mo 
95% CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 15. 3901071812. TQ 15.8662051558 


TOTAL # DATA POINTS = 4710 * DATA POINTS THIS UNLOADING = 208 
CORRELATION = .999677232181 
COD COMPLIANCE = 4.24377123766E-S5 mm/N 


AREA = 17933.7183158 Nemm 

J = 355.360006494 KJ/m*2 

MAX LOAD LINE DEFLECTION = 2.24213294335 mm 
CHANGE IN COD = .000508 mm ° 


---- UNLOADING # 14 ---- 


CRACK LENGTH = 16.1357477216 mm DELTA _a = 1.7030721668 mm 
95% CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 16.1463233402 T0 16.1251947229 


TOTAL # DATA POINTS = 5206 # DATA POINTS THIS UNLOADING = 239 
CORRELATION = .999696024845 
COD COMPLIANCE = 4.57508725867E-5 mm/N 


AREA = 19646.0705407 Nemm 

J = 385.620902249 KJ/m°2 

MAX LOAD LINE DEFLECTION = 2.47449867471 mm 
CHANGE IN COD = .0004826 mm 


---- UNLOADING # 15 ---- 


CRACK LENGTH = 16.4193583405 mm DELTA _a = 1.98668272568 mm 
95% CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 16.4364102859 TO 16.4022658462 


TOTAL # DATA POINTS = 5648 # DATA POINTS THIS UNLOADING = 194 
CORRELATION = .999324290353 
COD COMPLIANCE = 4.99261698493E-S mm/N 


AREA = 21264.009056 N*mm 

J = 412.0835098069 KJ/m°~2 

MAX LOAD LINE DEFLECTION = 2.70365472703 mm 
CHANGE IN COD = .0005334 mm 


---- UNLOADING # 16 ---- 


CRACK LENGTH = 16.9028626935 mm DELTA _a = 2.47018713865 mm 
95% CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 16.9171113874 TO 16.3888656585 


TOTAL # DATA POINTS = 6102 # DATA POINTS THIS UNLOADING = 162 
CMRRELATION = .999564028818 
COD COMPLIANCE = 5.8035G6094828E-5 mm/N 


AREA = 22950.1834392 N*mm 


J = 428.162914109 Ku/m"2 
MAX LUAD LiNE DESLECTION = 2.96034280C60S mm 
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DATA FROM J vs DELTA_a PLOT 


SLOPE = 164.148013927 VeTNTERGER? = 123. 935075018 
CORRELATION = .99667195294 FROM 9 DATA POINTS. 


Jie = 142.05774949 


TEST IS OVER! 
UNLOADING # DEL TA_a(mm) J(KS/m~2) 
.016075 14.8 


1 
2 038476 35.0 
3 054841 63.3 
4 . 108820 109.9 
=) . 308074 169.9 
6 "eG 2503 185.3 
7 436845 198°9 
8 982169 215.4 
3 -689001 23/.7 
10 829812 <62.4 
11 966679 294.0 
12 1.237324 325.4 
13 1.445466 355.4 
14 1.703072 385.6 
Je) 1.986683 412.1 
16 2.470187 428.2 


DATA FROM J vs DELTA_a PLOT 


SLOPE = 164.148013927 TOINTERCEPT 28 8122 958075018 
CORRELATION = .99667195294 FROM 9 DATA POINTS. 
Jic = 142.05774949 


feo! 1S OVER! 
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APPENDIX I 


J_DATA OUTPUT 


WHAT DATA FILE DO YOU WISH TO READ? 
J=INTEGRAL FEST OF AYSO-SB BY: KURT 
DANE tele -lo-os Tine 21035 


a0= 14.4326755548 WITH 95% CONDIFENCE INTERVAL 
FROM 14.4442813846 TO 14.42109582 


FROM COMPLIANCE = 2.84829954475E-5 
CORRELATION= .99984474937 FROM 176 DATA POINTS 
DOD YOU WANT A HARD COPY OF THIS INFORMATION? 


FROM UNLOADING # 1! RECORDED 128 DATA POINTS 
COD COMPLIANCE = 2.87555293817E-S mm/N CORRELATION = .999529434569 
CRACK LENGTH = 14.4487507815 mm DELTA_A = .016075226716 mm 
95% CONFIDENCE INTERVAL FROM 14.4693514291 mm TO 14.4282322251 mm 
AREA = 711.365042717 N«mm J = 14.8049806501 KPa*m 

DATA POINT INDEXES FROM 504 TO 632 WERE USED. 


FROM UNLOADING + 2 RECORDED 132 DATA POINTS 

COD COMPLIANCE = 2.89307913363E-5 mm/N CORRELATION = .999467111958 
CRACK LENGTH = 14.4711513219 mm DELTA_A = .0384757670308 mm 

95% CONFIDENCE INTERVAL FROM 14.4927078135 mm TO 14.4496847244 mm 
AREA = 1693.69098454 N*mm J = 35.2970225265 KPa*m 

DATA POINT INDEXES FROM 815 TO 947 WERE USED. 


FROM UNLOADING ¥ 3 RECORDED 157 DATA POINTS 

COD COMPLIANCE = 2.90942709989E-5 mm/N CORRELATION = .9939017363069 
CRACK LENGTH = 14.4875164202 mm DELTA_A = .0548408654168 mm 

95% CONFIDENCE INTERVAL FROM 14.5143346851 mm T0.14.460837178 mm 
AREA = 3033.0003923 N«mm J = 63.28381022 KPa*m 

DATA POINT INDEXES FROM 1143 TO 1300 WERE USED. 


FROM UNLOADING # 4 RECORDED 109 DATA POINTS 

COD COMPLIANCE = 2.94442988118E-5 mm/N CORRELATION = .999675611283 
CRACK LENGTH = 14.5414951598 mm DELTA_A = .108819604992 mm 

95% CONFIDENCE INTERVAL FROM 14.5601316376 mm TO 14.5229259715 mm 
AREA = 5208.32927275 Nemm J = 108.998439815 KPa*m 

DATA POINT INDEXES FROM 1561 TO 1670 WERE USED. 


FROM UNLOADING # S RECORDED 95 DATA POINTS 
COD COMPLIANCE = 3.10248861726E-S mm/N CORRELATION = .999731105207 
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CRACK LENGTH = 14.7407494376 mm DELTA_A = .308973882773 mm 

95% CONFIDENCE INTERVAL FROM 14.7587196422 mm TO 14.7228419911 mm 
AREA = 3216.7400787 Nemm J = 169.858689852 KPa*m 

DATA POINT INDEXES FROM 1965 TU 2060 WERE USED. 


FROM UNLOADING # 6 RECORDED 138 DATA POINTS 

COD COMPLIANCE = 3.17719219461E-5 mm/N CORRELATION = .998642018089 
CRACK LENGTH = 14.8152586049 mm DELTA_A = .382583050059 mm 

95% CONFIDENCE INTERVAL FROM 14.8481840789 mm T0 14.782543466 mm 
AREA = 9001.97373669 Nemm J = 185.342172083 KPa*m 

DATA POINT INDEXES FROM 2142 TO 2280 WERE USED. 


FROM UNLGADING # 7 RECORDED 147 DATA POINTS 

COD COMPLIANCE = 3.20855798284E-5 mm/N CORRELATION = .999677856393 
CRACK LENGTH = 14.8695206824 mm DELTA_A = .436845127588 mm 

95% CONFIDENCE INTERVAL FROM 14.8849566722 mm 10 14.8541311238 mm 
AREA = 9678.11102926 Nemm J = 198.876552682 KPa*m 

DATA POINT INDEXES FROM 2384 TO 2531 WERE USED. 


FROM UNLOADING + 8 RECORDED 37 DATA POINTS 

COD COMPLIANCE = 3.346/°3392487E-5 mm/N CORRELATION = .999364204197 
CRACK LENGTH = 15.0148445285 mm DELTA_A = .582168973634 mm 

95% CONFIDENCE INTERVAL FROM 15.0416778385 mm TO 14.9881516089 mm 
AREA = 10586.6765231 Nemm J = 215.448776655 KPa*=m 

DATA POINT INDEXES FROM 2664 TQ 2761 WERE USED. 


FROM UNLOADING + 9 RECORDED 119 DATA POINTS 

COD COMPLIANCE = 3.43343078736E-S mm/N CORRELATION = .99977636544 
’' CRACK LENGTH = 15.1216765395 mm DELTA_A = .689000984704 mm 

95% CONFIDENCE INTERVAL FROM 15.1358606171 mm TO 15.1075318754 mm 
AREA = 11709.4655627 N*mm J = 237.699077286 KPa*#m 

DATA POINT INDEXES FROM 2923 TO 3042 WERE USED. 


FROM UNLOADING # 10 RECORDED 167 DATA POINTS 

COD COMPLIANCE = 3.56785340653E-5 mm/N CORRELATION = .999700915055 
CRACK LENGTH = 15.2624877891 mm DELTA_A = .829812234314 mm 

95% CONFIDENCE INTERVAL FROM 15.2760292626 mm [10 15.2489823604 mm 
AREA = 12983.4083914 Nemm J = 262.394826944 KPa*m 

DATA POINT INDEXES FROM 3227 TO 3394 WERE USED. 


FROM UNLOADING # 11 RECORDED 204 DATA POINTS 

COD COMPLIANCE = 3.70351235759E-5 mm/N CORRELATION = .999734089719 
CRACK LENGTH = 15.3993548961 mm DELTA_A = .966679341294 mm 

95% CONFIDENCE INTERVAL FROM 15.4107739959 mm TO 15.3879615278 mm 
AREA = 14557.9277896 Nxmm J = 294.038000741 KPa*m 

DATA POINT INDEXES FROM 3663 TO 3867 WERE USED. 


FROM UNLOADING # 12 RECORDED 154 DATA POINTS 

COD COMPLIANCE = 4.00050024785E-5 mm/N CORRELATION = .999652892868 
CRACK LENGTH = 15.669999544 mm DELTA_A = 1.2373239892 mm 

95% CONFIDENCE INTERVAL FROM 15.6847004127 mm TO 15.6553416287 mm 
AREA = 16330.4160222 Nemm J = 325.442585526 KPa*m 

DATA POINT INDEXES FROM 4113 TO 4267 WERE USED. 
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FROM UNLOADING + 13 RECORDED 208 DATA POINTS 

COD COMPLIANCE = 4.24377123766E-5 mm/N CORRELATION = .9996772221381 
CRACK LENGTH = 15.8781418346 mm DELTA_A = 1.44546627973 mm 

95% CONFIDENCE INTERVAL FROM 15.8901071812 mm TO 15.8662051558 mm 
AREA = 17933.7183158 Nemm J = 355.360006494 KPa*m 

DATA POINT INDEXES FROM 4502 TO 4710 WERE USED. 


FROM UNLOADING ¥ 14 RECORDED 239 DATA POINTS 

COD COMPLIANCE = 4.57508725867E-S mm/N CORRELATION = .999696024845 
CRACK LENGTH = 16.1357477216 mm  DELTA_A = 1.7030721668 mm 

95% CONFIDENCE INTERVAL FROM 16.1463233402 mm TO 16.1251947229 mm 
AREA = 19646.0705407 Nemm J = 385.620902249 KPa*m 

DATA POINT INDEXES FROM 4967 TQ 5206 WERE USED. 


FROM UNLOADING + 15 RECORDED 194 DATA POINTS 

COD COMPLIANCE = 4.99261698493E-5 mm/N CORRELATION = .999324290353 
CRACK LENGTH = 16.4193583405 mm DELTA_A = 1.98668278568 mm 

95% CONFIDENCE INTERVAL FROM 16.4364102869 mm TO 16.4023658462 mm 
AREA = 21264.009056 N*mm J = 412.083508069 KPa*m 


DATA POINT INDEXES FROM 5454 TO 5648 WERE USED. 


FROM UNLOADING # 16 RECORDED 162 DATA POINTS 

COD COMPLIANCE = 5.80350094828E-5 mm/N CORRELATION = .999564028818 
CRACK LENGTH = 16.9028626935 mm  DELTA_A = 2.47018713865 mm 

95% CONFIDENCE INTERVAL FROM 16.9171113874 mm TO 16.888656585 mm 
AREA = 22950.1834392 Nemm J = 428.162914109 KPa*m 

DATA POINT INDEXES FROM 5940 TO 6102 WERE USED. 
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DATA FROM av rOR:hY80-58 BY: XURT 
TESTED ON: 12-13-83 TIME: 10:35 
Pl # LOAD «NM» DEFLECTIUN (mm) COD (mm) 
1 2262.5 .079424 .060575 
2 226.7 .074706 .060576 
3 fafd So Wate .076740 .060637 
4 eras .0729003 .061003 
5 2350.8 .093469 .063589 
6 2409.0 .094973 .065014 
7 2442.3 nUSSS27 .065395 
8 2463.4 .09649} .066429 
S 2478.9 .096968 .066954 
10 2496.2 .097378 .0672830 
11 2505 sc .097656 .067627 
12 255.4 .098165 .067772 
13 2526.1 .093554 ~068207 
14 2533.6 .099114 »063511 
15 2554.0 .099864 .068963 
16 2570.4 .100177 .069411 
17 259 50S .101411 paves 
1g 2624.1 .102889 SOP NAO 7 
19 2653.4 OSS .072045 
20 2633.5 .104678 .072903 
21 LI i b .1096054 .073863 
ee 2746.9 .106721 .074519 
23 2764.1 »107052 n07 5028 
24 LIEV cs peed . 108077 .075443 
2S 2799.9 s106725 .075956 
2 Coc . 109596 .076655 
Pat 2837.3 .110482 nO77 Ge 
23 2854.0 .110667 N077617 
2b | 2873.4 SicUS5 078313 
30 2939.56 .114805 .080284 
31 2984.9 -115695 .081471 
6085 5686.4 2.9$4239 2.460422 
6086 5692.7 2.954786 2.460701 
6087 5639.1 22955204 aoe) ais) 
6088 5704.9 Zi sSlejs ley 45: 2.461616 
6089 5716.5 2.956736 2.462454 
6030 5731.8 22957730 2.463294 
6091 5747.5 2. 2a s0gs 2.464359 
6092 5761.2 2.959789 2.465172 
6093 5768.6 2.960445 2.465730 
6094 S777.0 2.961089 2.466442 
6095 S/e5.5 22261750 2.467051 
6096 S7o2ec 2.962262 2.467534 
6097 5803.0 2.963643 2.468321 
6098 5812.6 2.964142 2.469007 
6099 Shel) si als: 2.964692 2.469540 
6100 5824.1 2.965503 2eo70175 
6101 5830.0 2.965801 2.470582 
6102 So5070 2.966848 Bea Ae 3S 


oS 





10. 


LIST OF REFERENCES 


Barret, CC. R., Nix, W. R., and Tetelman, A. S., 
Principles of Engineering Materials, pp. 208-2 
Prentice-Hall, 1973. 


Hertzberg, R. W., Deformation and Fracture Mechanics of 


Engineering Materials, 2nd ed., pp. 289-322, John Wiley 


and Sons, 1983. 


American Society for Testing and Materials, ASTM 
Standard E399-83, “Standard Test Method for Plane- 
Strain Fracture Toughness of Metallic Materials.” 


Rice, Js Re, "A Path Independent Integral and the 
Approximate Analysis of Strain Concentration by Notches 


of the American Society of Mechanical Engineers, pp. 


— SS —— p>) et ie eee ce oe mmm ama a A == ee ——- > a a ip ms eS ce —- = = oo oe 


379-386, June 1968. 


American Society for Testing and Materials, ASTM 
Standard E£813-81, “Standard Test Method for JfIe, A 
Measure of Fracture Toughness.” 


Mullican, J. N., Fracture Toughness Degradation in HY- 


tee SS a Se a ma => a = === —Et Se ot SEP 2 = SP Ga = 


80 and HY-100 After Prestraining, Mech. Engr. Thesis, 


Naval Postgraduate School, Monterey, Ca., 1983. 


Defense Research Establishment Pacific, Victoria, B. 
Oar Materials Report S80-A, JIc Measurement Point 


SE —_— aes EN ee a = ap Sa ap os 


R. Matthews and G. D. West, March 1980. 


Brown, WwW. F., Jr., and Strawley, J. E., Eds., “Plane 
Strain Crack Toughness Testing of High Strength 


-_— SES = Se SS = aE SS ap 


Testing and Materials, 1966. 


Bucci, J. R., Paris, P. C., Landes, J. D., and Rice, J. 
R., “J Integral Estimation Procedures,” ASTM STP 5314, 


American Society for Testing and Materials, pp. 4%0-69, 
LO72 « 


Rice, J.R., Paris, P. C., and Merkel, J. G., “Progress 
in Flaw Growth and Fracture Toughness Testing,” ASTM 


STP 336, American Society for Testing and Materials, 
pp. 231-245, 1973. 


118 





ake 


12. 


13. 


14. 


1S. 


16. 


17. 


Ernst, H. A. Paria, PP. C., and Landes, J. OD. 
"“Estimations on J-Integral and Tearing Modulus From a 
Single Specimen Test Record,. ASTM STP 743, American 
Society for Testing and Materials, pp. ¢76-502, 1981. 


Begley, J. A. and Landes, J. D., “The J-Integral asa 


of 
the 1971 National Simposium on Fracture Mechanics, Part 


Il, ASTM STP 514, pp. 1-23, 1972. 
Joyce, J. A. and Gudas, J. P., “Computer Interactive 
JIc Testing of Navy Alloys,” Elastic-Plastic Fracture, 


ASTM STP 668, American Society for Testing and 
Materials, pp. 451-468, 1979. 


Underwood, J. H., “JIce Results and Methods with Bend 
Specimens,’ ASTM STP 677, American Society for Testing 


and Materials, pp.463-473, 1979. 


Clarke, G. A. and Brown, G. M., “Computerized Methods 
for JIc Determination Using Unloading Compliance 


Techniques,“ ASTM STP 710, American Society for Testing 


and Materials, pp. 110-126, 1980. 

Clarke, G. A., “Single-Specimen Tests for JIc 
Determination-Revisited,"™ ASTM STP 743, American 
Society for Testing and Materials, pp. 3553-575, 1981. 


Wu Shang-Xian, “Crack Length Calculation Formula for 
Three-Point Bend Specinens,”“ unpublished, 1983. 


ako 





INITIAL DISTRIBUTION LIST 


No. Copies 


Defense Technical Information Center 2 
Cameron Station 
Alexandria, Virginia 22314 


Library, Code 0142 2 
Naval Postgraduate School 
Monterey, California 93943 


Department Chairman, Code 69 alk 
Department of Mechanical Engineering 

Naval Postgraduate School 

Monterey, California 933943 


Associate Professor K. D. Challenger, Code 693Ch > 
Department of Mechanical Engineering 

Naval Postgraduate School 

Monterey, California 93943 


LCDR W. Kurt Tritchler 2 
c/o Mr. William H. Tritchler 

3737 W. Miller Lane 

Milwaukee, Wi 53208 


Mr. Michael Vassilaros, Code 2814 Pd 
David Taylor Research and Development Center 
Annaplois, Maryland 21402 




















































































































edeeed © Wit Ara 4 wtad nie egh 4 Oto OH Pde Pe CA gk Crt 6 nh 6 in Ak Bes - pets doy Las Meo 
Tek Mek te Oe ke ee oe oe ee eee oe ee va 
BO PETS CEO BEEK © OM Ce COAG Ootid ab Wh CADE Bid Oe MLK se op R ACen * FE 

z 


prakanecsae das stre ere Elastic-plastic fracture toughness testi 





1a Ai FAs OR tet MERE MR OO bbe ge oan. 




































elt hg £2 A oh PMD A ON tee a CLR hd e ne F Cre hp Amp ADE EHH» f 

14D shed on WBirt o® of a Ohh eg if '5 012. dF MED BA J'. 2e'bsteeetsf a teegad tee tt 

CAA Ae rad 6 CAP UEO mm 6 Biel Rat A.OG Bid O68 Rt OMe Cod HeD Hd Labo a ROMEO 8 wa 4 Hy I nt) a CURE veneit ; 

CORR Ob: OO 8 bh PID! ere BA ONE Oh A UE LEE Me A0ehoe @ Wifinh i fon O a b tte ae f A fl I | yy AEE Ret ee 

Ed OO 10 OBOE HE A #0 OLA OM ye ME OE PALL OS Ch HO 4% ¥o08 ff ead : | | | Al! | | iH} | aa | 

A ata OO ORO OE e bon td wen bla OO aAIAT A ALPlLE be 19 re ae Gal HT \ | Wh } 

oh NE hh OME RG Ng 61k O18 deh Bem OF 66 OR mem of oe 2 4 ae a we {| Hi | }] } Wa dl 

ah tO td ht Ah. 6 Bach BO Be) weed tO MIE DIES \G H OMe wth 44 618 Bobs 4.8 6 ap One J ' HY | | } lf if | 

tg eM DMO eSB Mf Ae Os Hots Prk ee ee sae tS Hatt Wt | | | 
whatad 16 ora wd 00% 02 VON SAVAGE A INS Heh Bee SM eted FOE ELA yes Gm ROO Ew o> ss HHA Hil | | AT Hy 
oth het POG Rf) 6 bk T Min ee AAAS Oe eke e FAP CHP CLIT 8 hind Saw HEE Oem 7 aq ie al ' Hell UM ae {| | i | 

Mae Ste ate ORS te Ph OE UF ot git el 4 Ot a OH Rheem + Bt ew ee ee te oe ‘ & et Cn oe | ! a | IH | | 

et Aes A Beats = tinh Pao henn§ #9, r bie Mebae Foe ide revues tam Kagea he ie WP | *: ana «Big was vate bee aia 

Aber omet hws, 29:0 hlO 419 b. wdtrQ 6d Bee Fg th Sew. teins % Soe 4 ry Ar ‘ te 76 

nk. 0 0 wy en ELLEN. 0 Fd, fe 01918 Bem B ieee AAP aah Ape 1G sas Bi mde stom ORAL 8 47 

U Regine £2Apdpe y ta A apt ei 2582 UNC" 6 ting Gomi e AD OC RAIS Vt ge Hinf Oi COEDS tm Cee 6 6 Feolgel Ft be « 








ee eta Rhos ath, OF we Gabi syrh Pod ve PCr Ce eee ee ee ee ee ” Aer 


enkeDI AO Ow taset b's e chiimane ‘ses caend Pete de Ahn tie: of eek 6 i 16 Meo 6 ae DUDLEY KNOX Lf 
Pint AOR Aid 6 eA GA AAOALILE Set othe i Sb i a oe BOO AM Biaiys Ph Pea a Nhe Fy Cle 4 I aR ct RARY 





? 
bet Pe hed 8G shifate ee 8 eh Garde let el eed al Mel oF 2.5 pes © he MUL A tee Cet DOS a Fake” & 












































Ueki Be Arape, moti @ i odiad iat Oh et Me Ae A © Ai he it 298.8 Pw Bok wit he Sty a ie Fd 2 me 

Rae ek Eb ab sere Mh « AER MPV @ sedan. Geng a Af prin Grebe rdgh eh. 2 849 6 A ODAC 8g rete ' 

is pote f ML ORS O81 * why MEAT PeM nd OE CUA CG eM EL MRA es DRA ABD eer Dat om) 

sas ® sla Meat, 0) 6090 Bt ind A a Ftv eh ROME ik pa cdiet om ted Pe gicm Hwa A oe fb oneh Orit Bove 

very tt a) eee Fe ee et oe ce Pe aes COM OX bya le » bE RP aa & tee har pal # wt 

RAY) Bietiot fy! bow S..meyen 6,06) F4 alts DAV Hw ws UG et 2 ew hy een ea ie Oe ene) A ae eed a8 04 ane 

See ae ee eT ee wer rT ee roe oe ereein sy Sh ae mere a ee Sur mee ee ee 

sis re Bnet AM set) ghaoed &. Soh: AS pense mt Retran® Fin Mi O68-O bo ateien tel. 6. OS OR & ph t hod Mims os tube & 

1D Qed Meaeial? etised *, 210 eh TAD A wie Gite tet, bolt Or HA Vee ¢ CT er ad i AaKGs 6 
PTI ee 1 9 does o het. Os fas, Oy Woe en NO he A te a 

iA epdiln'ch Oe bee: « severieetl dF Apite pe nen as yuh Fhe we i Oe fe ADt ee OS 8 homnd re 

OReb et of Ade A wh | LW Be bo bef mM A om © , Wa 66 MME AUTEM 8 BAAR a! 0 HEGRE mgs helt Fuse! 

WF OO ah Got ee oed af bam Od Hi) gripe of of ake AS Monty bare 6. de Fee Ul os e * ae a Me 
eh-GrO A AAA OM Cb. e OS MO MLRM wt- ad uO Se Oke hal tre Pome Si aalne ad apa 6th weg ee Ee i 
ere Co Ph ee te eee eae es 2 ee ee ee fio) J ' ‘ Ate 








lt 2008 1. ymethukar SO set Unt 6 INF mi th Gree titind se Wo vee s Mime ALAR . js 
te Als a Go ed Be ae Porte wit) oe Meir at ee ee ee ee Cy ae CK bh Mad ay 8 athe MA beh 
Med Ag O10 mw, «tt ree 6 bid ahi argh ae ot Mm be 018 Med hEbeied AA. Oodeert og Hdved nH aahihnh 6 F hdak 6 48 Do 

























ADA 0 MOR Mat ee , IPM AEAMIM ell Sue hi e's lh @ “matt 6 led ie 4 ae mf od BE bbe F gel OOPS Kon F 
OeROT AMEE LAAT MS eH MAAR See A Mey FIRB a gthh UA. Ow COW a! oid OOOO ft #2 oie 
pn ee rt ae Ore a eee he eee ee ee Te eee mm, .9 ae Canrdels Mea sa SPO ahd ng 
Mrisiyé mag = Ww eh, Ash ee ALO .65* monn bak OAM cn trteice HO. eh Eh em oy SORA! CIP ahed WEL ot WME esate oy gs om 
nine ee Pein § get riamhaw qe Alt © ws wit Rieti, OOM die. b ain dtl OEE 1 EH yd eg H08 held 
aA OAD pha POE 4D .018y) caw adhsl me | OF 11 See whet aL fea fie iif # aT | 
2b P Kaas tedn ch A EL Cah wo meieed wp wre dln at, hop v0 for 
eh ones oh RM Ret LOM Aah A M00 20 AG eA AI AS Gs FOP 6 at um aw sé 





ali Pat aieh  ardee ht, 00 * AM Pb ume Qprn am Avbcoeees Nottia boo 2 018 6 















Ah Bet a Sen eWs rh Oe hd moO O18 etd ip Byte a Bie ah ew Me Ee ¢ 
W eof dah ee ee ate A Age tonne a 
Ri Meld «abot n diel @ 9h cbt A bid ep A DPA MEM £20 Kn HG tian 
ade Det ot 0 Bahr af Bo eG am AD Aehola A et dud O44 iw i ha Ye $n? 4 » ee 
CaP fs OAM hae OO) Anhed peer Are Te 4h 0 Fark @ OMUH OH Sam» LL a0 cae 
ee ee ee ee eee er ee ee oe aw) Tet) Ore ee 
Ot 008 Debate bo Bad Ad ate ari Ot A AE a me EF PuWehesint sib os ag wwe ween 
Owe ee re ee | ye of A mate ep GR Ge 8 Aceh ott? Coe end -™m* cA" 8 

6 ee NA AA Mies ah P Pree! Jip les S pl ete i 240 MRTG. 9 Ob pws Oe Var ?s aT le ‘ 1% 
ede et want em PA eh BAM AAG Bd ome AO MARSH OF we La A M 
Oe Oe MP eee ee yee er ee oe ee ee er ee ead i ee4s ate 
OP OP Tt Pa ot te Pl ae el oe ee Ver o-d FMEA OUt Ne ote aE me 
Pana @ AEE? CR Per 56 gtbiveod Aa wh & oh SA \s foo AA Ge wd a 
PREP “et Pa ebn AM TEES Ap bad SRE BF re phyS Sepetop es eh ee Oigynt we FAHY eras ee 


r] deat, 6 8 





Safer 8 6 re POU .ahe 4 ah 18 ss Com sh # Ot gt Tht et Ot od Gs 










‘ 
Phong £ as we da eee mtie at Got OC be eia fete =Ases he 0 Mas ot he ‘ 
dk oe te es oe Pe eT ee 2 aacee ' 
Ae Lk Bee iS PMMA BP Ot Fok bm Cams LPs Bin Gt tl apts t 
ee een a eee a ee hn 4 er ed 

watem g Fiat A.g Oop ot — PAMEO ES EB sre aes adaé Ae ff a 





ee Cot be OA Bet MMO See Ft vad bt shee oa ta! COS helt Is 








pitt aeh gg <P Pactweae dle Dp Ober Bead poms hOG a Og 8A ta! 
rad tb bh ei dit ge mtyel BAAD FOP oF tit AE MAA Behe AO 0 depot # 
PrEr AAG? fof A OF mit HG eee nehr? oil eanys 2 aad s10 argos o 
Cots cote! Aare ie ‘De Loe PR At A gh e } ee unk. 
g -@ 8A Oem o teecee oe Veg ar oe ye fee EC aera - 
al Pry | We eee Se vo we ap eden ches Aen dD 4 . Ter 
¥ a0. Soe a ata LARA PPM A 19 VOR OPS! Sf tee * eel UK, 
SMe Arised hide te) 6, Vg 48 eee At nt hee ae ot hat ef s ‘ ‘ & 





a Pat eGasaqs 94 8h vio! I Mates & 4 at gee * si As Rat ow a¢ 
yet 00 eter paentin dh "et tp PHOUISPA he Cran d aah gd itt A sfleg 
Fic gs oe Weswe ee VW gmetd FF al cp heres, De bw FO BN tw 6 ootem ‘ 















: 
ce ee ee eee Pe ee ee oe 4" aaeitinw rah hS 
¢ eet db Ae AOL Ram bf et DEX 8 Oo MeOH # oh ag nt Balt oh A aw eae ot a 
Popetttewr Sieve gbdewehe 1 6 Vt shame atk te Acs met es of 
eile oh gual od 6 6 eke Fd. wine Om oh te 0 200 bad 2d ait wet A “He 
oe nee be hihi wed eae, iow acs! woer tf eh Fee wd. > 86 Uyts ee 
Oat adh @& Fab ot pyhintte De shy MS Avi aon fe 39 @ + Om Bat ’ ‘ 
te tgif «Dts wy ee ee soe ee ee rel “es re itt) ,¢ 
PE whi Palin fy Ome 4 ar gtiah Mh af a el fn O18 Org MF FY es 
me Ta ml iste ae ee eee ee staid Jv 4 
awn? Ww Ae ay Jig’ t 1 hab net tit Aa eo! #2 
et Aol ge Ao th OM wegl nha dee oo tsa a ‘ : 




















alley ht nae Hs eu “Fae OF APU Ng Ge Mae re ere) Ce ne 
SD gaan twit SFP Alvret ol # . At © Aus 
eat © 8g, Wath - ayo a yt Jrae ? Has 
AdMIUAY Bouse 8 F ry 
aed fas Pes fos wot a: 
au feb ela a" ents vi 8 tet Nes 












te Peet BM. Kes se mate «6 Gh Anveet "£ « =) i 
ee et a i ee oy SY oe a et Fie eat oreg 
i Bante Fe 6h ttn Fo! £ ro oe 10 Sot i 
Pads aero ge oe Se. er ye ‘a yg 
elhas ten Se" laa ae wen tie Lees 
: 
‘ 





ab Wh Sat 
ats rad 






, a , 
aj we a, ee ”, 




















Pmgse wT = $s HH Y, 
es ws > wt hms Soh 
’ oe e » WE AAH oe 
4 = . ~~ tay + Gee Woy 1%) % -~ = . 
Big tt my cay VV www ey mo ‘emits M « =e len,» 
PASA AW Ye 4 ERE TH sale & 
Ere pode, ee SR we ee dilbine® = Set neg Ne he ae 
mm & ky Rg OKT] Oey of ~ ae wn ere . - coh 
& Tutt ot he A ph a ee 
wt” noe a vr NN. ‘= we eo eet Qari & 
% a. te ee vor OP or Aas ‘ 
ry wtwh . i Bh ea Sule “A GG nite ey 
me etal § BAe. BY? ye He, a | 
Se ateeP Oe Rh EQh © Cm ye SQ SHA QE Ne 
oy SAGE eS at Oh ie 
Miw & VSQrareRR™ AAak ue 1 & e« ee 
ee te eS ~~ of s ~ 
~ % wm, AUT TTY ree rm © he | « Sgae 
me bw tm te Py Se & A My let toh eter 9, “ss 
hey tee ayers, He om ivre Ketek 
Leh = me MD % ye ae ®& 














































































































































py Se oe ~~ r= poe 0 St Sa, 
fe Nee eral wR yee 26 mere | PY ys e™ “~~ 
re =o ~~ 1a, > ~e ~~ =e YS - oft iv 
ou od te el Be Ue ee ~» eee & Gfme 
. > Pe were OF OK Te AD bw =H Yt nS 
Pepe URN Se GS acer Bowe yh Ee aa v% : * : an . : a sk ‘ { " 
= fe Meee i hy we or Se 7 mu mye ta tig 4 . ag % 3 ’ * ‘ ‘ 7) % in *® » ' a Lae a 
; S * Fergie : + ‘deo <* ein ve a a ee a! ea "es 4) ’ ’ ‘ lie 1 « Te . . \ Heal DOG as 
. BWR AS we a, wee, 8 84 yee A, WHA ' F eS sat rk F . a ' ‘ ' ' ‘ t 
oe UG Sam Se wees syle eo y* iw Da ans a a ew oP PUA w As, a 3 8 a he. ae 1 fs, a iar . ane ’ 1 
COE qedy ere” He We hie e ‘se? ‘ ve ease t t te va ‘ '" Se ar" ie . sae +4 a bh oN 1% -* ' 
Wry evite? - me ee Ui efews bE Dk VE Gg RTE! el VOM AEE *" wg Re wie - vi Cae oe . a he | e “hw a3 a6 « ‘ aa & e 0 ¥ a ts. fe t 1¢ . ‘ ‘ 
PD ee) re mig OOF SH WH] 21a’ WES “Syerpa™ 4 Y y + er " vs “* an -- ee $ vi 4 Ani} a at * 
uty £8, “Sore AAW tw oe ® ag twee y Hs PSR « 4 ex VVC YW... % Ct ee » ASME CIA | 6 nae ‘ 4 ' he ea oa? a 
RFE ee Sy FA mele hk FG Ceo rh ® a Me ve) Bs A, euUieu sy “rey te en tO ys ee he) \ o. | m+ .* " er ‘ ‘ . 
ee et ee ee a oe i ce or eh ete Fe ao ~» PF) Teeny ¥ * Se a ee yaw 2 De © awee¥ . e ' ie | d 7 a> , cy 2 . 5 
9 Sy We EG WD ge Stern sw ale hh wu a i ee | é eet Rp SH Ce Uy yr i. ‘ *e Pukaea , Cir on . + ry . a Ve esi § ’ £ 4 id 
rh, Se arged A ee ot eee a a ewe The ok q's a 1.aQy asa te ae ER ’ u oe =e) ye eS on s Pt a e % ' et ‘ . 
wx o%, aa apres = a ed eke a a oe . Ve 'O'S | ae Ts » 4 = ar a ar wi 1 1, 6 F > vie s t ‘ 
VR OES Rh eS DEH ttt ale ™ eR ER ee we ie, Henesg wt. we'y a at ta Soot FU at -* % a . © Hous « a fe 7 
. WAP! SFE HiQnhOry * 4% ee wt * > ya * ah we, &Ty ‘oa ¥ * ' oy wee . + cag 44 t ve ’ ‘ 
Ut Ak wy art, ee bre @ Se eS Fea Ge VER re Ow A nk to, Od) ra 7 FAS F Bey ee ae P » 26 a %,¢ » € F ' 7 
We De EEE ym ewig we Pe SEGA G seh ny ea Ai ey ot Wipe Wife i es -h& «8 ws ' ‘ : ‘ ' : ae 
Py FE SR ee A SA’ & OE PA aries no ie ah m seve. B® “—" v Nw eé a ty i vt ° 2 “ose 8 , 3 ' at Ns ' » , 
ade 0 tele Ue 2 eg hf SD TRS OPER ate Sth te 1m S te He Beelmye we Rse? OY yi, wry % - ae P v, i ¥ 7 ‘ . a6 
wah eee . Reese eas Vig Pee R’UVUEN we S*EE ce AP Y teu tm etek eye - Be 7 f a ae 4a F ey 1s . 1] ' « ‘ « 4 a 
Pe TRAE EM see YN RMT SN we ee Ow EO WE) Ww, a ~ay o4an. CAR vraea , t ake oe sate, a ‘ ‘ 
SePe ol ye Hare Dwele ay Aen. > OLR OH Tt EPR A ENS * BEd Pee] me Ae gid CN Lar x ’ “Wa Yon, ad A nt eo oa ° al gy . ie " 
i Hare, Ty Re lee rays te Hee LOM Qt. 6 UPR AD Bq: oh WN yh G6 ee i Le We id On J 7 . ave a 2 y ' 4 1 52 ‘ ‘ Uy 
1 My Dp aE ey res | SS ORS HHS E CEs She CUBE WI ae oD Wy oer oe <ep, oy x 9 ri er en i | 1 4 as Da Oe ' a : : 
i wh Fas ay ay a, pe Om 8 Oe OS. SE OEM EPL 8 HR, ON BA Gly we Se Serv BO BRE" 7) ie . ‘ e a y 4 
nr or toe th en Fa Te ee 2 oe ev Y ¥h'e ¢,'9 et ¢hhy "i ee ? a os 1 ‘ Oe | e ‘ . ae 
eT a ee ee OL Rin ehh he tT me, eh ree She ae tee eters t Ne oa be ' Wars 
ach th ke aay! HE EE'S SAMY TRH WUT HA UH VAN tod VULN FOUTS 8, ate es % . ¢ t * ae 
Re BAe TOD Ee ep oye | Abe 4 A STOPS Bee Owe Pils wa By > kva 8 & . + vets ® ae 1? ‘ ‘ “ &« U 
PE lg te Se ep Bee tere ig ett 6A, CEG REE, he UR ou SPARE ye BT vag, hehe Ba Gite ‘ 1% i aa a | 
eS BSB Ok CVC eM lyre ie wy IER OTR ENS. | AACE ET OT EMM CHW Meg Ue KR EER VA! at i's ’ al tae) a 
ey tested) © Ee PW Ce, ear SO EY Gt as TA R'E iG ee FH HR ABS att ' :* » “te aa x ' “ eo , 
ED BSN Veaws Br QA ler, ew Gin Deemed EY EQ ae oe oe oe i ee ey vn 6 S's % rf tec af © ug? ‘ a be \ ' 
Dedhiedin Letelen, ton, fo Fe, t Sch Pie tehath “ete ton 2 Bes Th dees the hh Ae dh nd Le de ee ee Cr re oh i ee is é ear a 2 a | ' GC. t - ' ‘68 a ae ‘ ee 
TS OEY Ee yk BIER MP @ POS BAS re ree Eu Ge, VA ,PRREVED EWE Ta, . e . U.P boa] D ‘f° J y «© Ee ae a 
at one YAY Sree ate oe 'Y UK ED UDR TRAE Oe dU LESAN EN ARAGON OLY f . as 8 a Wet er 2 . * +t 5 1 Fee e a 
PEWS A. Cg Ute Are oN UCR CW WEE oe e ® % SFk EE WET Foe be Ow, VIR A SE Oy Ft 2° gov We te a™uia % ‘ 24 asad . : a 
Pr ae lg he ee SAL SY Va mew args ag ey' & a vAwERS, ae, eh Ae ghee we Gr, ay rue ¥ oe 7 ‘ ° x ‘ * ‘ ba 
Qay Ay AE Ah Asst ew @ UN Ue wae pe'wead ae a oC ee, oh #46 i | UPPED T SS hE ’ FS « ts f a ' 
ee tee te A ee ee bode ae Se ee ere ee cs ker a | eka tewEal. CYS We eure tw ' La a ' i eat * ‘ ' 
AR A, BO1O S| Be WH RA BD SOG 02 Re, 4 WTO el tryre ate CHE Y Vi He OD HUE a te Uw Rw 8 ° . wa a oY . ‘ a! t “ e ' 
pre tap eee ee he HEY we Pa ere SG Ob A TE CWE Se Gey Ate TEAL ME ew Nee Et ase- es t sh Y, 08 oat + 4 ¢ 
fre Te TEE aL, B/E Mew. WHATS Fay YS WEVA SS P FU REM CUBS A VFGAA Es WUE Uae DR, . 6.4% ~ . : ¥ Ht, ava 6 . Jes . : 
| Sh Sg Ha Ube Ee Ds kee ok We aH ut! & 6 ose hike EV EHO YES RH RM Bye ew eR, Me GS zs « ar i i ’ 2 4 
rok oe eee eh he ee et Re. Pe ee ee here he te ene ee Vevige S 601% qa wy erties Ff NM cM, “ ri Cee i ar ' s e 
eae th ede eee me ok oe a 2c ENS a ee ee ee ee a Te iuyw 6 ¥ ‘ " " ia , ee « .? ‘ ' 
UA te, & ATT Rey ee ae ON EQ QE, © OR | BP SD OREM EWA DD Ko uis, VV GUO BE hr OT u, Te etanr. te p e 6 tC me «4 ‘ ‘ ' : 
OX sieves “Sa -s RR SNS! € wy! 22-8 Os * Le a Oe ee ie @ Riegel ee eWHA wee, an ee y Le ee | ‘ mn wrdsewm wa i * 6 he i thie ‘ . ‘ 
Se WR S PIAL EG Clean et-g on TRE, 8 ROR bee ye WOR HOEY, | oo ‘ ie a Se Ti en a F ' ius P Vo” ' : 
Lae Ls th oe ee ee de ee De Oe or si AG ss oh eee 8M 1 We oo 6 t.¢ a ‘ ‘ aes 
er sD ay & +h mer be eh Bh OFS TE & ew RAHI Se VIVT VHeyvwin OG watawe © RG Ban Ss, rt a. oo We a t Le iy wef ry wea . ' aa ‘ ' 
PRN ae ay ee he eS ee SF BERD VGH EULA NOM oS VOM, 6, OST YD, ny » PR eval es awe we oak Dr FY wt he wr 1 8 * va ‘ 
sy Oe Sy ey gg i BTR ow Gere UAE Cb EV eH Awh we OR ‘a CCR RE HVA ES © ee sae v “ ott . ’ z ‘ a ove ' . 
100 Te Coreg a ye ie CO he BS A Fe BESO) SUD bere. Gee g © 0 4 » ee Ve WIRD FOE ge eG ee 6 " ‘ ve ei ’ Was ‘ ' ‘ n ¢ ‘ ot 
baptist halos Latebda 1. okt EAR AB Uh OO HUT: EY i Ae te a eh Se ek i Ys? , a 8% a 8 ' ' i ‘ ' 
Be ee Sie Rigi. Bey oh MWA ATA CRD GUE Owe RS TWh, TOMY PIL VY WEY PEP PWRE ea, br roe ‘ 7 4 i ‘ er pe ‘ re ' .. ' 
ie 8d Py emp rig Ae ee ee ee ht Le ee | wee rd 6 EVD UK YY VOGT. PASH SAE 0724 Ss OW i r ~ ¢# i v 1% 4 ' , ’ 
ere Cee Qe © tO Cem re 4) eo ere iy 4h @ RE HY AES? © Om UR Wik weyst “y Be 4 ee . rast ThA a a *aq & ie a bos ' a 
Me OER 8 Ye May hy a ie BPR Py em) ORE 6 @ 7G ETB ag He). DP | 8. Kore eave vi emtat ee a « % t ‘ a ¢ ‘ 
On ily Ut Wh EOE SS 8 OS OR OUTS BO AY PT DDH Be OU EL he el WY & ye HT Se Ps vi euy ’ x a) . 8 ' 
Tee OE EG ww) Fy ety ye Tye a Be, ere EW RE Eh HAY DE MGrw TNE UT RAY OD Bes | we a =a ' : J oe ‘ 
re 4 OE ES Rie TS ESO OORT A © UTD & htt POS Ewe Yr dy & & Fw tH 6 O° US Ue RO we ¥ ‘ ' © te t i « . 
ae ASUS Oy 9 a wiles Ter ME OH AW Ven a we e@'4 tinea Uy © yh HAS > S-wW { Wrede @ . ‘ Bay, ba % ‘ P, yy Ce s or. % 8 i | 
VO | BARU way ag Ag eG Ry BA SOP Fin Ae > De 6 wR Py VPM Hoy hs HUY |e RW Sy: ©% *¢ >Gek ' ‘ ' | ‘ e ' ‘ : 
KAU Ecce ra SAD AMY SED AS Be 4 BE | DO Qe AN el ¢ orton FD, oe oe ee yO : y ey vo ' ° r 
NCR ote ey OOM we OG BLT Reb 6 ee ROU CUTE WA WHE DaAe GB SP eres Rhee Ce a CTD vs .* ‘ ’ 4 ae? is 
Ng Ee AAW, OH ERE eS OATH SEE RRY PY PMT SH WHAT HRAKne” J Et-5 SEP KEWY*AHY 9 v * e be . , * 
we ae Rew O'S Re ath ely Ate WMAPE PMG ioryy  : , dm 2a Es we, OA mh e'¥ be, ee weBOSE wwt.¢ YW, ‘ + t v ' a ’ 
i te ke ee ee a i be Oe Oe | Rene +h BER" ciety oy: VR WAM AE POA. wey Ate We Qe OD > & ne F at ‘ ' 1 1¥t . . ’ 
Leg RW BOGS EA ORC ET § rh: MH yy 6 ON Wee. SOU CHF Beh ey +. = 9% ' ' as ws ¥ % bd ‘ 
eT Lt he ee CO lh eee a he Ge en oer ee a ¢ a ty, a “ 
Be a eT te ie ee eo ee ee 2 ee ee « Cr a a we vt a at 4 ‘ 
St ee OQ wert ese Din tree PUTS Y Se, 2 ? . ‘33 ' ' 
Pe RE BS 6G FON TES RAE SY TD FB » 4 +, # 14 » * » ¢ 
Ain hp ee ee SE ” Seley 7 Ow ee HD eS eey 5 ‘ t ’ J 
WS eee we Sere" A 6 OG ES PVH MPAA VE 88 or BET YW C4 e. ‘che mn ’ ' ' ‘ ‘ 
te et SO St BRE OR CYP RY SE Oy OWED hae Ge 3 " y ‘ . ‘ . t : a 
We Seat teed Bes 2:8 VHRR GP Gwe 34 QW eQiQ Ey Pome. rly gb © . e ant a ' 2 ¥ * 4, t 
SS Ee, tn Me Wii Sy A 4 On Be DAT. FPR IL OY Te YO ERY Ay é % “eh . ose 
ht ee haa 7723S +S a, BA « : + «2 a y ® ‘ ‘ ' ‘ “ 
ee -ed el we AY B42. 1F- S10-8 30. e >y a oe ee z a ® a hi ‘ 








‘ 
’ 
‘ 
‘ 
: 
. 0 . 
‘ 
to 
reat Tal 
' 
' ‘ 
' t 
e 
' 
1 ' 
oo ‘ 
’ 
' ‘ 
’ ' 4 
‘ 
‘ 
’ ‘ . 
' 
‘ a t ‘ 
fee ' 
’ ’ 
' r ' oo 
t ‘ “ ' ' 
ae ' : 
' 
' -d ‘ 
8 f 
: t= a ' 
° 
' ‘ ' ’ 
20 offre + ‘ ' ? $ 
‘ 1 6 
“poo ' 
ao 
a ‘ . 
f 
‘ ' 
' . ’ 
‘ ’ 
’ ‘ ‘ - ® 
' 
’ ' 
‘ . a 
. 8 F] 
‘ ’ ' . 
soe 
Pa) 
» 8 er 
» 4 
'e 7 
' ' . 
ni . 
Cee : . ' 
Ue ee wee ‘ s ae 
’ a 8 ‘ 
8 , 8 ‘ ' 6 
' 
« ' : . 
' 
, ‘ e : 
ae ’ oo ' 
‘ oa ‘ 
' se ‘ se 
. 4 ' ’ 
’ oo ‘ 
; ‘ 
‘ toa 
: ' a) 
’ 
' ’ poo 
: e 
mm. : * 
; 1 
5 ' 
. 8 ' 
CAG 
' 
fer 
4 ' e . ‘ 
“8 . 
‘ 1 
on a 
1 
e P so a ‘ 
' tae 
’ ' oo 
fos 7 
' a 
‘ a 
' ' ' 
. 
' " 
on ' ' 
1 . . 
: ‘ 
1 ‘ 
‘ ' 
' ‘ . 
a 
a 
» ¢ 44 ‘4 
; a 
. . 
‘ . 
‘ 
v ‘ ' 
a ‘ 
6 . . 
’ 
' 
. . 
a 
1? fos 
at 
. . a 
’ ‘ ee | 
‘ 
’ 
a 
' . . 
a : 
a 
' a 
“ 
' 
a a 
’ 
lc . ¢ 
« 
‘ 
' 
' 
. 
a 
‘ a 
a 
‘ 
a 
' ' ' 
' 
| ’ 
' 
. 
a 
. 
' 
a4 
' e 
a 
‘ ' . 
‘ . 4 
4 . 
a 
a 
a 
» oo 
® . 
4 
’ « ‘ 
' 
. 
' 
a | 
' 
4 
e 
. ' ' 
' . 
‘ 
e . 
: . 
4 
' 
‘4 
' 
‘ aa 
a 
2 
. a 
at t 
a 
‘ 
a ‘ 
‘ 
' ' 7 
Fy 
ie . 
a 
' 
a 
ny 
a a 


